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FOREWORD 


This document updates the Shuttle digital flight control design 
developed by the JSC Guidance and Control Systems Branch and the Draper 
Laboratory Control and Flight Dynamics Division. The flight control 
algorithms of the Space Shuttle from launch through touchdown, and for 
powered horizontal flight are presented. It is anticipated that this docu- 
ment will be utilized in establishing an initial baseline for the digital 
flight control software design. This report supports the generation (by 
Rockwell) of the Level I "Orbiter Avionics Functional Subsystem Software 
Requirements" documentation series. 

Revision A of the Digital Flight Control Software Design document 
incorporates updates from the design activities performed over the past 
six months. In addition, Honeywell has completed an extensive review of 
the original version and many changes in Revision A resulted from the 
constructive review comments to NASA. Additionally, treatment of control 
modes in this revision reflect technical design exchanges with Rockwell 
and the Flight Crew personnel. 

There are several important design areas which have not been covered 
in the revision. It is anticipated that they will be addressed as base- 
line additions in the future. These areas include the treatment of redun- 
dant sensors and redundant actuators , interface with pilot type displays , 
and roll-out control. All design to date has assumed rigid body models; 
however, the filter structure has been specified to facilitate the handling 
of bending and slosh dynamics. 

In this document, a number of design elements have been addressed in 
detail. The following paragraphs summarize the scope of the document. 

The unified structure which is presented for the DFCS (digital flight 
control system) is a flexible framework for the evolving Space Shuttle con- 
trol system. 
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Two designs for jet selection are included. The first, which is 
based upon linear programming, is considered flexible. The second, which 
is based upon a table look-up scheme, is more restricted. The latter is 
recommended for the baseline design until the Shuttle requirements are 
found to go beyond its capabilities. 

The design for the cruise phase addresses the requirement for a 
variety of manual, semi-automatic and automatic control modes during the 
developmental horizontal flight tests as well as during the final stages 
of operational vertical flights and during ferry flights. The control 
gains have been calculated for the Rockwell OO89B vehicle, but only for a 
single flight condition; the scheduling of gains for a full range of con- 
ditions is currently in progress. 

A number of features have been added to the design for the entry 
phase. The applicable flight regime has been extended to subsonic, thus 
eliminating the separate transition phase, and a state estimator has been 
introduced. Rudder augmentation has been added to the control during the 
early part of entry. The gains have been calculated for the Rockwell 
OO89B vehicle. A capability for higher maneuver rates in response to 
large bank commands has been added. The elevator- trim algorithm during 
early entry has been changed from open- loop to closed-loop. 

The design for control using the reaction control jets during the 
on-orbit mission phase has been developed and implemented, as has been the 
design for the attitude maneuver routine during the phase. Attitude 
control design presented herein does not reflect the availability of 
vernier jets or the unavailability of a clock for on-off commands other 
than once per minor cycle; these factors must be considered in designing 
firing policies. The design for control during major maneuvers in orbit 
(including insertion) is in development; only a functional outline is 
presented in this document for that segment of the D1US. 
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There are areas where the various design efforts have not yet been 
fully integrated. As an example, the phase plane developed for on-orbit 
AC PS control may well be usable in the entry phase . 

A version of the DFCS which incorporates the designs for the cruise, 
entry, and boost phases has been implemented in FORTRAN for testing on the 
SSFS (Space Shuttle functional simulator). A second version has been pro- 
grammed in HAL using structured programming techniques. It incorporates 
the designs for the on-orbit phase, and runs on the functional simulator 
at the Charles Stark Draper Laboratory. 



ABSTRACT 


The objective of the integrated Digital Flight Control System is to provide 
rotational and translational control of the Space Shuttle orbiter in all phases of 
flight: from launch ascent through orbit to entry and touchdown, and during pow- 
ered horizontal flights. The program provides a versatile control system structure 
while maintaining uniform communications with other programs, sensors, and con- 
trol effectors by using an executive routine /functional subroutine format. The 
program reads all external variables at a single point, copies them into its 
dedicated storage, and then calls the required subroutines in the proper sequence. 
As a result, the flight control program is largely independent of other programs 
in the GN&C computer complex and is equally insensitive to the characteristics 
of the processor configuration. 

Section 1 of Part I describes the integrated structure of the control system 
and the DFCS executive routine which embodies that structure. Section 2 
deals with the input and output. The succeeding sections show the specific 
estimation and control algorithms used in the various mission phases. 

Part II describes attitude maneuver routines that interface with the DFCS. 
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SECTION 1 


INTEGRATED STRUCTURE; EXECUTIVE ROUTINE 


1. 1 Integrated Structure 

The DFCS program consists of an executive routine, several major (functional) 
subroutines, and a number of utility subroutines. The executive routine performs 
" housekeeping' 1 which is required regardless of the mission phase and calls the major 
subroutines in the appropriate order. There is a major subroutine for each of the 
following functions: Backup Control System (BCS) implementation, "pad load" manip- 
ulation, input manipulation, filter and parameter initialization, filter update (or meas- 
urement incorporation), control logic, command transmission, filter pushdown, (or 
state extrapolation), and parameter estimation. Most of these routines are divided, 
at least in part, into separate segments for the various mission phases. 

This structure of executive routine and functional subroutines gives modularity 
to the program; it allows complete flexibility in designing the state estimation and 
control algorithms for each of the mission phases, designating new sensors or con- 
trol effectors, and indeed, redefining the mission phases themselves. At the same 
time, however, the unified structure which groups together the coding for similar 
functions encourages economy both in computer storage, and in engineering design 
and programming effort because it is so easy for different mission phases to use 
common coding whenever it is appropriate to do so. Another advantage to this struc- 
ture is that the calls to the functional subroutines indicate the completion of one set 
of calculations and the beginning of another; hence, they are suitable break points 
for external interrupts, and they facilitate the testing of new coding. For example, 
if modification to the orbit thrust vector control law is contemplated while other seg- 
ments of this phase's computations are unchanged, the new coding can be added as a 
totally separate subroutine, leaving the old coding intact. During testing, it is then 
possible to switch between old and new control laws following completion of the first 
segment of state filtering. This feature also could be of utmost importance for 
emergency or backup coding that is required to be functionally separated from the 
primary coding. 

The unified DFCS structure has additional, managerial advantages that are 
valuable in an enterprise as complex as the Space Shuttle program. The use in dif- 
ferent mission phases of the same assumptions about input /output, timing and other 
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interfaces is encouraged if not enforced. Furthermore, commonality in conventions, 
nomenclature and formats is encouraged, which simplifies and clarifies both testing 
and documentation. 

Not all operations of a control system need to be performed with the same 
frequency. The current version makes allowance for 3 sampling intervals, the long- 
er intervals being multiples of the shortest interval. In the current cruising flight 
phase, for example, the "fast” rate (10 sps) is used for angular control, the "medium” 
rate (2 sps) is used for velocity control, and the "slow" rate (1/2 sps) is used for 
altitude control and parameter updating. The medium and slow computations are 
offset, in order to prevent them from occurring on the same control cycle. Allow- 
ance is made for different sets of sampling intervals in different mission phases. 

1.2 DFCS Executive Routine 

An overview of the logic flow through the DFCS Executive Routine is shown in 
Fig. 1-1; Fig. 1-2 illustrates the implementation of this logic. The major subrou- 
tines of the DFCS are denoted in these figures by the hexagonal blocks. Although the 
executive serves all the mission phases, many of the algorithms within the subrou- 
tines are specifically designed for a single phase. Table 1-1 defines the mission 
phases for which the DFCS makes provision. The inputs to the Executive Routine 
are listed in Table 1-2; the outputs are listed in Table 1-3, The operations performed 
by the Executive Routine are discussed in the following paragraphs. 

The phase indication is changed from "entry" to "cruise" when the Mach number 
drops below 0. 9. All other determination of phase is currently assumed to be done 
outside the DFCS. 

When the Backup Control System (BCS) is requested, almost all the control 
system functions are performed using algorithms in a special subroutine that is sep- 
arate from the rest of the DFCS. The executive transfers control to this subroutine 
if the appropriate flags have been set. The BCS mode has been designed only for the 
cruise phase at this time. 

If there is an external restart request, or the phase or mode has changed since 
the previous pass, the executive sets flags to 

1) call for a DFCS initialization pass 

2) call for initialization of parameters for any modes that are new 

3) cause the computations that are done at the slower sample rates to be 
performed this pass. 

If the reason for initialization was an external restart request or a phase change, 
the Executive Routine establishes the counters and timers that reflect the desired 
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Fig. 1-1 Logic Flow of DFCS Executive 




rom External Scheduler 



Fig. 1-2 


Detailed Logic 


Flow of DFCS Executive 
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Table 1- 1 DFCS Mission Phases 


1 

2 

3 

4 


5 

6 

7 

8 


Rollout Transition from aerodynamic control to steering and braking 

with landing gear; after touchdown; not currently implemented. 

Cruise Approach and landing following entry phase (currently below 

Mach 0. 9) and horizontal powered flight (take-off through land- 
ing). 


Entry 


On- orb it, 
TVC 

On- orb it, 
AC PS 


From after the deorbit maneuver to cruise phase; from total 
reliance on jets (ACPS), high angle of attack and hypersonic 
to primary control with the aerodynamic surfaces (ACS), low 
angle of attack and subsonic. 

Thrust vector control during maneuvers with the Orbital 
Maneuvering System. 

Control with the jets during coasting flight. 


Insertion 


Boost 


Insertion into orbit following the Boost phase; not currently 
implemented, probably identical to On-orbit TVC, 

Flight with solid Rocket Motors and/or Main Propulsion 
System. 
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Table 1-2 Inputs to DFCS Executive Routine 


ISTRTR 


External request for DFCS initialization (1); none (0); 
cleared in Executive Routine 


m 


Mach number 


MOD A, MODE, 
MODR, MODI! 


Current control mode - aileron, elevator, rudder, 
speed loops 


MODAG, MODEG, Current guidance mode - aileron, elevator, rudder, 
MODRG, MODUG speed loops 


MODAGR, MODEGR,] Requested guidance mode - aileron, elevator, rudder, 
MODRGR, MODUGR speed loops 


MODAR, MODER, 
MODRR, MODIIR 


Requested control mode - aileron, elevator, rudder, 
speed loops 


PHASER 


Flight phase flag (external) 


Table 1-3 Outputs of DFCS Executive Routine 


TM, IS 


IS TART 
ITM, ITS 
MODB 


Medium- and slow- sampling- rate indices: 0 = perform 

calculations this pass 

Initialization pass (1), not initialization (0) 

Number of fast passes per medium (slow) pass 
BCS mode (1), not BCS mode (0) 


MODACH, MODECH, Indicate mode change - aileron, elevator, rudder, 
MODRCTT, MOD UC T I speed loops: 1 = no change; 2 = guidance mode change; 

3 = control mode change; 4 = phase initialization 

PHASE Flight phase flag (internal copy of PHASER) 

t , T,,, T c Fast-, medium- and slow- sampling- rate sample 

F * M S 

periods for current phase 
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sampling rates for the current phase, and calls the major subroutine which transfers 
the set of pad- loaded constants required for the current phase from the pad- load lo- 
cations into working storage. The "pad- load” denotes those parameter values that 
are determined prior to a flight and are loaded into the flight computer on the launch 
pad . 

Regardless of whether or not it is an initialization pass, the Executive Routine 
calls the Input Interface Routine, which reads all the dynamic input required by the 
current phase and mode: sensor measurements, navigation estimates, guidance 
commands, and manual commands. All the dynamic data is read at the same point 
to ensure that all the values are registered with respect to the other values. The 
parameters are differenced, rotated into other coordinate systems and otherwise 
manipulated so as to be in the proper form for use in the other routines. During 
the appropriate phases, an Air Data Routine is called which manipulates some aero- 
dynamic parameters. These two routines are described in Section 2. 

At this point, the logic flow is detoured in an initialization pass to go through 
the Filter and Parameter Initialization Routine and the Parameter Estimation Routine. 
The former assigns values to those variables that need to be initialized when the 
phase (or mode) is entered. The latter determines the current values of parameters 
that will be regularly updated in the course of this phase (or mode). 

The flow then returns to the normal path. The subroutines entitled Filter - 
Part 1, Control Algorithms, Control Output and Filter - Part 2 are called in sequence. 
These subroutines comprise the core of the control system. After they are executed, 
the Parameter Estimation Routine is called; it is not called at this point in an initial- 
ization pass because it was executed earlier in the pass. The state filters are divided 
into two parts in order to minimize the amount of computation between data input and 
control command. The first part is limited to updating the filter estimate according 
to the most current measurement. The second part, which occurs after the control 
is exercised, is used to pushdown data (in the case of classical digital filters) or to 
propagate the state (in the case of discrete- time modern control estimation). The 
Control Algorithm Routine scales and combines the error signals from the filter rou- 
tine to form control commands for the aerodynamic surfaces, reaction control thrust- 
ers, etc. The Control Output Routine, which is described in Subsection 2. 2, pre- 
pares the control commands for the DFCS interface with the control actuators. The 
Parameter Estimation Routine updates the gains required for estimation and control. 
The functions performed in this routine may include scheduled or adaptive gain adjust- 
ment. Much of the Parameter Estimation Routine is executed only at the slow sam- 
pling rate. 

After these routines have been executed, the Executive Routine increments the 
indices that indicate the "medium-" and "slow- sampl ing- rate" DFCS passes. If it is 
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the end of an initialization pass, the "slow- sampling - rate" index is given an extra 
increment. This extra increment will smooth the load on the computer by offsetting 
future "slow" passes from "medium" passes, provided that one sampling rate is an 
integer multiple of the other. The Executive Routine then returns control to the ex- 
ternal executive. 

The algorithms incorporated in the routines pertaining to input and output are 
presented in the following section of this document. The contents of the other rou- 
tines are presented phase by phase in the succeeding sections. 
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SECTION 2 


DFCS INPUT AND OUTPUT 


This section describes the interface between the DFCS and other components of 
the Space Shuttle avionics. Subsection 2. 1 discusses the DFCS inputs. It lists all the 
inputs; i. e. , state measurements, guidance commands, and pilot commands. It in- 
cludes a description of a routine which processes some of these quantities prior to 
their use in the filtering and control algorithms. It also describes a routine that fills 
the function of an air data computer. 

Subsection 2. 2 lists all the outputs from the DFCS to the control actuators and 
describes the limiting performed on these commands. Finally, Subsections 2. 3 and 
2.4 present two candidate jet selection schemes; the first is powerful and flexible, but 
the second may suffice for the Space Shuttle requirements. 

Much of what is presented in this section depends heavily on many features of 
the Shuttle avionics that have not yet been specified, much less modelled in a simu- 
lator. Therefore, in many cases convenient assumptions have been made for the 
sake of proceeding with the DFCS design. These assumptions include the following: 

1) The inputs indicated in Subsection 2. 1 are, in fact, available to the DFCS. 

2) The actuators indicated in Subsection 2.2 are, in fact, available to the DFCS. 

3) The DFCS does not deal with redundancy management of sensors or actuators. 

4) The DFCS does not create outputs for pilot displays or telemetry. 

5) The DFCS can schedule software tasks to turn ACPS jets off and a 
sufficiently fine clock granularity is available for such scheduling. 

It is expected that each of these assumptions will, in time, be replaced by more 
realistic specifications, at which time these designs will have to be revised. 

2. 1 Inputs Required by the DFCS 

Table 2- 1 lists all inputs currently required by the DFCS in the various mission 
phases. They are presumed to be available from sensors (or sensor- interface modules), 
the guidance and navigation module, and the pilot (or command- interface modules). 

In this interface, and throughout the DFCS, English units (viz feet, pounds, and seconds) 
and radians are used. 
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Table 2- 1 Inputs to the DFCS 


Phase 

Mnemonic 

Symbol in Text 

Definition 

2,4 

ACC BOD 

a b 

Linear accol measurement, (body coord) 

2,4 

A G 

a G 

Angle-of-at^ack guidance command 

8, 6 

AT V (1 ) 

0 

Platform ou^er qimbal angle 

8,6 

AI V (2) 

e 

Platform inner qimbal angle 


A I V ( 3 ) 


Platform middle gimhal angle 

8 

A I V C ( 1 ) 

0 c 

Platform outer gimhal angle command 

8 

A I VC (2) 

i^C 

Platform inner qimbal angle command 

8 

AT VC (3) 


Platform middle gimhal anqle command 

2 

A Y y 

ay 

fl ea sured lateral acceleration 

2 

AZ y 

az 

Measured normal acceleration (tzbody) 

4 

PAN KG 

0 VG 

Bank guidance command 

2 

8G 


Sideslip guidance command 

2 

CG 

C G 

CAS guidance command { QG or CS T A RG ) 

2,4 

CRTS 

C b 

r 

Ref er e nce-t o- body direction cosine matrix 

2 

DAG 

6a G 

Aileron guidance command 

2 

DDG 

6 G 

T A EM guidance interface variable 

2 

DEG 

6e G 

Elevator guidance command 

8 

DELOPP 

6 orb 

Prestored orbiter engine deflection command 

8 

D E I.S F M 

6 srm 

Prestored SR M engine deflection command 

2 

DKG 

6r G 

Rudder guidance command 

2 

DSG 

68 G 

Speed brake guidance command 

2 

DTG 

6t G 

Throttle guidance command 

2,4 

DYNP 

q 

Dynamic pressure estimate 

2 

HG 

h g 

Altitude guidance command 

2 

Hr 

H 

Measured altitude 

2 

TAGAIN 


AGAIN switch position (1 thru 4) 

2 

IATRTM 


Roll stick *rim activity: 0 , * 1 (rgt) , - 1 ( 1 efi 
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Table 2- 1 Input to the DECS (Cont) 


Phase Mnemonic Symbol in Text 


2 rip; AIN 

2 IKTP T r 

2 ir.EAiy 

INGUIP 

2 T P H * 

2 IPCHAI 

2 I PC hi: I 

2 I PC AIN 

2 I^TPIM 

T S T RT ? 
2 HODAGF 

2 * O D A F 

2 ^ 0 D H F 

2 HODKGR 

2 M OD v P 

2 or PGP 
2 lOCPP 

2 10DHGP 

2 'lopur 

2 PG 

PHASFP 


INGUID 


ISTRTR 


P G 

PHASER 


2 PH SEL 

P*1 


SEL 


Definition 


EGA IN switch position (1 thru 4 ) 

Pitch stick trim activity: 0,«-1(fwd) r -i(aft1 

Extend gear (1), retract q^ar (0) 

Flaq specifies guidance command set 
PHfi valid (1) , invalid (0) 

Request change I A I N T state (1) , no request (01 
Request change IETNT state (1), no request (0) 
RGAIN switch position (1 thru 4) 

Pedal trim activity: P,*1(lft),-1(rgt) 

Request ini t i al izat ion ( 1 ) , no request (0) 


Requested 

va lues 

for 

I10DAG 

Reg nested 

values 

f or 

'I OD A 

Requested 

va 1 ues 

for 

*ODB 

Requested 

va lue s 

for 

10DEG 

Reg uest ed 

va 1 ues 

for 

10DE 

Requested 

va lues 

for 

'lODRG 

Requested 

va 1 ues 

for 

1 0 DR 

Requested 

va 1 ues 

for 

*0DUG 

Requested 

values 

for 

^ 0 DU 

Roll rate 

guidance command 


Flight phase flaq (e xt erna I ) : 1 =rol lou t phase 
2=cruising phase 
4=entry phase 
5=orbit T VC phase 
6=orhit PCS phase 
7=insertion TVC phase 
P=booster TVC phase 
Pi lot- sel ected values for PHRFF 
Measured roll rate 
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Table 2- 1 Input to the DFCS (Cont) 


Phase 

Mnemonic 

Symbol in Text 

2 

PSG 


2 

PS M 

$ 

2 

P S S 2 1 

^SEL 

2 

QG 

q G 

2,4 

() M 

q 

2,4 

PM 

r 

2,4 

ROLL 



ROLLG 

0 bG 

4 


— rel 

2,4 

T E MPF K 

tk a 

2 

THC 

®G 

2 

THM 

e 

2 

T MS EL 

®SEL 

6,P 

T NOW 

t 

n 

2 

UG 

U G 

2 

US EL 

U SEL 

2,4 

UX V, 11 Y V , II Z V 

UV , uv / 
‘ x “ y 

4 

VRODY 


4 

VMACH1 


4 

V? 

v _ 

—rel 

8 

WR (1,1) 

RRATE. 

i 

8 

WB (2,1) 

PRATE. 

i 

8 

WP (3,1) 

YRATE. 

i 

2 

X AM 

xa 

2 

XF M 

xe 

2 

XFM 

xt 

2 

XRM 

£r 


Definition 

Healing guidance commani 

Measured heading 

Pilot-select ed values for PSPFF 

Pitch rate guidance command 

Measured pitch rate 

Measured yaw rate 

Computed roll angle 

Guidance commanded roll angle 

Position vector in reference coordinates 

Input, air temperature in deg, Kelvin 

Theta guidance command 

Measured theta 

Pi lot -se lect ed values for THPEF 
Present time 

Velocity guidance command 

P i lot- se lect ed values for URFF 

Unit vectors along the velocity coordinate 

system, in reference coordinates 

Air-relative velocity in body coordinates 

Lateral switch parameter 
Air- relative ve loci tv in reference 
coord 1 n itns 

Individual roll rate gyro signal 
Individual pitch rate gyro signal 
Individual yaw rate gyro signal 
Measured roll stick deflection 
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Measured pitch stick deflection 
Measured flap controller deflection 
Measured pedal deflection 



Table 2- 1 Input to the DFCS (Cont) 


Phase 

Mnemonic 

Symbol in Text 

Definition 

2 

XS M 

o«/ 

XS 

Measured speed brake controller deflection 

2 

XT* 

xt 

Measured throttle controller deflection 

9 

YACIN (I) 

Y. 

L 

Individual y accelerometer signal 

8 

ZACIN (I) 

Individual 7 accelerometer signal 
Table 2-2 Inputs to the Input Interface Routine 

Phase 

Mnemonic 

Symbol in Text 
~b 

Definition 

2,4 

ACCBPD 

a 

Linear acc^l measurement (body coord) 

2,4 

AG 


A n g le- of-a t t a ck guidance command 

4 

BANK DO (2) 


Bank rate filter extrapolated value 

4 

BANKG 

0 VG 

Bank quidanco command 

4 

CA T 

C 09 Qf rp 

Cos (AMEAN) 

2,4 

C PT3 

c b 

r 

Reference- t o- body direction cosine matrix 

2,4 

DYNP 

q - 

Dynamic pressure estimate 

2 

M M 

H 

Measured altitude 

2,4 

PM 

P 

Measured roll rate 

2 

OG 

q G 

Pitch rate guidance command 

2,4 

QM 

n 

Measured pitch rate 

4 

FcDOTC 

; c 

Bank rate command 

2,4 

RM 

r 

Measured yaw rate 

2,4 

POLL 


Computed roll angle 

2 

ROLLCt 

^bG 

Guidance commanded roll angle 

4 

FR 

R , 
— rel 

Position vector in reference coordinates 

4 

SAT 

s inofrp 

Sin ( A M E A N ) 

4 

VP 

—rel 

Air-relative velocity in reference coordinate 
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2.1.1 Input Interface Routine 


The purpose of the Input Interface Routine is to preprocess measurements and 
compute state- estimator inputs. Inputs to the routine are listed in Table 2-2. Out- 
puts are listed in Table 2-3. Interface- routine operations are governed by the dis- 
cretes INGUID, ISTART, PHASE, and MODA in the manner summarized in Table 2-4. 
The remainder of this subsection is a description of the logic flow of the Input Inter- 
face Routine, which is illustrated in Fig. 2-1. 

The routine is called by the DFCS executive. If the phase is not either entry 
or cruise, it returns immediately. If it is a DFCS initialization pass, some variables 
internal to this routine are set to nominal values. 

The section of the routine for preprocessing of measurements is then entered. 

The angle of attack, the sideslip angle, and the air speed are computed using the 
air-relative velocity and the direction-cosine matrix. The Air Data Routine is called 
if the phase is in the aerodynamic regime. Some aerodynamic "measurements ,f are 
generated in the Air Data Routine, which is described in the next subsection. The 
computation of bank angle (under the indicated conditions) concludes the measurement 
preprocessing segment of the interface routine. 

The next section of the routine computes some nominal state values. Nominal 
values are computed for those state components of the entry phase state estimator 
(Section 3) for which nominal values are not specified by guidance. Nominal body rates 
are computed to provide turn coordination in the entry phase. A nominal (commanded) 
bank rate is also computed for use in the entry phase control logic in the specified 
flight regime (Ij ^ 0, q > q 3 ) when the bank error is sufficiently large. 

The filter input computation section is now entered. When in cruise phase, 
either the current roll angle or the bank angle measurement is put into the cell PHI 
for processing by the cruise phase filters. When in entry phase, the state estima- 
tor inputs (namely, errors in the angles of attack, sideslip and bank, and errors 
in the body rates) are computed. Additionally, estimator extrapolated values are 
adjusted for changes in the state- variable nominal values. 

Finally, a check is made to ensure that INGUID is negative. This avoids re- 
processing non-updated G&N inputs in subsequent DFCS passes. Control is then re- 
turned to the DFCS executive routine. 

2.1.2 Air Data Routine 

The Air Data Routine is called by the Input Interface Routine during the atmos- 
pheric flight phases. It computes several parameters required for aerodynamic 
control, serving the purpose of an air data computer. 
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Table 2-3; Outputs from Input Interface Routine 


Phase 

Mnemonic 

Symbol in Text 

4 

AT f 1 ) 

L a 

2 

A M 

a 

4 

AO (2) 


4 

BANK! (1) 


4 

BT {1 ) 

A? 

2 

n M 

& 

2 

HM 

H 

2 

PHS 


2 

PH M 

0 

4 

PT (1) 

Ap 

4 

01 (1) 

A*T 

4 

PI O) 

A r 

2 

f J* 

tr 


Definition 

Anqle-of-attack filler input. 

Measured angle-of -a ttack 

Ang le- of-a t ta ck filter extrapolated 

Bank filter input 

Sideslip angle filter input 

Measured sideslip 

Measured altitude 

Bank angle guidance command 

Measured bank angle 

Boll rate filter input 

Pitch rate filter input 

Yaw rate filter input 

Measured velocity 


value 
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Table 2-4 Discretes Governing Logic Flow in Input Interface Routine 


ISTART 

DFCS initialization pass flag. 

= 1 : initialize variables this pass. 
= 0 : do not initialize this pass. 

PHASE 

Mission phase . 

= 2 : cruise phase. 

= 4: entry phase. 

(Routine not entered in other phases. ) 

INGUID 

Flag from guidance. 

= ± 2 : guidance commands are pitch rate and 
roll (about body x-axis); terminal part of cruise 
phase. 

= ± 1 : guidance commands are angle -of- attack 
and bank (about velocity vector); entry and early 
cruise. 

< 0 : as above, but no new guidance commands to 
process this pass (set negative internally). 

MO DA 

Mode flag for cruise phase aileron channels. 

= 4: automatic mode; compute roll (or bank) 
filter inputs. 
t 4: not automatic mode. 

J LAT 

Mode flag for entry phase lateral control. 

= 0 : early entry, compute desired bank rate. 
= 1 : late entry 


Note: The parameters q 3 , r^j, r ( j > 2’ r <^ 3 are used in the Input Interface 
Routine as limits. They are introduced in the discussion of the entry phase 
control logic in Section 5. 
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The inputs to this 

routine are the following measurements which 

from outside the DFCS 

or are computed by the Input Interface Routine: 

tk a - 

temperature in degrees Kelvin 

a - 

angle- of- attack (radians) 

H 

barometric altitude (ft) 

U 

air-relative velocity (ft/sec) 

The outputs are: 

O' >-p 

mean (trim) angle- of- attack (radians) 

sina T | 


sine (cosine) of a 

cosarp ; 

U M * 

mean airspeed (ft/sec) 

m - 

Mach number 


The means (which are the outputs of first-order filters with time constants 
equal to the slow sampling period) and the sine and cosine are updated every pass. 

The means are set equal to the values of their respective measured inputs in an 
initialization pass. 

The other outputs are only computed on slow- sample- rate passes. The com- 
putations are shown in Fig. 2-2. 

It is expected that, in the future, when more realistic models of the sensors 
are available, this routine will compute the dynamic pressure and perform other 
computations . 

2.2 Outputs to Control Effectors 

The DFCS computes commands for the control effectors listed in Table 
2-5. The question of how the commands are transmitted to these effectors is con- 
sidered to be outside the scope of this document; the following simple implementation 
is currently assumed: 

1) Commands are transmitted by the DFCS to the actuators (or a software 
module that interacts with the actuators) at the end of each execution of 
the control routine. 

2) The DFCS limits the magnitude of each output, but not its rate of change. 

3) Timing of jet firings is obtained by means of a computer clock that can 

be set to call a DFCS routine which will turn jets off after the required time. 
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Fig. 2-2 Computations in Air Data Routine 


2 - 








Table 2-5 Control Effeelors Used by the DISC'S 





Phases in 
Which Psed ' 

Aerodynamic Control System (ACS) 

elevens, rudder, speed 
flaps 

brake. 

1,2,4, .8 

Attitude Control Propulsion System (ACPS) 

40 jet thrusters 


4 , 5 , 6,7 

Air Breathing Engine System (ABES) 

throttle 


1 , 2 (horizontal 
Flight only) 

Solid Rocket Motors (SRM) 

2 engines; gimballed at 
body axes 

45 from 

8 

Main Propulsion System (MPS) 

3 engines; gimballed 


7, 8 

Orbital Maneuvering System (OMS) 

2 engines; gimballed 


5,7 


sjt 

Phase designations defined in Table 1-1. 


Table 2-6: DFCS Limits for Commands to Control Effectors 
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The assumed sign conventions for the aerodynamic control surfaces and the associated 
pilot controller motions are illustrated in Fig. 2-d. 


All the commands are limited within the DFOS prior to being issued. 1 he 
limits magnitudes are shown in Table 2-fi. The aileron and elevator commands are 
treated separately -ven though they apply to the same pair of actuators, the elevens, 
and are, therefore, interrelated; this question will be dealt with in the future. 


The limiting and blending (for the elevens and M PS gimbals) of the actuator 
commands is performed in the Command Output Routine, which is called by the DFCS 
executive after the Control Routine. The Command Output Routine serves as the- 
single point from which the actuators (or the software modules that interact with the 



5a - (5EL - 5ER)/2 



Fi'g. 2-3 Sign Conventions for Aerodynamic Control Surfaces 
and Manual Controllers 
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2. 2 ACPS Jet Selection Subroutine 


This section describes a scheme for converting angular -rate-change and 
translation requests into timed jet firing policies for the ACPS (Attitude Control 
Propulsion System). This routine has been tested with the OFCS program. 

2 . 2 . 1 Thruster Geometry 

The thruster configuration currently incorporated in the jet selection routine 
is the North American Rockwell ATP design. 

As shown in Fig. 2-4, the thrusters for the ATP configuration are grouped 
in 16 locations. There are 40 thrusters in all, and these are distributed among 
the 16 locations as shown in Tables 2-7 and 2-8. There are 16 thrusters in the 
nose of the fuselage and 24 in the "tail". 

2.3.2 Jet Selection Logic 

The DFCS Control Routine determines a commanded velocity change and calls 
the jet selection subroutine, which handles angular velocity components and trans- 
lational velocity components simultaneously with a single algorithm. A flag can be 
set instructing the jet selection subroutine to ignore translation, simplifying and 
speeding up the jet selection procedure. Furthermore, the jet selection logic can 
be instructed to deal with jet clusters or with the individual thrusters. 

Figure 2-5 is a functional flow diagram of the jet selection logic. A dis- 
cussion of this logic follows. 

The ACPS jet control authority subroutine is only fully executed when the 
jet selection routine is called for the first time in a mission phase. The control 
authority subroutine calculates the ACPS control authorities on the basis of total 
vehicle mass, the vehicle inertia matrix, and the center of mass location. The 
control authorities are calculated for the jet clusters and for the individual thrust- 
ers. The control authority is the angular and translational acceleration in vehicle 
coordinates generated by firing a jet or cluster of jets. In addition, a cost per 
unit firing time is assigned to each jet. Currently, all these costs are equal. 

When a jet is disabled (or restored) and the jet- selection logic is operating on 
clusters, the jet control authority subroutine is partially executed, appropriately 
altering the corresponding acceleration vector and cost element. 

The jet selection subroutine next checks its inputs to see whether or not 
there is a new velocity change request. If not, the jet selection subroutine call 
must be for the purpose of turning off jets which have completed their firing 
intervals. If the velocity change request is a new one, the linear programming 
subroutine is called, which expresses the velocity change vector as a linear 
combination of jet acceleration vectors, the linear coefficients being the firing times. 
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Fig. 2-4 Thruster Locations - NR ATP Baseline Orbiter 


Table 2-7 Distribution of A CDS Thrusters 
Nit ATP i baseline Orb iter- 


Thruster Location 

Numb-T of Thrusters 

Thruster 1 1J) 

1 

2 

10, 12 

2 

2 

14, 16 

3 

2 

13, 15 

4 

2 

9, 11 

5 

2 

6 , 8 

6 

2 

5, 7 

7 

4 

17, 19, 21, 23 

8 

4 

18. 20, 22, 24 

9 

3 

32, 34, 36 

10 

3 

26, 28, 30 

11 

3 

25. 27, 29 

12 

3 

3 1 , 3 3 . 3 5 

13 

2 

2, 4 

14 

2 

1, 3 

15 

2 

38, 4 0 

16 

2 

3 7, 3 9 


40 total 
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Table 2-8 Parameters of Individual Jets 


ID Jet Location — Body Station (inches) Thrust Components 

Fabrication Frame Vehicle Frame 


1 

318. 

, 0 

- 96. 0 

374. 0 

-1. 

,0 

0, 

, 0 

0. 

0 

2 

318. 

0 

96. 0 

374. 0 

-1, 

0 

0. 

0 

0. 

0 

3 

318. 

0 

- 96. 0 

349. 0 

-1, 

, 0 

0. 

0 

0. 

0 

4 

316. 

0 

96. 0 

349. 0 

- 1 . 

0 

0. 

0 

0. 

0 

5 

318. 

0 

-109. 0 

368. 0 

0. 

, 0 

1 . 

0 

0. 

0 

6 

318. 

0 

109. 0 

368. 0 

0. 

0 

-1. 

0 

0. 

0 

7 

318. 

0 

-109. 0 

355. 0 

0, 

0 

1. 

0 

0. 

0 

8 

318. 

0 

109. 0 

355. 0 

0. 

0 

- 1 . 

0 

0. 

0 

9 

318. 

0 

-109. 0 

346. 0 

0. 

0 

0. 

0 

-L 

0 

10 

318. 

0 

109. 0 

346. 0 

0. 

0 

0. 

0 

- 1. 

0 

11 

330. 

0 

-109. 0 

346. 0 

0, 

,0 

0. 

0 

- 1 . 

0 

12 

330. 

0 

109. 0 

346. 0 

0. 

0 

0. 

0 

- 1 . 

0 

13 

318. 

0 

- 80. 0 

367. 0 

0. 

0 

0. 

0 

1. 

0 

14 

318. 

0 

80. 0 

367. 0 

0. 

0 

0. 

0 

1 . 

0 

15 

330. 

0 

- 80.0 

367. 0 

0. 

0 

0. 

0 

1 . 

0 

10 

330. 

0 

80. 0 

367. 0 

0. 

0 

0. 

0 

1. 

0 

17 

1536. 

5 

-158. 0 

425. 0 

0. 

0 

1. 

0 

0. 

0 

18 

1536. 

5 

158. 0 

425. 0 

0. 

0 

- 1 . 

0 

0. 

0 

19 

1532. 

5 

-158. 0 

43 8. 0 

0. 

0 

1 . 

0 

0. 

0 

20 

1532. 

5 

158. 0 

438. 0 

0. 

0 

-1. 

0 

0. 

0 

21 

1528. 

5 

- 158. 0 

450. 0 

0. 

0 

1. 

0 

0. 

0 

22 

1528. 

5 

158. 0 

450. 0 

0. 

0 

- 1 . 

0 

0. 

0 

23 

1525. 

5 

-158. 0 

464. 0 

0. 

0 

1 . 

0 

0. 

0 

24 

1525. 

5 

158. 0 

464. 0 

0. 

0 

-1. 

0 

0. 

0 

25 

1541. 

5 

- 130. 0 

480. 0 

-0. 

234 

0. 

0 

0. 

980 

26 

1541. 

5 

130. 0 

480. 0 

-0. 

234 

0. 

0 

0. 

980 

27 

1555. 

5 

-130. 0 

484. 0 

-0. 

234 

0. 

0 

0. 

980 

28 

1555. 

5 

130. 0 

484. 0 

-0. 

234 

0. 

0 

0. 

980 

29 

1567. 

5 

-130. 0 

487. 0 

-0. 

234 

0. 

0 

0. 

980 

30 

1567. 

5 

130. 0 

487. 0 

-0. 

234 

0. 

0 

0. 

980 

31 

1546. 

5 

-130. 0 

416. 0 

0. 

407 

0. 

0 

-0. 

915 

32 

1546. 

5 

130. 0 

416. 0 

0. 

407 

0. 

0 

-0. 

915 

33 

1559. 

5 

-130. 0 

420. 0 

0. 

407 

0. 

0 

- 0. 

915 

34 

1559. 

5 

130. 0 

420. 0 

0. 

407 

0. 

0 

-0. 

915 

35 

1572. 

5 

-130. 0 

428. 0 

0. 

407 

0. 

0 

-0, 

915 

36 

1572. 

5 

130. 0 

428. 0 

0. 

407 

0. 

0 

-0. 

915 

37 

1579. 

5 

-149. 9 

468. 0 

0. 

980 

0. 

0 

0. 

234 

38 

1579. 

5 

149. 0 

468. 0 

0. 

980 

0. 

0 

0. 

234 

39 

1579. 

5 

-149. 0 

455. 0 

0. 

980 

0. 

0 

0. 

234 

40 

1579. 

5 

149.0 

455. 0 

0. 

980 

0. 

0 

0. 

234 

c.g. 

1096. 

5 

0. 1 

367. 9 
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Yes 



Fig. 2-5 ACPS Jet Selection Functional Flow Diagram 
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Generally, many different jet combinations could satisfy a given velocity change 
request; if there is more than one such combination, the total cost is minimized 
in the selection process. Any jet which is flagged because of failure or mission- 
related constraints is omitted from consideration. 

The selected jets are turned on (if they are not already firing), and all others 
are turned off and removed from the internal jet- off list. The interrupt clock is 
set to return control to the jet selection subroutine when the next jet- off is due, 
and the "new velocity request” flag is zeroed. 

The interface parameters for the jet selection package are: 

INPUT: 

a) Velocity change command vector. 

b) Flag for new velocity vector. 

c) ACPS jet exclusion list. 

d) Translation inclusion flag. 

e) Phase flag. 

OUTPUT: 

a) Jet- on^off commands. 

b) Interrupt clock settings. 

c) Zeroed flag for new request. 

Figure 2-6 a presents features of the linear programming subroutine, CCWUJS, 
and Fig. 2-6 b describes its logic flow. 
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PURPOSE : 

INPUTS : LROW - 

LCOL - 
W - 

COST - 

B - 

LALARM - 

JF AIL - 
BI, CF, LQ - 
PARAMETERS - 


Fig. 2-6a 


Linear Programming Subroutine 

Represents input vector W as a linear combination of linearly 
independent columns from matrix B, choosing columns from B 
to minimize cost of solution. 

Number of components in W. (Dimension of requests) 

Number of columns in B matrix. (Number of jets or clusters) 

Vector (size LROW) to be analyzed as a linear combination of 
columns from B. 

Array (size LCOL) of strictly positive cost coefficients, one 
associated with each B column. 

Matrix (LROW by LCOL). Each column represents the accelera- 
tion associated with a one-second firing of the corresponding 
jet in vehicle coordinates. Components may include rotational 
or translational acceleration elements, or both. 

Non-zero indicates that previous jet selection may not be taken 
as a starting basis, for one of the following reasons: 

1) Jet exclusion list has been changed 

2) B, LROW, or LCOL have been altered since the preceding 
pass 

3) Preceding delta- V request was not successfully resolved into 
jet firings. Initialized to 1 on first pass. 

Array (size LCOL) of flags, J^ 1 element non-zero to indicate 
jet failed. 

Can be inputs as well as outputs. See definitions below under 
OUTPUTS. 

SAVING, MAXPAS, MAXITR, TOLRNC, URGNCY can be chosen 
to match the precision of the computation and the accuracy of the 
data . 


Features of Linear Programming Subroutine, CCWLJS 
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linear Programming Subroutine (Cont) 

OUTPUTS : BI - Transformation matrix (LROW by TROW) going from vehicle 

coordinates to coordinates based on selected columns of B. 

V - Solution vector (size LROW)_of coefficients for selected columns 
of B matrix, satisfying: KI W = V 

C F - Vector (size LROW) of cost coefficients selected from COST 
array, corresponding to chosen columns from B matrix. 

Z - Cost of solution, satisfying: Z = CF • V. 

LQ - Array (size LROW) giving the indices and sequence of selected 
B columns. 

M COUNT - Number of basis element substitutions performed in current 
execution . 

LA LA RIM - Non-zero if W was not successfully resolved into jet firings before 
MCOUNT exceeded MAXITR . 


ASSUMPTIONS: 

D 

Cost array entries are all strictly positive, i.e., greater 
than zero. 


2) 

Negative IXJLIST outputs and negative V components are to 1 
excluded from any implementation of the solution. 


3) 

COST array entries are small compared to URGNCY. 


4) 

BI matrix is never altered by any other program. 


5) 

LALARM will be input as non- zero whenever changes occur 
in LROW, LCOL, B, or the jet exclusion list. 

ENTRY: 

Call 

CCWLJS (B, W) 

EXIT: 

RETURN to caller. 


Fig. 2-6 a Features of Linear Programming Subroutine, CCWLJS (Cont) 
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For I 


For I 



1 to LROW, then go to 840 



(840y 

1 to LCOL, then go to 100 



MOOUNT counts iterations on this pass 


If BI cumulative roundoff is not bad, 
if previous solution was valid, and 
if LROW, LCOL are unchanged, 
then start with solution from previous 
pass. Otherwise initialize completely. 

If W is unchanged from previous pass, 
use same solution. 


NPASS counts iterations since last 
complete initialization. 


Negative LQ values identify vehicle 
frame coordinates, which are non- jets. 


Positive LRCORD values correspond 
to non- selected jets. 


Fig. 2-6b Flow Chart of Linear Programming Subroutine, CCWLJS 
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BI W - V 


V contains latest version of solution 
in terms of current basis. 


I = 1 to LROW, then go to 110 


1/5(1) > CT 
and 

sy(I)^0. . 


1 V ( T )|^0 


Cost target vector is_initially a 
weighted version of V (or a signum 
(V) vector). 


loop to 
v 101 v 


URGNCY*V(I) 

-CF(I) 


URGNCY factor places premium 
on jets which approximate W or 
which fire in directions represented 
so far by non- jets. 


URGNCY - CF(I) 


:vd) = 


URGNCY-s ignum(V(I)) 
- CF(I) 


loop to 
V 101 > 


MCOUNT > MAXIT1 


GO TO 
. 445 v 


Y-CF BI 


SAVING - II 
1 + LROW - J 
1+LCOL -» K 


Y (I) is weighted cost of generating 
the I th original (vehicle) basis 
vector as a linear combination 
of current basis vectors. 

Initialize the search for the next 
jet to be included in basis. 


GO TO 
, 200 > 


Fig 2-6b Flow Chart of Linear Programming Subroutine, CCWLJS (Cont) 








For I - 1 to LCOL, then go to 300 



Search for next jet to be 
introduced to basis, bypassing 
any which are failed or already 
in basis. 


X is the weighted price penalty 
for simulating a one-sec firing 
of I th jet, using the basis rather 
than the jet itself. 


Greatest saving per second is 
H; corresponding jet is K. 



If no jet was selected for 
inclusion, then exit. 


K jet is to be included. Y 
is K>h jet acceleration vector 
described in current basis. 


Fig. 2-6b Flow Chart of Linear Programming Subroutine, CCWLJS (Cont) 







For I 



= 1 to FROW, then go to 330 



Choose .T basis member for 
replacement so that the new 
V(J) is ^ 0. , and no other V 
component changes sign. 


Fig. 2-6b Flow Chart of Linear Programming Subroutine, CCWLJS (Cont) 
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Bad result. 


Fig. 2-6b Flow Chart of Linear Programming Subroutine, CCWLJS (Cont) 
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EXIT 



Calculate _cost of solution: 
Z = CF • V 

Save W for comparison 
in next pass. 


Negative firing times represent 
requests which cannot be met 
with currently non-failed jets. 


NPASS reflects quality of BI 
matrix. 


Return to caller. 


Fig. 2- 6b Flow Chart of Linear Programming Subroutine, CCWLJS (Cont) 
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2.4 AC PS Jet Selection Logic for the ATP Conf iguration 


A jet selection routine is presented in this section that is more restricted 
than the design of the preceding section. It is essentially a table look-up; 
the firing policy for each possible command, as well as the individual jets 
and the control authority corresponding to each jet cluster are pre-determ ined 
and stored for access by this on-line routine. The North American Rockwell 
ATP design for the ACPS, which has been presented in Fig. 2-4 and Tables 
2-7 and 2-8 was assumed in the off-line computations. 

The flow chart of Fig. 2-7a describes in detail the jet selection logic 
for the ATP configuration. Figure 2- 7b defines the flowchart symbols. A 
discussion of some of the features of the logic follows. 

General - To provide attitude control with no translational disturbances, jet 
couples must be used. However, pure jet couples exist in only one axis, the 
roll axis. Jet couples to be used for pitch control (e. g. 3 & 12) and yaw 
control (e. g. 5 & 7), also effect roll. To achieve pure pitch or yaw two 
thruster couples must be fired to null the roll disturbance; for example, 3 and 1 
12 with 2 and 10 for pitch and 5 and 7 with 3 and 1 for yaw. Hence, couple 
control can require as many as four thrusters to fire whereas non- couple 
control may require only one thruster firing (e. g. yaw control). Additionally, 
the multiple thruster firings required for couple control produces large minimum 
impulse rate changes and these cause fast limit cycles and inefficient propellent 
usage. This logic has the option of attitude control with either coupled or 
non- coupled jets. 

Nose versus Tail Thrusters - for Non- coupled Jet Operation - The logic allows 
the crew to select either nose or tail thrusters for control of each axis. 

This serves two purposes: 

1) Comparing nose thrusters to tail thrusters, the nose thrusters have 
shorter roll moment arms, but longer pitch and yaw moment arms. Hence, 
the crew can select the preferable moment arm length for each axis. 

2) This is a simple means of accounting for thruster failures. Should 

all the thrusters in any single location fail, control can still be maintained 
(Note: Individual failures in any jet location are accounted for by placing another 

jet of the same location in ''common” with the location, see page 16 of the 
flowcharts ). 

Phase Plane Jet Selection Logic Interaction - (Non- cou pled Operation) - 
Thrusters in locations 1 through 4, and 9 through 12, produce both pitch and 
roll torques. Hence, it takes two thrusters to provide either pure roll or pure 
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pitch. However, smaller minimum impulse rate changes are possible by firing 
only one thruster. Therefore, even if there is a phase plane command from' 
only one axis (roll or pitch), the logic will attempt to do the latter (fire only 
one thruster) if the resulting d isturbance in the "non-command" axis reduces 
the rate error (called RR^ or RR^ 7 in the logic). 

Use of the tail thrusters for yaw control causes a roll disturbance; should 
the roll disturbance exceed a rate limit, the logic will consider the disturbance 
part of the phase plane command and cause corrective roll thrusters to fire. 

Translations - No combination of thrusters will produce Y or Z translations 
with negligible rotational disturbances. Thrusters must be modulated to null 
rotational disturbances. This modulation could be performed either open 
loop in the jet selection logic, or closed loop through the phase plane logic 
(the crossing of a switching line) and jet selection logic. The latter method 
was chosen because it will not produce unnecessary modulation (thruster 
firings) during small Y or Z translations (i. e. , when rotational disturbance 
is not long enough to cause switching line crossing). 

During + X- translations , it is possible to control yaw by modulating the + X- 
thrusters; this option is available in the logic. 

Jet Coupling and Logic Flow - With respect to thruster control actions, there 
are two basic types of thrusters: 

1) Those that produce Y translation and yaw. 

2) Those that produce Z translation and roll and pitch. 

Sequentially, the logic first calculates commands for Y translation and yaw and 
then issues commands for Z translation, roll, and pitch control. 
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Commands for translation, either positive, none, or 
negative, depending on +1, 0, or -1 

Commanded rate changes from switching logic. 

Angular Acceleration for i ‘ jet for either roll, 
pitch, or yaw. 

Value of 1 (not 0) signifies request for long moment 
arm thrusters 

Indices for jets 

Minimum impulse time 

Roll rate disturbance for yaw r, nd/or Y translation firings 
Roll rate limit before RDFY is added directly to 0 C 
Roll rate error from phase plane logic when 0 ,, is zero 
Index indicating presence (l) or absence (o) of roll command 
Sample time 
Jet indices 

Maximum roll compensation time available by unbalancing 
pitch firing times 

Small number to account for round off in comparisions 
(say 10-4) 

One jet roll correction time 
Net roll carmand, i.e., RR + 0^ 

Roll rate limit before RR^is added to 0^ during translation 
Sample time plus E 

Pitch rate error from phase plane logic when 0^ is zero 
Couple option 


Fig. 2- 7b Definition of Symbols for Jet Selection Logic 
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SECTION 3 


STATE ESTIMATOR (ENTRY PHASE) 


3 . 1 Introduction 

The function of the state estimator is to obtain estimates of controller state 
variables in an uncertain environment. This uncertainty arises from noise in the 
measurement process and noise driving the states directly, e.g., due to wind gusts. 

The state estimator design must meet certain prescribed requirements in the 
face of externally-imposed constraints. These requirements, many of which are 
not yet specified quantitatively, pertain to the degree of accuracy and smoothness 
in the estimated state variables, the time delay between reading the sensors and 
issuing computed control commands, and the effective phase lag associated with 
filter dynamics. External constraints are imposed by fixed sensor hardware and 
configuration characteristics and computer timing and sizing limitations. 

Meeting all requirements while satisfying all constraints may necessitate 
design compromises. The design should be comprehensive to ensure prescribed 
accuracy by making best use of available information, but design sophistication 
is limited by computer constraints. Noisy inputs may require heavy attenuation to 
ensure prescribed smoothness, but the low gains associated with heavy attenuation 
must not cause excessive phase lag or excessive sensitivity to modeling errors. 
These compromises are central in the development of the design presented below. 

The current state estimator is designed specifically for the entry flight phase 
(and some of the parameters used by the estimator are defined in the entry control 
description, Section 5). The structure used, however, is applicable to any flight 
phase. Extension of the state estimator to other flight phases is under study. 

3 . 2 General Description 

The design presented for the state estimator is linear and therefore implies 
approximation, namely linearization of the state dynamics about a set of operating 
points and linear modeling of control effector and sensor operation. The linear 
assumptions facilitate both design and analysis. Nonlinear effects will be noted as 
development continues and possibly included. 
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The design is cast in the state-space formulation and uses more than one 
measurement to produce each state-variable estimate (multi-input, multi-output). 

It is comprehensive in the sense that it uses all pertinent information. However, 
approximations are made to decrease computer requirements; namely, only a 
reduced number of state-variables are estimated, and decoupling of longitudinal 
and lateral dynamics is assumed. An important feature of the design is that control 
effects are fed into the estimator’s state-extrapolation step, keeping the estimator 
response fast even while noisy inputs are being heavily attenuated. 

Time delay between reading sensors and issuing control commands is mini- 
mized in the sense that it provides for incorporating the current measurements and 
extrapolating the state 1 in two separate steps with the latter performed after the 
commands are generated. Furthermore, the slow-varying filter parameters are 
updated separately at a slower repetition rate. Details of these computer timing 
considerations are given in Implementation, Subsection 3.3. 

The estimator gains are constant with respect to time at any particular flight 
condition; however, their on-board computation is dependent on the current estimates 
of dynamic pressure and Mach number. When the dynamic pressure is below a 
prescribed level, ft the estimator gains are set equal to predetermined fixed values. 
When dynamic pressure is above this level, the gains are interpolated from pre- 
determined tabulated values with Mach number as the sole independent variable. 

The fixed gain values used at very low dynamic pressure were determined 
classically on the basis of desired filter bandwidth and damping ratio . tft The 

jjt ^ 

scheduled gain values used for high (above or q ) dynamic pressure were found 
by analogy to the steady-state Kalman (Wiener) filter for the reduced, decoupled, 
linearized state assuming locally stationary noise statistics at each operating point. 
Measurement noise strengths were chosen to reflect relative sensor quality; driving 
noise strengths were chosen to keep the gains at a selected level, sufficiently low 
to achieve desired smoothness but sufficiently high to exclude extreme sensitivity 
to modeling errors. 


* Analogous to "filter update" and "filter pushdown" in classical filters. 

ft The prescribed dynamic pressure level is^ for the ^longitudinal channel and q^ 
for the lateral channel. Current values for q* and q| are 20 and 30 lb/ft^, 
respectively. 


ttt 


Currently, the low-dynamic-pressure estimator bandwidth is 
damping ratio is 0.7 sec”^. 


3 rps and the 
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3 . 3 Implementation 


3.3.1 Input-Output 

Estimator state variables (outputs) are linearized variations about nominal or 

a 

commanded states. These state variables are the error angle-of-attack Aff and its 

A a ^ 

time derivative Ad, the error sideslip angle A3 and its time derivative A/1, and the 

A ^ 

error bank angle A 0^ and its time derivative A<P v » (The symbol [ A ] signifies an 
estimated quantity.) 

To obtain variational outputs, the measurements must be differenced with 
their respective desired values to form variational inputs. These variational inputs 
are the measured errors in angle-of-attack A o', body pitch rate Aq, sideslip angle 
A /3, bank angle A0 y , body roll rate Ap, and body yaw rate Ar. (The tilda [ ] 
symbolizes a measured quantity.) The measurements are either read directly from 
sensors or computed (preprocessed) from the sensor data. The body rate measure- 
ments are assumed to be directly accessible from body- mounted rate gyros, but 
the angle measurements are derived from IMU (Inertial Measurement Unit) and 
body-mounted accelerometer information. Commanded values corresponding to each 
of the measurements are either received from the guidance subsystem or computed 
internally. A detailed description of estimator input computation is presented in 
Subsection 2.1. 

3.3.2 Measurement Incorporation 

Current measurements are incorporated into the state estimate immediately 
after they are obtained and preprocessed. This is immediately before the control 
subroutine is called, in order that control commands computed using the most recent 
information may be issued as soon as possible after the DFCS is called. Subsequently, 
the state estimate is extrapolated in preparation for the next DFCS pass, as shown 
in Subsection 3.3.3 below. 

Measurements are incorporated according to the equations 
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where the prime (') indicates a state variable propagated from the previous sample 
instant, is the 2x2 longitudinal gain matrix, and K^ at is the 4x4 lateral gain 

matrix. The gain matrices are updated at the slow repetition rate in the parameter 
estimation routine, which is performed after the state extrapolation step in the 
normal DFCS cycle. 

The flow charts in Fig. 3-1 and 3-2 show the implementation in the DFCS 
program of the measurement incorporation equations and the computation of the 

♦ A 

estimator gain matrices. When dynamic pressure is below q for the longitudinal 
channel and q 2 for the lateral channel, the constant values KC1 and KC2 are loaded 
into the estimator gain matrices as shown; the current values for these parameters 
are 0.45 and 2.5 respectively. For dynamic pressure greater than q^ and q 2 for 
longitudinal and lateral channels respectively, the estimator gain matrices are found 
by linear interpolation on Mach number; the values are given in Tables 3-1 and 3-2. 

3.3.3 State Extrapolation 

After control commands are computed using the current state estimate (from 
Eq (3-1) and (3-2)) and issued, the estimator state is extrapolated to the next sample 
according to the equations: 



a ✓ 

A0 


A <p 


* , 

A <f> 

v 

V J 


k+1 


> - «e < 


lat 


A0 


A <t> 

v 

V J 


(sin <? T ) (-cos 
0 0 


(cos Q' t ) (sin a ) 



(3-4) 


CL CL 

where $ lon and are the longitudinal (2x2) and lateral (4x4) state transition 

matrices 1 ; Ap J , Aqj, Ar J are the average changes in body angular rates about the 
x, y and z axes due to ACPS jet activity over the previous sample interval, and 
is the mean angle-of-attack. 


t C T 

The transition matrices are designated (i.e., <t> closed-loop) to indicate that 

closure of the aerodynamic control loop is taken into account in their derivation: 

.CL x OL , aero ^ aero 
cp = + r l 

where (£> OL is the open-loop state transition matrix, T aero is the aerodynamic- 
control-to-state matrix, and c aero are the linear-feedback aerodynamic control 
gains . 
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Fig. 3-1 Measurement Incorporation for State Estimation 
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Fig. 3-2 Computation of Estimator Gain Matrices 
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Table 3*1 Longitudinal Estimator Gains Scheduled vs Mach Number 



Mach =0.9 

Mach = 2 

Mach = 5 

Mach = 10 

Mach = 20 

(1. 1) 

0.0731 

0.0769 

0.0880 

0.0568 

0.0401 

(1, 2) 

0.0408 

0.0411 

0.0354 

0.0523 

0.0586 

(2, 1) 

0.0338 

0.0353 

0.0329 

0.0207 

0.0116 

(2, 2) 

0.4320 

0.4289 

0.2351 

0.2616 

0. 1909 


Table 3“2 Lateral Estimator Gains Scheduled vs Mach Number 



Mach =0.9 

Mach = 2 

Mach = 5 

Mach = 10 

Mach = 20 

(i, i) 

0.2119 

0.0578 

0.0427 

0.0314 

0.0222 

(1, 2) 

-0.0297 

-0.0794 

-0.1495 

-0.1098 

-0.0942 

(1, 3) 

-0.0115 

-0.0282 

-0.0234 

-0.0210 

-0,0203 

(1, 4) 

-0.0337 

-0.0456 

-0.0637 

-0.0605 

-0.0659 

(2, 1) 

0.2410 

0.0238 

0.0214 

0.0115 

0.0065 

(2, 2) 

0.0257 

0.0155 

0.0016 

0.0177 

0.0075 

(2. 3) 

0.0298 

0.0696 

0.0701 

0.1256 

0.0900 

(2, 4) 

-0.4127 

-0.4213 

-0.2223 

-0.2261 

-0. 1650 

(3, 1) 

-0.0027 

-0.0008 

-0.0015 

-0.0011 

-0.0009 

(3. 2) 

0.4924 

0.4372 

0.3569 

0.3213 

0.2440 

(3, 3) 

0.0148 

0.0187 

0.0198 

0.0224 

0.0256 

(3, 4) 

0.0010 

0.0030 

0.0072 

0.0122 

0.0142 

(4, 1) 

-0.1164 

-0.0212 

-0.0181 

-0.0133 

-0.0070 

(4, 2) 

1.355 

1.130 

0.8698 

0.6984 

0.4032 

(4, 3) 

0.3703 

0.3846 

0.1853 

0.2043 

0.1473 

(4,4) 

0.0387 

0.0725 

0.0783 

0.1229 

0.0903 
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The state transition matrices are computed to be the discrete versions of the 
respective matrices which describe the desired continuous closed-loop aerodynamics 

(3-5) 


$ CL « I + f cl it 


Cl 

where I is the identity matrix, AT is the sample interval, and F is the system 
dynamics matrix; i.e.. 


Ax = F CL Ax 


(3-6) 


where Ax is the variational state. Equation (3 - 5) is the first order power series expan- 

C I 

sion to the solution of Eq (3-6). A more accurate representation of $> J would 

probably be necessary only if AT were not small with respect to the natural time 
CT 

constants of F . 

CT 

The longitudinal state transition matrix is found by Eq (3-5) from Fj^ , which 
takes on either of two forms: 


p CL 

Ion 


(3-7) 


„CL 

*lon 




1 

■ 2 ^d w d 


(3-8) 


where the parameters and specify the desired longitudinal response. (See Sub- 
section 5. 3. 1 . 1 ). Equation (3 -7) is used for dynamic pressure below q^ , when the 
elevator is essentially used only for trim. Equation (3-8) is used for dynamic pres- 
sure above q ^ . 

C I 

The lateral state transition matrix is found by Eq (3-5) from F^ ^ , which takes 


on either of two forms: 


p CL 

lat 


0 0 
0 0 
0 1 
0 0 


(3-9) 
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where the parameters a n through t> 22 specify the desired lateral response 
(see Subsection 5. 3.2. 1.2). Equation (3-9) is used for dynamic pressure below 
q 2 , when the aileron is essentially not used. Equation (3-10) is used for dynamic 
pressure above q ^ . 

The flow charts in Figs. 3-3 and 3-4 show the implementation in the DFCS 
program of the state extrapolation equations and the computation of the state transi- 
tion matrices. The numerical values of the parameters used in the state extrapola- 
tion are indicated in Table 3-3. The control gains C ' , , C t ' C n wL1] 

6a 6a j3 

be obtained from the entry control module in the future. Currently, however, they 
are independently scheduled as functions of Mach number, as shown in Fig. 3-5. 
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Table 3-3: Parameter Values for State Extrapolation 


Aerodynamic and Inertial Parameters 


s 

2690. 

b 

78. 056 

I 

XX 

805000 

I 

zz 

6067000 


Fixed Control Gains 



22 


'll 


12 


21 


22 


0. 2401 
6 . 0 
2. 94 
2. 9008 


See Table 5-1 


Scheduled Control Gains 


See Figs. 5-16 through 5-21 





See Fig. 3-5 
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MACH NO. 
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Fig. 3-5c C* vs Mach Number 






0.28 
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Fig. 3-5d 





3. 4 Preliminary Performance Summary 


Simulated maneuvers have been performed on the Space Shuttle Functional Sim- 
ulator using the state estimator in conjunction with the entry DFCS and the NK00089B 
(February version) vehicle. Both large and small step changes in angle of attack and 
bank were commanded under a wide range of flight conditions (Mach = 2, q = 150 lbs/ft 
to Mach = 25, q = 0. 25 lbs/ft 2 ). 

Comparisons were made between simulations in which the DFCS controller oper- 
ated on state estimator outputs (FILTER = 1) and corresponding simulations in which 
the DFCS controller operated directly on measured data (FILTER = 0). When random 
noise sources were not simulated (perfect measurements), comparisons between 
corresponding simulations were very good in every case. When random noise sources 
were simulated, corresponding simulations exhibited similar mean behavior, but aero- 
dynamic surface commands were smoother for those cases which used the estimator. 
Figure 3-6, for example, shows actual angle of attack, roll angle, and aerodynamic 
surface position time histories for corresponding cases at Mach 14 and a dynamic pres- 
sure of 108 lbs /ft 2 . 
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Fig. 3 -6b Longitudinal and Lateral Responses with Noise Sources 
Simulated and FILTER = 1 (estimator in loop) 
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SECTION 4 


DFCS - CRUISE PHASE 


4. 1 Introduction 

The cruising flight phase of an orbital space shuttle flight begins at about 
Mach 0.9, ends at touchdown, and follows the high-angle-of-attack/high-Mach- number 
entry phase. During the cruise phase the control is accomplished by moving the aero- 
dynamic control surfaces: elevons, rudder, rudder flare (speed brake), and body 
flap. In addition, the cruise phase DFCS has been extended to deal with shuttle pow- 
ered horizontal flight, which encompasses developmental testing and ferry missions. 
During powered flights, thrust is an additional control variable. 

For the cruise phase the DFCS contains three primary modes of operation: 
manual, semiautomatic and automatic. Manual mode consists of two types of control 
systems: Stability Augmentation System (SAS) and Command Augmentation System 
(CAS). By means of a side arm controller and rudder pedals, the pilot uses either of 
these two systems to manually control the orbiter. The pilot can reduce his workload 
while controlling the orbiter by selecting certain automatic functions provided by outer 
loops. This is the semiautomatic mode of operation. In the automatic mode of opera- 
tion guidance algorithms provide the commands to the digital flight control system, 
and the pilot's control role is reduced to monitoring the overall performance. One 
further system available is the Backup Control System (BCS) which is nearly identical 
to the SAS. This system is completely manual, and following its initial flight test, 
is intended to remain unchanged throughout the Space Shuttle Program. 

Speed control is accomplished by using either speed brake or throttle. The 
settings can be controlled in an open loop fashion either by the pilot in the manual 
mode or by the guidance algorithms in the automatic mode. Alternatively the settings 
can be obtained from closed loop speed control loops. In this case, depending on the 
mode of operation, the speed commands come either from the pilot (select/hold mode) 
or from the guidance. For the purpose of pitch trimming, the body flap's setting is 
scheduled as a function of Mach number, which can be supplemented by a manual flap 
command. 

Stabilization of the vehicle in the cruise phase is accomplished with compensa- 
tion filters which appear in forward loop, feedback, and crossfeed paths. It was 
found that for control of the rigid body many of these filters are not needed and good 
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control performance can be obtained with simpler compensations than allowed for. 
However, these filters will become necessary and may, in fact, become insufficient 
in size as control of the flexible vehicle is considered. A table of gains suitable for 
control of the rigid body 89B vehicle at a particular flight condition (Mach No. 0.6, 
dynamic pressure 300 psf, altitude 6000 ft) is provided (Table 4-1). The gains which 
must be varied as a function of some flight parameter (such as dynamic pressure), 
are marked by an asterisk. 

The description of the cruise phase autopilot is presented in the order of 
increasing detail, starting with a discussion of simplified block diagrams, proceed- 
ing to detailed functional block diagrams, and then finally to detailed flow charts. 
Subsection 4.2 describes manual, semiautomatic, and automatic control modes as 
well as speed and body flap controls. The description of the BCS is followed by that 
of mode selection and mode change trim initialization. Filter implementation is 
discussed in Subsection 4.3 where a generalized n^-order digital filter implementa- 
tion is detailed along with low -pass and washout filters, and integration functions. 
Subsection 4.4 describes preflight initialization of gains and filter coefficients, as 
well as limits, deadbands and other items. Subsection 4.5 presents summaries of 
input/ output quantities for various channels such as elevator, aileron, etc. The 
detailed functional block diagrams in Subsection 4.6 are preceded by a list of conven- 
tions in order to facilitate their presentation. Subsection 4. 7 contains functional 
flow charts for initialization subroutine of the cruise DFCS, and then detailed flow 
charts are presented in the Appendix, Subsection 4A. 

4. 2 Control Modes 

4.2.1 Manual Modes of Operation 

For manual operation the two types of control systems are a command augmenta- 
tion system (CAS) for normal operation, and a stability augmentation system (SAS) for 
lower performance downmode operation. The pitch rate and roll rate feedback loops, 
with simple gain compensations in the feedback paths, are employed in SAS (Fig. 4-1). 
The manual commands to the elevator and aileron are summed with these feedback 
signals and sent out to the control surfaces. The manual rudder command is summed 
with two feedbacks and one crossfeed signal from the aileron before it is sent to the 
rudder. Washed-out yaw rate feedback augments the dutch-roll stability. The turn 
coordination in SAS is achieved by either a combination of aileron to rudder crossfeed 
and lateral acceleration, or sideslip feedback. For selecting the best control response, 
pitch rate, roll rate, and yaw rate gains can be scaled by the pilot from his control 
panel. 

A simplified block diagram of CAS is shown in Fig. 4-2. The longitudinal and 
lateral control systems are pitch rate (or C ) loop and roll rate loop closures having sim- 
ilar forward path compensation filters followed by a proportional path and an integral 
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* Gain-switch-sensitive gains 


Fig. 4-1 Simplified Functional Block Diagram of The Stability Augmentation System 
(also Backup Control System) Mode 


4-3 





Fig. 4-2 Simplified Functional Block Diagram of the Command Augmentation System Mode 







bypass. The integral bypass provides automatic trimming and assures unit-gain rate 
tracking. This function can be switched in or out from the control panel, and the pitch 
rate /C loop has an inner damping loop and an outer feedback loop enabling a selection of 
either C or pitch rate. C is a blend of pitch rate and normal acceleration. The 
pitch rate feedback can be fed back directly, passed first through a washout filter, or 
fed back with an added steady turn compensating component. This component is pro- 
portional to the yaw rate and tangent of the roll angle which is measured or derived 
from the measured roll rate. Depending on the type ot feedback (C or pitch rate), 
the manual command is scaled according to either C or q for the control variable. 

The roll rate loop does not have an inner damping loop and thus uses roll rate feed- 
back. The rudder channel is the same as in the SAS. 


4.2.2 Semiautomatic Modes of Operation 

Certain selectable outer loop closures provide semiautomatic modes of opera- 
tion, the construction of which is shown in Fig. 4-3. For longitudinal control the in- 
puts are either pitch attitude or altitude references with the controlled variables being 
pitch attitude and altitude. For lateral control the references and the controlled 
quantities are the bank angle and heading. In both cases, engaging the control stick 
opens these outer loops and enables the control stick steering lunction. Therefore, 
returning the stick to null reengages the selected semiautomatic mode once the body 
rates fall below a threshold value. 


4.2.3 Automatic Modes of Operation 

For one particular cruise phase guidance scheme the automatic mode is sub- 
divided into two regions of operation according to the iollowing table: 


Name of 
guidance phase 

Approach 

Terminal 


Approximate 
altitude range 

40000 to 7000 feet 
7000 feet to touchdown 


Guidance 

commands 



Du ring the approach phase the longitudinal control is of the angle-of-attack. In the 
terminal phase, longitudinal control is accomplished by the pitch rate loop which is 
identical to the CAS loop described in Subsection 4. 2. 1 (including pitch rate compensa- 
tion during steady turns). In both cases the lateral control is of the roll attitude. 
However, the roll angle commands are different in the two guidance cases: roll 
angle is taken about the velocity vector in the approach phase, and about the body 
(vehicle x-axis) in the terminal phase. For both cases the form of the lateral control 
system and the gains are the same. 

As illustrated in Fig. 4-4, guidance commands other than those shown above 
can also be accommodated. Thus the longitudinal control system can be configured 
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Control Stick Steering 
Pilot Pitch Rate 
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Fig. 4-3 Simplified Functional Block Diagram of Outer Loops for Semiautomatic Modes 
of Operation 

















as either an altitude loop or a pitch attitude loop, both of which are outer loops closed 
around the inner CAS loops. Similar to the CAS, an additional longitudinal control 
mode has C feedback. As well as the roll angle loop already discussed, the lateral 
control system is provided with a heading loop closed around the roll attitude loop. 

The rudder loop is implemented similar to the SAS or CAS mode, except that 
a provision is made to accept sideslip guidance commands and direct or washed-out 
yaw rate feedback for automatic decrab maneuver. Rudder pedal commands enable 
this maneuver to be made manually. 

It should be mentioned that, although not shown on the simplified block diagram, 
additional guidance commands can be applied in the automatic mode. These are 
elevator, aileron, and rudder commands applied directly to the appropriate summing 
junctions, and roll rate guidance command to the roll rate summing junction. 

4.2.4 Speed Control 

Control of indicated airspeed for the cruise phase is accomplished by two control 
effectors, the speed brake (also commonly known as "rudder flare”), and the air 
breathing engine system (ABES). Use of the ABES is anticipated only for purposes 
of powered horizontal flight operations, but the appropriate control algorithms re- 
sulting in throttle command are instituted in the current comprehensive cruise DFCS. 
Simplified functional block diagrams of the speed control modes are presented in 
Fig. 4-5. 

In BCS and manual modes, which are functionally identical, the pilot’s manual 
commands are passed directly to the control actuators. Although not indicated in 
Fig. 4-5, speed brake and throttle commands in BCS and manual modes, as in all 
speed modes, are crossfed to the elevator for the purpose of trimming out resulting 
pitching moment. Note that BCS and manual control may employ speed brake or 
throttle or both. 

The IAS-hold mode is functionally similar to the automatic guidance mode which 
provides IAS commands, and is a velocity loop closure with forward loop compensa- 
tion. Currently, the speed brake loop is implemented with proportional plus-integral 
compensation, and the throttle loop has lead-integral lag compensation. In this mode 
either speed brake or throttle is utilized and the nominal elevator mode selection is CAS. 

The automatic (direct guidance) mode issues commands to the speed brake 
and throttle based directly on guidance issued speed brake and throttle commands. 

This may occur for either speed brake or throttle or both. 

Finally, the automatic Terminal Area Energy Management (TAEM) guidance 
mode receives the TAEM guidance rate command and first passes through switching 
logic to decide whether the command is to be issued to speed brake or throttle. If 
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BCS AND MANUAL MODES 
(Speed Brake or Throttle or both) 


Manual 

Speed Brake /- 
Command 


Speed Brake 
Command 


Manual 

Throttle 

Command 


Throttle 

Command 


LAS- HOLD OH AUTOMATIC MODES 
(Speed Brake or Throttle) 


IAS or 
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Velocity 
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Fig. 4-5 Simplified Functional Block Diagrams - Speed Control 
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an indicated airspeed decrease is called for and the throttle is presently opened 
above its minimum, the command is issued to decrease throttle, or otherwise, the 
speed brake is opened. On the other hand, if an indicated airspeed increase is 
called for and the speed brake is presently opened above minimum, the command is 
issued to close the speed brake, or otherwise, the throttle is opened. These com- 
mands are scaled appropriately, integrated, and issued to the controller. 

4.2.5 Body Flap 

Control of the body flap, an aft pitching moment generator, is achieved via 
commands issued according to a piecewise linear body flap versus Mach number 
profile, biased with a direct pilot manual command (see Fig. 4-6). In addition, 
appropriate command biases are instituted to provide for full control surface deflec- 
tion range . 

4.2.6 Backup Control System (BCS) 

The BCS mode is an autonomous, digitally implemented, multi-axis manual 
control capability which serves as pilot requested backup in the event of primary 
DFCS failure. Although the design presented here is only of a preliminary nature, 
it is intended to demonstrate the capabilities currently envisioned for BCS mode 
plus techniques for their implementation. 

Whereas the primary DFCS is likely to be subject to modifications throughout 
the flight testing and orbital flight operations, the intent for BCS is to obtain an 
acceptable design prior to the first flight and to leave it unmodified thereafter. As 
the only BCS coding which is integrated into the primary DFCS structure is the BCS 
call flag from the executive routine, possibility of subsequent faulty modifications is 
virtually eliminated, thus providing the pilot with a continuously reliable alternative. 

BCS mode is functionally similar to primary manual mode and shares the same 
simplified functional block diagrams (Figs. 4- 1, 4-5, and 4-6). However, there are several 
fundamental differences which are not indicated: The axis- by-axis, multi- routine 
format of the primary DFCS is replaced by a single routine which processes all 
axes in succession; and the filter update, control and filter pushdown (wrap-up) 
functions are imbedded in the proper sequence. Whereas the primary system will 
have some gains scheduled as a function of dynamic pressure, a single set of fixed 
gains which provides suitable performance over as large a cruise flight spectrum 
as possible will be utilized for BCS. Finally, pilot signal preprocessing is not in- 
cluded for elevator,, aileron, and rudder, and there is no Mach number profile for 
body flap. 
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Fig. 4-6 Functional Block Diagram of Body Flap 
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4.2.7 Mode Selection 


Mode select /indicate interfaces with the pilot are helpful in understanding 
the capabilities of the cruise phase DFCS, especially where manual control and a 
multiplicity of modes complicate the picture. Figure 4- 7a is compatible with simul- 
taneous multi-axis mode selection in contrast to axis-by-axis mode selection. Thus 
it is compatible with earlier versions of the cruise phase DFCS wherein axis-by-axis 
mode changes were not considered. Multi-axis mode selection is not unlikely for the 
ultimate cruise phase DFCS version, considering such things as control loop design 
failure modes, pilot training, and computer input /output simplifications. The gain- 
select switches and integrator moding select /indicators are on a separate panel. 

This panel might be associated with the ABES console since these capabilities are 
intended primarily for the flight test phases of development. 

Figure 4- 7b, on the other hand, suggests an axis-by-axis mode selection 
capability compatible with the present cruise phase DFCS design, including all 
peripheral moding capability. All indicated functions are implemented, except for 
the pilot selectable pitch turn-coordination-override capability. Gain-select switches 
are active only in the SAS mode. Pilot-selectable integrator-moding-override capa- 
bility is provided only in the CAS mode (pitch and/or roll). 


4 . 3 Digital Filters and Their Implementation 

4.3.1 Generalized n^- Order Filters 

Control filters for the cruising mode are classical, recursive, n^-order 
digital filters, whose Z transform is D(Z): 


D(Z) 


a Ca 2 Z~ 1+ ••• + a 2 n Z~ n 
1 + a 3 Z _1 + ... +a 2n+1 Z' n 


The general form of their implementation is presented in block diagram format 
in Fig. 4-8. The symbol x represents the controlled variable and the a's are either 
fixed or scheduled gains (currently fixed) which correspond to the gains of form 
Gx T which appear in the flow charts and block diagrams presented later in this 
section. 
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Fig. 4 -7a Mode Select/ Indicate Panels, Type 1, Cruise Phase DFCS 
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Fig. 4- 7b Mode Select/Indicate Panels, Type 2 (Axis-by- Axis) Cruise Phase DFCS 
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Fig. 4-8 N -Order Digital Filter 











nci 

The current filter implementation provides for no greater than 2 -order filter 
capability, but extension to higher order filters (if needed) is straightforward. The 
Z 1 calculation blocks are achieved in the DFCS programming by the relative sequence 
of instructions, which permits the DFCS cycle time lag to occur at the point so desig- 
nated. In the coding, each 2 nd -order state filter is divided into two parts in order to 
minimize the computational transport lag between data input and control command. 

The first part, which occurs in the Filter Update Routine is limited to updating the 
filter estimate according to the most current measurement: 

x , = a„ x. + x^ T , 
out 1 in N1 

Where represents the n th -node, or summing junction, of the filter. Note that 
exactly n nodes are reouired to implement an n^-order filter. 

The Filter Update Routine is sequentially followed by the Control Routine which 

issues control commands based upon the most current filter output value x . Filter 

out 

data "pushdown", which is actually the tabulation of the current value associated with 
each node, is then accomplished in the Filter Pushdown Routine by the following two 
equations for a 2nd-order filter: 

x tvti = x. - a n x + x x n 
N1 2 in .1 out N2 

x tvtd = a /i x - “ a c x , 

N2 4 in 5 out 

This sequence of instructions results in the following equivalent difference 
equation for the filter output: 

x , “ a i x -^ + a o x ; + a , x. -a 0 x , -a r x 
out 1 in 2 in_i 4 m_ 9 3 out_ 1 5 out_ 2 


Where subscripts -1 and -2 denote past values, one and two samples previous. 

Control loops which utilize 2 nd -order filter capability are available for control 
of the following quantities: 


BCS 

Manual 

Modes 

Automatic and 
Semi-Automatic 
Modes 

Feedback Quantities 

Others 

roll rate 

roll rate 

roll angle 

measured pitch rate 

turn coordination 

pitch rate 

pitch rate 

pitch angle 

measured yaw rate 

rudder crossfeed 

yaw rate 

yaw rate 

yaw angle 

measured lateral 
acceleration; 




angle -of-attack 

sideslip angle 

indicated air- 
speed 

altitude 

normal acceleration 
(with gravity) 
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4.3.2 Specialized l St -Qrder Filters 

Recursive, l st -order washout filters are utilized in the pitch rate and yaw 
rate feedback paths of the primary system and in the yaw rate feedback path of BCS. 
The Z transform of the filter is: 


D(Z) = 


1 - z ' 1 


1 + G „,Z 
xW 


-1 


A block diagram and the update and pushdown equations for the washout filter 
are presented below 



Fig. 4-9 Digital Washout Filter 


Filter Update : 

Filter Pushdown: x. 


x , = x. + x, , 
out in N 1 


N 1 


-x. - G , Tr x , 

m xW out 


. st 


Recursive, l^-order low pass filters are utilized for the filtering of the fol- 
lowing quantities: 


1) pilot stick inputs 

2) crossfeeds to elevator (from speed brake, throttle, and body flap) 

3) elevator trim tracker 

4) yaw rate for pitch turn coordination 


The Z transform of the low pass filter is 


D(Z) = 


1 + G 


1 + G Z 
x 
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and the block diagram and update and pushdown equations are presented below 



Filter Update: x , = x. + G x X7i 

out in x N1 

Filter Pushdown: x XT1 = -x , + x. 

N1 out in 


4.3.3 Integrators 

Integration is achieved by both trapezoidal and rectangular integrators. Trape- 
zoidal integrators are employed for the integral bypass functions in the elevator and 
aileron loops, while rectangular integrators are utilized for functions in which trape- 
zoidal integration proves to be of no additional advantage. Thus, rectangular inte- 
grators are utilized for the integration of pilot trim discretes (elevator, aileron, and 
rudder in primary system and BCS) and for the integration of TAEM guidance rate 
commands (speed brake and throttle). The Z transforms of these integrators are 

Tp T - 1 

Rectangular: D(Z) = j- Trapezoidal: D(Z) = — — — 2 

1-Z' 1 2 2 . z -l 


and the block diagrams and update and pushdown equations are as follows: 

limiter limiter 




Fig. 4-11 Digital Integrator Filters 

Rectangular Integrator Trapezoidal Integrator 



Rectangular 

T rapezoidal 

t f 

Filter Update: 

x out X out + T F x in 

x out 2 X in f X N1 

T 

Filter Pushdown: 

none 

F 

x .T i 3 — r— x. + x 
N1 2 in out 
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Note from the block diagrams that the output of the integrators is limited; in 
practice this occurs in the filter update routine. Also notice that rectangular inte- 
gration requires no pushdown equation, thus substantiating its use whenever possible. 

4 . 4 Preflight Parameter Specifications 

Basic parameters for the cruise phase DFCS include filter coefficients, gains, 
limits, deadbands, logic-control flags, etc. Each of these basic parameters must 
be numerically specified prior to the cruise phase initialization pass. A representa- 
tive preflight parameter specification is presented alphabetically: Table 4- la for 
non-subscripted parameters and Table 4- lb for subscripted parameters. All numeri- 
cal data is based upon a single-point design at Mach 0.6 and 300 psf (about 6000 feet) 
utilizing aerodynamic data for the February 1973 AEROIB version of the 89B orbiter . 
The RC AH version of the CAS elevator mode is indicated (ICSTAR = 0); the rudder 
sideslip feedback option is chosen (IBFDBK = 1); pitch turn coordination is bypassed 
(IQTRN = 1) . 

Most of the basic parameters are utilized directly, but a number of ’’working” 
quantities are employed, such as sample-interval-multiplied trim rates or ^-guidance 
values of CAS filter coefficients, limits and loop gains. These are computed on the 
initialization first pass (Fig. 4-24) or upon mode initialization as required (Fig. 4-25 
through 4-28). 

Earlier orbiter cruise phase DFCS designs have utilized scheduled gains and 
coefficients to provide desirable response characteristics over the cruise phase 
flight envelope (including Powered Horizontal Flight). The scheduling parameter 
has been dynamic pressure. The scheduling operations (when enabled) are carried 
out every slow pass by the Parameter Estimation Routine (Fig. 4-12), which operates 
on preflight- specified tables of parameter value vs dynamic pressure. Typically 
seven data points are used: 100, 150, 225, 300, 425, 525, and 650 psf. The thirteen 

parameters corresponding to the scheduled quantities of earlier designs are identi- 
fied by asterisks in Table 4-1 . Tabular specifications for these gains and coefficients 
will be added to the preflight parameter specification tables, as they become avail- 
able . 

4.5 Input /Output Summaries 

Figures 4- 13a through 4-13e detail those quantities which are input to and 
output by the five basic control channels (speed, elevator and body flap, aileron, 
rudder and BCS). Input includes information transferred across the various inter- 
face boundaries — executive, guidance, and pilot, as well as miscellaneous moding 
information, crossfeeds from other loops, and sensor measurements. Output con- 
tains control effector commands, pertinent mode status for display, and crossfeeds 
to other loops. 
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Table 4- la Preflight Parameter Specifications for the Cruise Phas 

(Non-subscr ipted Parameters) 
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Scheduled Gains 



Table 4-la Preflight Parameter specifications for the Cruise Phase (Cont) 
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Scheduled Gains 



Table 4- lb Preflight Parameter Specifications for the Cruise Phase (Cont) 


(Subscripted Parameters) 


OFCMCH(I,J) (-.2094, .1), (-.2094, .2), (-.2094, 

I-.OH73, .52), (-.0873, 6), (.1745, 
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.01, o . , n . 
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Fig. 4-12 Parameter Estimation Routine for Cruise Phase 
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Fig. 4- 13b Elevator and Flap Input/Output Interface Variables 
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^ ig. 4- 13c Aileron Input/Output Interface Variables 
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Fig. 4- 13d Rudder Input/Output Interface Variables 




GUIDANCE 
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Fig. 4-13e Backup Control System Input/Output Interface Variables 



4 . 6 Functional Block Diagrams for Cruise Phase DFCS 


This section presents a detailed control-engineer’s description of the cruise 
phase DFCS, using functional block diagrams. These include all the feedbacks, 
crossfeeds, limiters, switches, and compensation filters currently available with 
preflight parameter specification options. 

A number of simplifying conventions have been employed, particularly in the 
areas of filter and limiter descriptions. These conventions are summarized in 
Fig. 4-14. 

Tables 4-2 through 4-4 indicate the realizable cruise phase DFCS modes of 
operation, and the corresponding settings of all flags. Note for example, that while 
MODS = 3 (speed brake usin^u- filter for IAS-hold) and MODT = 3 (throttle using 
u-filter for IAS-hold) are both valid modes of operation, they must not occur simul- 
taneously. The mode change Initialization Subroutine (Fig. 4-25 through 4-28) 
handles all the setting of all derived flags, according to the logic indicated in these 
tables . 

Pilot input signal processing for speed brake throttle and flap controllers, for 
pitch and roll sidearm stick, and for pedals, is shown in Fig. 4-15. Pitch, roll and 
yaw manual trim circuits are shown in Fig. 4-16. Outer loop (semi-automatic or 
automatic) error signal generation is shown in Fig. 4-17. Note that the pilot '’Select" 
option is implied (except for altitude) simply by "writing” into the appropriate storage 
location, e.g., H^pp* 

Speed control, including throttle and speed brake options, is shown in Fig. 

4- 18a, which indicates an input from TAEM-guidance interface signal processing as 
shown in Fig. 4- 18b. Elevator control loops are shown in Fig. 4- If). Also shown 
in Fig. 4-10 are the generation of for turn coordination (pitch and yaw rate 

command circuits), the generation of body flap commands (manual override of a 
Mach profile), and an elevator trim-tracker network akin to conventional stabilizer 
trim-tracking. (This trim indication is being studied for trim- initialization on 
emergency downmodes.) 

Figure 4-20 presents the aileron control loops, and Fig. 4-21 presents the 
rudder control loops, including provision for both 0 feedback and "0" feedback (a 
combination of yaw rate command, yaw rate feedback without washout filtering, 
lateral acceleration, and aileron) . Note that pedal and trim are available even in 
auto mode (through the gain GRXAUT) for manual decrab maneuvers. 

Figure 4-22 indicates the override capability of the BCS mode, indicated by 
MODB = 1. Details of BCS loops (essentially the corresponding MANUAL or SAS 
loops with independently specified gains, coefficients, and nonlinearities) are given 
in Fig. 4-23. 
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Fig. 4-14 Simplifying Conventions for Block Diagrams 
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Fig. 4-14 Simplifying Conventions for Block Diagrams (Cont) 
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Table 4-2 Logic for Speed Control Flag Settings 


Input Flags 

Derived Flags 


MODU 

MODUG 




MODUFP 

BCS, UBCS 

1 

X 

1 

0 

X 

X 


1 

X 

1 

1 

X 

X 

UMAN 

2 

X 

2 

0 

0 

1 

(IAS - MAN) 

2 

X 

2 

2 

0 

1 

UASB 

3 

X 

3 

0 

1 

2 

(IAS - SB) 

3 

X 

3 

2 

1 

2 


L 


1 

5 

0 

0 

1 




1 

5 

5 

0 

1 




2 

4 

0 

1 

3 




2 

4 

2 

1 

3 




3 

4 

5 

1 

3 




4. 

2 

4 

2 

3 

AUTO, UAUT 



5 

5 

4 

2 

3 




6 

6 

0 

0 

4 




6 

6 

2 

0 

4 




7 

6 

5 

0 

4 




8 

2 

6 

0 

4 




9 

5 

6 

0 

4 




10 

6 

6 

0 

4 

UATH 

5 


X 

2 

3 

2 

2 
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Table 4-3: Logic for Elevator Control Flag Settings 



BCS, EBCS 


SAS, ESAS 


CAS, ECAS 


Input Flags 


MODEG 


EATT 


EALT 



Derived Flags 


ITHC 




Table 4-4 ; Logic for Aileron Control Flag Settings 


Input Flags 


Derived Flags 


MODA MODAG IPC IPHC 


ABCS 

1 

X 

A SAS 

2 

X 

ACAS 

3 

X 


4 

1 

AUTO 

4 

2 


4 

3 

AATT 

5 

X 

AHDG 

6 

X 
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Fig. 4-15 Block Diagram, Pilot Input Signal Processing 
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Fig. 4-16 Block Diagram, Pilot Manual Trim Circuits 
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Fig. 4-17 Block Diagram, Outer Loop Error Signal Generation 
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Fig. 4-18b TAEM Guidance Processing for Speed Control 
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Fig. 4-20 Aileron Control Loops 
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Fig. 4-21 Rudder Control Loops 





















NOTE: 

BCS mode selection 
is "across-the-board" 

(MODB = 1: MODX = 1 
where X = U, E, A, and R) 
rather than axis-by-axis. 

BCS is activated by 
selecting the BCS mode 
in any one axis (MODXR = 1, 
where X = U or E or A or R) 
or by selecting 4-axis BCS 
(MODUBR - 1). 


combined to generate 
elevon commands 
5 EL c and 5 ERq 
in the control output 
subroutine 
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Fig. 4-22 Control Output Signals 
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Fig. 4-23 Backup Control System 
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Fig. 4-23 Backup Control System (Cont) 














4 . 7 Functional Flow Charts - Cruise Phase Initi alization 

Figures 4-24 through 4-28 represent preliminary attempts to specify the details 
of the cruise phase DFCS at the engineering/functional level rather than at the coded 
software level presented in Appendix 4A . Implementational freedom (especially with 
respect to programming language syntax) should not be restricted by software over- 
specification. Hence, the flow charts of this section will be augmented, until such 
time as a total engineering/functional specification is achieved. Figures 4-25 through 
4-28 interface with the skeletal structure shown in Fig. 4a- 2. 
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Prepare working registers for integrator 
input gains (premultiplied by TF for 
rectanguLar or TF2 for trapezoidal), and 
for rate limits (premultiplied by TF): 


SPEED 


TRIMS 


T F 4s ma , X T p *6t max , T p * xt max (where xt max = 6t max /GTX), 
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(load 1 1 * for automatic integrator 
IEINOK = IAINOK = 1 activation at gear-retract; 

load *0* for manual activation) 


END) 


Fig. 4-24 Functional Flow Diagram, Load Manipulation Routine 
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Fig. 4-25 Functional Flow Diagram, Speed initialization Routine 
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Fig. 4-26 Functional Flow Diagram, Elevator Initialization Routine 
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Aileron Initialization Routine 
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Fig. 4-28 Functional Flow Diagram, Rudder Initialization Routine 





APPENDIX 4A 


DETAILED FLOWCHARTS FOR CRUISE 
PHASE DFCS 


Detailed specifications for the cruise phase DFCS, in terms of FORTRAN 
flow charts of an operating version of the DFCS, are presented in Figs. 4A-1 through 
4A 20. Minute functional details are covered, but unfortunately they are often over - 
whelmed by implementational aspects (e. g. , higher-order-language syntax), and are 
thus less-than-ideal for specification purposes. 

Figure 4A-1 depicts the cruise phase preflight (or first-pass) load manipulation 
operations. F'.gure 4A-2 presents the skeletal structure into which the remaining 
figures fit. In particular, there are five basic subroutines. These subroutines 
interface directly with the Dt CS executive as shown in Section 1. 2. Their functions 
for cruise phase are as follows: 

Initialization Subroutine (Fig. 4A-3 through Fig. 4A-6) 

handles all first-pass (preflight turn-on or inflight phase-change) and 
mode change initialization activity. 

Filter Part 1 Subroutine (Fig. 4A-7 through Fig. 4A-10) 

handles all preparatory calculations for the Control Output Subroutine, 
including pilot input signal processing, filter output node computations, 
control stick steering override detection, etc. 

Control Law Subroutine (Fig. 4A-11 through Fig. 4A-13) 

handles all control output calculations in preparation for the Control Out- 
put Subroutine (Subsection 2. 2). 

Filter Part 2 Subroutine (Fig. 4A-14 through Fig. 4A-18) 

handles all filter wrapup processing and miscellaneous activity, such as 
body flap command calculations, bank-angle filter calculations, and re- 
engaging semiautomatic operation following a control stick steering override. 

BCS Subroutine (Fig. 4A-19 through Fig. 4A-20) 

handles sequentially, the Filter Part 1, Control Law, Control Output, and 
Filter Part 2 calculations for a 4-axis BCS mode of operation. 


4-52 




(*) For automatic integrator -activation at gear-retract; 
for manual reactivation load IEINOK = LAINOK = 0. 


Fig. 4A-1 Load Manipulator Routine Cruise Section 
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Fig- 4A-3 Initialization Subroutine, Cruise Phase Section, Speed Subsection 
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Fig- 4 A- 4a Initialization Subroutine, Cruise Phase Section, Elevator Subsection 
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Fig. 4A-4b Initialization Subroutine, Cruise Phase Section, Elevator Trim Subsection 















Fig. 4A 5a Initialization Subroutine, Cruise Phase Section, Aileron Subsection 
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Fig. 4A-6b 
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Fig. 4A-6a Initialization Subroutine. Cruise Phase Section. Rudder Subsection 
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Fig. 4A-6b Initialization Subroutine, Cruise Phase Section, Rudder Trim Subsection 







Fig. 4A-7a Filter Part I Subroutine, Cruise Phase Section Speed Subsection 


4-62 





















Fig. 4A- 7b Filter Part I Subroutine, Cruise Phase Section Speed Subsection, 
(TAEM Guidance Routine) 
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Fig. 4A-8a Filter Part I Subroutine, Cruise Phase Section, Elevator Subsection, 
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Fig. 4A-8b Filter Part I Subroutine, Cruise Phase, 
Elevator Subsection, Outer Loops 
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Fig. 4A-8c Filter Part I Subroutine, Cruise Phase, /A 

Elevator Subsection, Feedback (END 
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Fig. 4A-9a Filter Part I Subroutine, Cruise Phase Section, Aileron Subsection 
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Fig. 4A-9b Filter Part I Subroutine, Cruise Phase Section, 
Aileron Subsection, Outer Loops 
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Fig. 4A-10 Filter Part I Subroutine, Cruise Phase Section, Rudder Subsection 
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Fig. 4A-11 


Control Law Subroutine, Cruise Phase Section, Speed Subsection 
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Fig. 4A-12 Control Law Subroutine, Cruise Phase Section, Elevator Subsection 
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Fig.4A-13 Control Law Subroutine, Cruise Phase Section, 
Aileron and Rudder Subsections 
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Fig. 4A-14 Filter Part II Subroutine, Cruise Phase Section, Speed Subsection 





Fig. 4A-15 Filter Part II Subroutine, Cruise Phase Section, Elevator Subsection 
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Fig. 4A-17 Filter Part II Subroutine, Cruise Phase Section, Rudder Subsection 
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Fig. 4A-18 Filter Part II Subroutine, Cruise Phase Section, 
Gear /Integrating Moding 
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Fig. 4A-19a Cruise Phase BCS Subroutine, Speed, Elevator, and Flap Subsection 
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Fig.4A-19b Cruise Phase BCS Subroutine, Aileron and Rudder Subsection 
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XTPPB = XTPB 


QMNB (1) = QMNB (2) + GQM 
QMNB (2) = + GQM] 

B (2) * QM - GQMB (3) QMFB 
B (4) ' QM - GQMB (5) * QMFB 



PM MB (1) = PMNB (2) + GPMB (2) * PM - GPMB (3) * PMFB 
PMNB (2) = + GPMB (4) * PM - GPMB (5) * PMFB 

r — 



RENB (1) = RENB (2) + GREB (2) * RM - GREB (3) -:= REFB 
RENB (2) = GREB (4) * RM - GREB (5) * REFB 

RW1B = - REFB -GRVVB * RWFB 



Pig. 4A-20 Cruise Phase BCS Subroutine, Wrap up Subsection 
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SECTION 5 


ENTRY PHASE 


5 . 1 Introduction 

The DFCS design presented in this section is intended to provide stability-axis 
attitude control for the vehicle during entry, where this mission phase begins outside 
the sensible atmosphere (e.g., at approximately 400,000 ft altitude) and ends in the 
high subsonic flight regime. The control logic makes use of both the ACPS and the 
aerodynamic control surface torques, blended together in such a manner as to mini- 
mize ACPS jet activity, while providing capability for simultaneous parallel use 
when necessary. 

The basic design features include gain scheduling and control submode switch- 
ing (both as functions of selected trajectory parameters) to account for variations in 
the vehicle dynamics and to allow for reasonably constant closed-loop transient 
response characteristics. In addition, except for very low dynamic pressure opera- 
tion, maximum use is made of the aerodynamic control surfaces to minimize ACPS 
fuel. Finally, the controller synthesis, and thus the design presented here, makes 
advantageous use of the vehicle's aerodynamic torques resulting from trim attitude 
deviations, allowing the vehicle's rigid body aerodynamics to aid in attitude maneu- 
vers . 

This section is organized into three subsections. Subsection 5.2 provides an 
overview of the entry DFCS organization and operational characteristics, while Sub- 
section 5.3 is a detailed description of the entry DFCS, both functionally and as 
implemented in the present FORTRAN simulator. Finally, Subsection 5 . 4 is a brief 
summary of simulation experience to date. 

5 . 2 Overview 

Figure 5-1 provides an overview of the entry phase control logic structure. 

The interface, state estimation, and air data routines are described in Section 2 
(Subsection 2. 1), Section 3, and Section 2 (Subsection 2.2), respectively. The pur- 
pose of the first two is to provide estimates of the stability axis angles and rates, 
along with the guidance commands, while the purpose of the latter is to provide 
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Fig. 5-1 Overview of Entry Control 





estimates of the trajectory parameters (trim angle-of-attack, Mach number, and 
dynamic pressure). 

longitudinal control of angle-of-attack and pitch rate is provided by elevator 
(tandem elevon) deflections and pitch ACPS jet firings. The elevator logic uses con- 
ventional rate and position feedback (with feedforward trim integration) through 
scheduled gains to maintain relatively constant closed-loop pitch dynamics, while the 
ACPS logic uses phase-plane control to maintain a desired limit-cycle operation. 

lateral control of sideslip and bank angle is provided by aileron (differential 
elevon) and rudder deflections, and roll and yaw ACPS jet firings. The aileron/ 
rudder logic uses rate and position feedback of the lateral variables (with feed- 
forward trim integration for center- of- gravity offset compensation) to provide atti- 
tude and rate control of both sideslip and bank. The ACPS logic uses phase-plane 
control to maintain a desired limit-cycle operation when there is low dynamic 
pressure, and sideslip rate damping to maintain lateral stability when there is high 
dynamic pressure. 

In both the lateral and longitudinal channels, blending of control by the surfaces 
and the ACPS is provided through a logic which determines ACPS inhibition on the 
basis of the magnitude of the current surface command, and on the basis of the 
vehicle's Mach number and angle-of-attack. 

5 . 3 Description of Control Laws 

This section presents the control system block diagrams, equations, system 
parameters, and detailed flow charts defining the entry DFCS design. For conven- 
ient reference, this section is divided into four subsections. Subsection 5.3.1 and 
5.3.2 describe the longitudinal and lateral control laws, while Subsection 5.3.3 
specifies the fixed and scheduled parameters used in the present design. 

5.3.1 Longitudinal Control 

Longitudinal control uses both elevator (tandem elevon) and pitch ACPS jets 
to maintain the commanded angle-of-attack. As shown in Fig. 5-2, which illustrates 
a typical sequence of control activity during entry, there are two control submodes for 
the elevator and one for the pitch jets. 

Early in the entry, when the dynamic pressure is low, the elevator is used 
strictly for trim, so as to minimize pitch disturbance torques and thus fuel usage. 

The trim control law operates in a closed-loop fashion by slowly driving the elevator 
in the direction which nulls the pitch disturbance. The estimate of the disturbance 
torque is obtained by observing the pitch ACPS firing history. During this flight 


5-3 









regime, the pitch ACPS jets, controlled by an angle-of-attack /pitch- rate phase- 
plane logic, provide attitude and rate control in response to guidance commands. 

There is no jet inhibition in this regime. 

When the dynamic pressure gets high enough, (q^), the elevator is used for both 
trim and transient control. Trim is initially maintained through the use of the trim 
logic of the previous submode. At a slightly higher dynamic pressure, (q^ ), trim is 
maintained through the use of an integrator acting upon the angle-of-attack error; 
the initialization value of the integrator is obtained from the value generated by the 
previous trim logic. Elevator control of transient errors is maintained through 
feedback of angle-of-attack and pitch rate, with appropriate gains. During this 
flight regime, pitch jet control identical to that of the previous regime may be uti- 
lized together with the elevator, or control with the jets may be inhibited. This is 
done by the blending logic which is essentially a two-sided deadbanded relay with 
hysteresis. The commanded elevator deflection is compared with preset fractions 
of its maximum and minimum values; should the command be ’’small 1 ’, it is assumed 
that the elevator has sufficient control authority, and thus the pitch jets are not 
required. Conversely, should the elevator command be ’’large”, the pitch jet control 
is not inhibited, but is allowed to assist the elevator in controlling the vehicle. This 
blended control continues throughout the rest of the entry, so that the pitch jets are 
always available for control assistance. 

The logic provides for a gradual turn-on of the elevator control so as to avoid 
a switching transient when the second control submode is entered. 

The remainder of this section describes the longitudinal control laws, and their 
digital implementation, in detail. The block diagrams and equations defining the 
control logic for the elevator and pitch jets are given in Subsections 5. 3. 1. 1 and 
5. 3. 1.2. Subsection 5.3. 1.3 consists of the controller coding summarized in flow 
chart form. 

5.3. 1.1 Elevator Control 

Figure 5- 3a is a block diagram of the elevator control loop, showing the two 
control submodes, with switching determined by q^. 

When the dynamic pressure is less than q^, only an elevator trim setting is 
sent to the elevator actuator. This trim setting is determined by integrating the pitch 
ACPS commands (u^), or » equivalently, by slowly driving the elevator so as to null 
the longitudinal disturbance torques. 

When the dynamic pressure reaches q^, closed-loop control of the angle-of- 
t See Subsection 5.3.3 for numerical values of the parameters 
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attack is provided by conventional position (o') and rate (o’) feedback through a net 
work which attempts to maintain a constant transient response throughout the flight 
envelope. Specifically, the desired closed-loop transfer function is chosen to be a 
well-damped second-order system, or 



desired 


S 2 + 2i; d <o d s + « d 2 


(5-1) 


where (£ d , u> ) are specified for acceptable performance (k^ is not explicitly speci- 
fied). The gains shown in the diagram may be separated into one of three types: 

(1) prespecified by the desired closed-loop characteristics; (2) dynamic pressure 
dependent (l/q); or (3) scheduled gains which are functions of the aerodynamic flight 
regime (e.g., 4^)! Trim control is accomplished in this submode either by the 
trim logic of the previous submode (if q < q’j) or by a conventional clamped trim 
integrator acting on the angle-of-attack error, thus providing a steady- state trim 
elevator command in the absence of transient angle-of-attack errors. 

As can be seen from Fig. 5-3a and 5-3b, in order to allow a gradual blending 
in of active elevator control, the elevator command limits are functions of dynamic 
pressure which gradually open up from the trim setting to the physical deflection 
limits . 


5. 3. 1.2 Pitch ACPS Control 

Shown in Fig. 5-4 is the ACPS phase-plane switch logic representation used 
for control of the pitch jets. For longitudinal control, the phase-plane coordinates 
are angle-of-attack error and angle-of-attack rate, or 

(e i , e 2 ) = d) (5-2a) 

The phase-plane is separated into three regions (by the switch curves), in which the 
commanded jet firing is either positive (u = +1), negative (u = -1), or zero (u 0). 
Thus, for the particular application of pitch jet control 

(u, U) = Uy U y ) (5- 2b) 

where u and U symbolize the jet command and the available acceleration level, 
respectively . 

The equations for the switch curves themselves are based on the fuel-time 
optimal solution to the double integrator control problem and are modified by 


t See Subsection 5.3.3 for numerical values of the parameters 
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deadband incorporation to ensure practical limit-cycle convergence (see typical 
trajectory of Fig. 5-4). The switch curves are defined by: 


y±: e 1 =F (6 i ' 2U e 2 * 
r± : e 1 = ± + 2U C 2 ^ 


(5- 3a) 


where, in this case, the deadbands are pitch attitude deadbands, or: 

6 2 > = 6 « 2 > (5 ' 3b) 
and cr is a fuel-time weighting constant greater than unity given for pitch control by: 

a = a (5-3c) 


The interface (or blending) logic between commanded pitch jet firings and the 
commanded elevator deflections is shown in Fig. 5*5. The approach taken is to 
inhibit firing when there is "sufficient" pitch control acceleration from the elevator. 
This measure of sufficiency is obtained by comparing commanded elevator (Se^, 
shown in Fig. 5-3a) with threshold values to determine whether or not "excessive" 
elevator deflection is being called for. The hysteresis path included in the blending 
interface avoids chatter in elevator and ACPS activity due to cross-coupling effects. 
Note that this design includes an additional switch to preclude jet inhibition when the 
dynamic pressure is low and the elevator is used solely for trim. 


5 . 3 . 1 . 3 Digital Implementation 

Coding implementation for longitudinal control during the entry phase is 
summarized in flow-chart form in Fig. 5-6. 


5.3.2 Lateral Control 

Lateral control uses both aileron (differential elevons) and rudder, and roll 
and yaw ACPS jets to maintain small sideslip angles and follow commanded bank 
maneuvers. As shown in Fig. 5-7, which illustrates a typical sequence of control 
activity during entry, there are three control submodes for the aerodynamic control 
surfaces, one for the yaw jets, and three for the roll jets. 

Early in the entry, when the dynamic pressure is low, the aileron and rudder 
commands are zero. During this regime, the yaw jets, controlled by a bank angle/ 
yaw- rate phase-plane logic, provide attitude and rate control of the bank angle in 
response to guidance commands. The roll jets are controlled by one of two logics, 
depending on the yaw jet activity. When there are no yaw commands, the roll jets, 
controlled by a sideslip /sideslip- rate phase-plane logic, maintain small sideslip 
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Fig. 5-6a Logic Flow for the Longitudinal Control Channel 
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Calculate Elevator Command 



Fig. 5-6b Logic Flow for Elevator Control (Cont) 
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Fig. 5-6c Logic Flow for Pitch ACPS Control 
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Fig. 5-7 Controller Activity during Entry (Lateral) 






attitudes and rates. When the yaw jets are commanded to maneuver the vehicle, 
the roll jets are controlled by a sideslip/sideslip-rate phase-plane logic, which is 
similar to the first submode except that compensation for yaw jet torques is in- 
cluded. Again, the purpose of this logic is to maintain small sideslip attitudes and 
rates. There is no jet inhibition in this flight regime. 

Later in the entry, when the dynamic pressure reaches a certain prespeci- 
fied value (q 3 ), the aileron and rudder commands are still zero. The yaw jets are 
controlled by the same logic as in the previous regime. The roll jets are controlled 
by one of two logics, again depending on yaw jet activity. When there are no yaw jet 
commands, the roll jets are commanded by a deadbanded relay logic which simply 
provides rate damping in sideslip. When the yaw jets are commanded to maneuver 
the vehicle, the roll jets are commanded by the same corresponding roll jet logic of 
the previous regime. There is no jet inhibition in this flight regime. 

Still later in the entry, when the dynamic pressure is higher (>q 2 >, the aileron 
is used for both trim and transient control. Trim is initially maintained at zero. 

At a slightly higher dynamic pressure (q 2 ), trim is maintained through the use of an 
integrator acting upon the bank angle attitude error; the initialization value of the 
integrator is zero. Aileron control of transient errors (in a reverse sense, as it 
has come to be labeled ) is maintained through feedback of bank and sideslip attitudes 
and rates with appropriate gains. During this regime, the rudder command is null 
until the Mach number reaches a certain level (77?^, after which it is proportional 
to the commanded aileron deflection from trim. The proportionality variable is 
a constant. Also during this regime, roll and yaw jet control identical to that 
of the previous regime may be utilized together with the aileron and rudder, or 
control with the jets may be inhibited. This is done with the same type of blending 
logic used longitudinally; here, the commanded aileron deflection determines jet 
inhibition. Also as in the longitudinal case, there is a gradual turn-on of the allowed 
aerodynamic control authority. 

Still later in the entry, when the vehicle has crossed a prespecified curve in 
the Mach-alpha plane, only the aileron and rudder are used for transient control, 
the jets being inhibited for the remainder of the entry. Aileron and rudder control 
(in a "proverse" sense to distinguish it from that of the previous submode) of 
transient errors is maintained through feedback of bank and sideslip attitudes and 
rates (to both surfaces), with appropriate gains. Lateral trim for this regime has 
not yet been implemented. 

The remainder of this section describes the lateral control laws, and their 
digital implementation, in detail. The block diagrams and equations defining the 
control logic for the aileron and rudder are given in Subsection 5.3.2. 1, while those 
for the yaw and roll jets are given in Subsection 5. 3. 2. 2, Subsection 5. 3. 2. 3 con- 
sists of the controller coding summarized in flow-chart form. 
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5.3. 2.1 Aileron /Rudder Control 


Figure 5-8 is an overview of the aileron/rudder control logic, illustrating the 
two control submodes, with switching determined by the flag I j ^ . This flag is 
specified by vehicle controllability considerations (see Subsection 5.4) and is defined 
by a curve in the Mach-alpha plane as shown in Fig. 5-9. 


5. 3. 2. 1.1 Early Aileron/Rudder Control 

Figure 5- 10a illustrates the "reverse 11 control submode . When the dynamic 
pressure is less than aileron trim is set to zero, as are the aileron and rudder 
commands. When the dynamic pressure reaches q^, aileron trim and rudder com- 
mand are still set to zero, but closed-loop control of bank and sideslip (by the 
aileron alone) is provided by position (ft <t>) and rate (ft <j ) ) feedback of the lateral 
state variables, through a network which attempts to maintain a constant response 
throughout the trajectory. Specifically, the desired closed-loop transfer function 
for bank angle response is chosen to have four specified complex poles, or: 

k <d s ~ ^ > 

2 22 9 (5-4) 

desired (s + 2^^+ Uj ) (s + 2^u 9 s + w 4 ) 

where (£., <*>.) are chosen for acceptable transient response <k^ is not explicitly 
specified and to^ is a function of the airframe dynamics.). This response is ob- 
tained, in the absence of rudder effectiveness, by making use of sideslip mistrims 
to bank the vehicle. 



As with the longitudinal controller, the gains of Fig. 5- 10a may be separated 
into one of three types: (1) prespecified by the desired closed-loop characteristics 
(e.g., a t ): (2) dynamic pressure dependent: or (3) scheduled gains which are functions 
of the aerodynamic flight regime (e.g., C^). convenient reference, the gains 

of the first category are defined in terms of desired response characteristics with 
the following equation set: 


a 

a 

a 


1 

2 

3 


a 


4 


2( h U 2 + ^2“l )/ <“l“2 ) 
2 2 

to 4- + 4 it ft to Lo 

1 2 1 2 


(5-5) 


where, as noted above, the (£., co^) are chosen for acceptable roll response. 
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Fig. 5-8 Overview of Aileron/ Rudder Control Logic 
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Fig. 5-9 Control Submode Switch Parameter Definition 
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and it 


Aileron trim is commanded zero until the dynamic pressure reaches , 
is then determined through the use of a conventional clamped integrator in the bank 
angle error channel. This allows the vehicle to trim to a non-zero sideslip angle 
(by commanding a non-zero trim aileron deflection) in the presence of body-axis 
roll disturbance torques due to lateral displacements in the vehicle's center-of- 
gravity. To avoid interpreting large transients (such as occur during a bank rever- 
sal) as lateral cg-offset induced mistrims, integration of large magnitude errors 
is inhibited by the X<£ switch. The gradual expansion of the aileron command limits 
during this submode as dynamic pressure increases is shown in Figs. 5- 10a and 
5- 10b. 


PARAMETER 



DYNAMIC PRESSURE, q 


Fig. 5-10b Aileron/ Deflection Limit Parameter 


When the dynamic pressure reaches q g , later in this "reverse" submode, a 
rate command loop becomes effective when large attitude errors are detected ( as 
in a bank reversal). As shown in Fig. 5-10a, when the attitude error is greater 
than r^j, the logic is no longer concerned with maintaining a bank attitude, but in- 
stead is configured to drive the vehicle to a desired bank rate. This allows the 
vehicle to complete the maneuver at higher rates than normally obtained with the 
attitude control loop providing the feedback. 


Finally, when the Mach number drops below JJly the rudder is switched in to 
augment the aileron effectiveness. The aileron to rudder crossfeed gain is a 
constant . 
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5.3.2. 1.2 Late Aileron/Rudder Control 


Figure 5-11 illustrates the "proverse" control submode. Closed-loop control 
of bank and sideslip is provided by position (j3, <j>) and rate (fi, <b) feedback of the 
lateral state variables, through a network which attempts to maintain a constant re- 
sponse throughout the remainder of the trajectory. Specifically, the desired closed- 
loop transfer matrix describing bank and sideslip attitude response to a guidance 
commanded step is chosen to be: 


<S 2 + a u s + bjj) 

< a 12 s + b t 2^ 


P 


0 

(a 2 i s + b 21 ) 

< s2 + a 22 S+ b 22 ) 


* - *G 


0 


where the a^. and b.^ are chosen for acceptable transient response. It may be noted 
that in this proverse" submode, the availability of the rudder allows for con- 
siderably more freedom in shaping the transient response, so that both sideslip and 
bank characteristics may be specified, in contradistinction to the previous submode 
(see Eq (5-4)). 


As shown in the block diagram, there are eight gains (r.. and p..) which are 
derived functions of both the desired closed-loop parameters (in Eq(5-6)) and of 
four selected scheduled gains. They are given by: 


r l t a l 1 7 1 1 + a 21 G 1 7 12 


P 11 b ll 7 U + b 21 C l 7 1 2 


r 12 a 12 7 11 + a 22°l 7 1 2 


"21 a ll 7 21 + a 21 G 1 7 22 


P 12 b 12 7 11 + b 22 G 1 7 1 2 


P 2 1 b ll 7 21 + b 2 1 G 1 7 22 


(5- 7a) 


r 22 a 12 7 21 + a 22 G 1 T 22 


P 22 b 12 7 21 + b 22 G 1 7 22 


where G^ is a constant and the 7^ are scheduled gains. Also shown in the block 
diagram are a constant gain, G^, two additional scheduled gains and v^> and one 


additional derived gain, a function of the 7 . .: 
& ij 


* = 


7 11 7 22 ' 7 12 7 21 


(5- 7b) 


5. 3. 2. 2 Yaw and Roll ACPS Control 

Yaw and roll jet activity is completely inhibited while in the ’'proverse" aileron/ 
rudder control submode (I^ T = 1). This section describes the lateral ACPS 
control logic while in the "reverse" surface control submode. 

Yaw jet ACPS control synthesis is based on double integrator modeling of the 
bank angle dynamics, so that the switch logic presented for pitch ACPS control may 
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be used. Roll jet logic consists of three models, one based on double integrator 
modeling of the sideslip dynamics, a second based on oscillator modeling, and a third 
which accounts for control axis cross-coupling. 


5. 3. 2. 2.1 Yaw Jet Control 

A double integrator model of the bank angle dynamics, combined with the 
assumption of turn coordination provided by the roll jets, allows for the use of the 
A CPS pitch jet control logic, with the appropriate redefinition of control parameters 
Specifically, yaw jet control is defined by the phase-plane switch logic of Fig. 5-4 
and Eq (5- 3a), where the coordinates are given by: 

(5- 8a) 

and the switch curve parameters are given by: 

(u, U) = (u z , U z ) (5 - 8b) 

(6 1' 6 2 ) = U2 ) ( 5 - 8c) 


(e 


r 


C 2 ) 


■ 


P n . sin a r 
lim t* 


The blending interface logic is discussed below. 


(5-8d) 


5 . 3 . 2 . 2 . 2 Roll Jet Control 

Roll jet control of vehicle sideslip is accomplished by the use of three control 
modes. Two of these use a double integrator model for sideslip: the remaining one 
uses an oscillator model. 

The purpose of the first mode is to provide turn coordination when the yaw jets 
are firing (u z * 0). This is accomplished by the phase-plane logic of Fig. 5-12, 
which is a generalized (0, 0) plane. The figure illustrates a typical state trajectory 
and the resultant limit cycle (with switch curve overshoots). The equations for the 
switch curve y. are: 

W : SB = 2 ~ - /3 2 j if s/5 > li m 

(5- 9a) 

= -2U“fd m ) 2 - S' J if SB < 8 m 
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where 


Ui = U z COS CUrj. 

IT = U sin a - U„ cos a 

2 x T z T (5- 9b) 

S = sgn(u z ) 

0 m = \/2U 2 6 1 

and where <5^ is defined below in Eq (5-10). 

The purpose of the second mode is to provide direct control over sideslip when 
the yaw jets are inactive (u = 0) and when the dynamic pressure is low. With double 
integrator modeling, the phase-plane logic of Fig. 5-4 and Eq (5- 3a) may be used, 
where the coordinates are given by: 

(e i, e 2 ) = (0, P ) (5" 10a ) 


and the switch curve parameters are given by: 

(u, U) = (u x , U x sin at T ) 

(6 r « 2 ) = « i32 ) 



(5- 10b) 
(5- 10c) 
(5-10d) 


Fig. 5-12 Roll Jet Control when u^ * 0 

It is convenient to label this control mode (used at low q and when u z = 0) as attitude 
hold, since both rate and attitude error are driven toward zero. 


5-23 


The third control mode takes advantage of the vehicle's dynamic stability in 
sideslip, which becomes significant once the dynamic pressure becomes sufficiently 
high. What is required is simply rate damping of the sideslip angular velocity as 
shown below, where 0 m is given by Eq (5- 9b). 


0 


- 

t 

■ U x 

j8 m 

u x 


-/3 m 

~ U x - 

- 




Fig. 5-13 Roll Jet Control when = 0 (q > q^) 

The control modes described above, defining lateral control with yaw and roll 
jets, are summarized in block diagram form in Fig. 5-14. Several points should be 
noted. First, the mode select block which defines the roll jet control mode to be 
used is defined by: 

If u z ^ 0 then Mode = I 

u z = °, q < q 3 II (5-11) 

u z = 0 q > q 3 III 

where the modes are coordination, attitude hold, and rate damping, respectively. 

The three roll control blocks correspond with the above modes. A second point to 
note is the gating by the blending parameter X which is defined by logic that is 
identical to that of the pitch channel. Finally, it may be noted that the logic dis- 
allows ACPS inhibition during aileron trim-only (null aileron command) operation. 

5 . 3 . 2 . 3 Digital Implementation 

Coding implementation for lateral control during the entry phase is summarized 
in flow-chart form in Fig. 5-15. 
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5.3.3 DFCS Parameter Specification 


All parameters used in control computations for the entry phase are either 
fixed or are scheduled on the basis of Mach Number, 7 7[. Table f>-l tabulates all 
entry- phase fixed control parameters. The scheduled control gains ^ through 
and 1 / and v 2 are plotted in Figs. 5-1G through 5-23. The scheduled control param- 
eters through y. ?2 are plotted in Figs. 5-24 through 5-2 7. As noted in Subsection 
5. 3. 2. 1.1, the lateral gains r„, p^, and ri are calculated as functions of the y.. # 
according to Eq (5-7). 

In an initialization pass through the entry control logic, several parameters 
which are functions of pad-loaded values are calculated and stored for future use. 

The Toad Manipulation Routine for the entry phase is shown in Fig. 5-28. 

5 . 4 P erformance Summary 

Shown in Fig. 5-29 are three curves in the Mach-alpha plane which define re- 
gions of lateral controllability using differing control modes. Also shown is a 
typical late entry Mach-alpha profile. The curves were derived by conducting simu- 
lated maneuvers with the entry DFCS and the NR 00089B (Feb. version) vehicle at 
various Mach numbers and angles-of-attack. The curves separate desirable from 
undesirable simulated vehicle response. 

Curve A indicates the lower boundary (in the Mach-alpha plane) to which the 
vehicle may be flown using ailerons alone for lateral control. By driving the rudder 
proportionally to the aileron command (with a crossfeed gain of 2), the region of 
"reverse" control may be considerably extended farther down the entry profile 
(curve B) . Finally, the region of controllability using the proverse control mode 
is bounded (above) by curve C. The resultant overlapping region, defined by curves 
B and C, provides the motivation for the switch curve choice shown in Fig. 5-9. 

It should be noted that an on-going effort in determining design sensitivity to 
aerodynamic coefficient uncertainties will reduce the area of the overlap region. 
However, the control design used in deriving the curves of Fig. 5-29 was not tuned 
so as to maximize the resultant overlap area. It is anticipated that future efforts 
in this direction will enlarge this region for the nominal aerodynamic coefficient 
values, while counteracting the shrinking effect of off-nominal aerodynamics. 
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Table 5-1: Fixed Parameters 


Symbol 

Value 

Dimension 

U 

X 

1. 55 

deg / sec^ 

U 

y 

1. 32 

deg/ sec 2 

u 

z 

0. 66 

deg/ sec 2 

6 *1 

0. 75 

deg 

6 «2 

1. 0 

deg 

6 i3 1 

0. 2 

deg 

6 /32 

0. 4 

deg 

%1 

0. 75 

deg 

6 *2 

1. 5 

deg 

a 

a 

5 

- 


2 

- 

% 

5 

- 

f i 

0. 8 

- 

f 2 

0. 6 

- 

f 3 

0. 6 

- 

f 4 

0. 8 

- 

q l 

10 

lbf/ft 2 

CM 

Icr 

20 

lbf/ft 2 

q 3 

5 

lbf/ft 2 

q l 

20 

lbf/ft 2 


30 

lbf/ft 2 

q 3 

80 

lbf/ft 2 
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Table 5-1: Fixed Parameters (Cont) 


Symbol 

Value 

Dimension 


6 

-1 

K «i 

- 1 . 0 

sec 

-1 

K a2 

0. 2 

sec 

-1 


0. 01 

sec 

-1 

w d 

0. 7 

sec 

-1 

W 1 

0. 7 

sec 

-1 


0. 7 

sec 

Cd 

0. 7 

- 

Cl 

1. 4 

- 

<2 

0. 7 

- 

r l 

5 

deg 

r 8 

25 

deg 

r 0 l 

5 

deg 

r 02 

25 

deg 

r 0‘3 

8 

deg/ sec 

C RUU 

2 

-1 

a ll 

1. o 

sec 

-1 

a 12 

0. 0 

sec 

-1 

a 21 

0. 0 

sec 

-1 

a 22 

1. 3 

sec 

-2 

b ll 

1. o 

sec 

-2 

b 1 2 

0 . 0 

sec 
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Table 5-1: Fixed Parameters (Cont) 


Symbol 

Value 

Dimension 

b 21 

1.0 

-2 

sec 

b 22 

0. 8 

-2 

sec 

G 1 

0. 133 

- 

G 2 

29. 0 

lbf-sec' 

. min . max 
6e , 6e 

see 


. max . max 
6a , 6r 

Subsection 2-2 
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Fig. 5-20 Gain Schedule for 



Fig. 5-21 Gain Schedule for 


MACH NO. 


MACH NO. 


5-37 













START 



Fig. 5-28 Entry Load Manipulation 
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Fig. 5-29 Controllability Regimes 



SECTION r> 


ON-ORBIT PHASE 


6.1 Introduction 

The DFCS (Digital Flight Control System) will control the attitude and the 
translation of the Space Shuttle during the various orbital activities by using either 
the A CPS (Attitude Control Propulsion System) or the OM S (Orbital Maneuvering 
System). For the purpose of control, there are two orbital phases: the ACPS phase 
when the lets are the control effectors, and the TVC (Thrust Vector Control) phase 
when the gim baled engines are used for control. The ACPS phase is presented first, 
in Subsection 6.2, followed by a preliminary description of the on-orbit TVC phase in 
Subsection 6.3. Digital implementation of the ACPS phase, presented in Subsection 
6.4, and preliminary testing was accomplished using structured programming techniques 
in the HAL language. All symbols relevant to the on-orbit phase may be found in a 
separate glossary at the end of this section. 

6.2 On-orbit ACPS Phase: Engineering Description 

A block diagram of the DFCS during the on-orbit ACPS phase is 
illustrated in Fig. 6-1. The on-orbit attitude maneuver routine or another inter- 
face routine (see Part II) supplies to the DFCS: 

1) a set of reference gimbal angles, [0j] 

2) a set of incremental reference gimbal angles, [A 9^ ] 

3) the desired angular velocity vector, 

4) a set of bias gimbal angles, [ AS ] 

5) a bias acceleration vector, 

Every t f seconds the DFCS adds [A0j] to [0j] to give a new set of reference gimbal 
angles. The DFCS then compares these to the measured set of gimbal angles. [# GA ], 
and computes the angular errors in control axes. The bias angles are added td these 
errors in control axes and fed to the phase plane switching logic. Also fed to the 
phase plane switching logic are the angular velocity error. U e> the net estimated 
disturbance acceleration, a, , the magnitudes of control accelerations, [<* c ]. and 
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Fig. 6-1. On-Orbit ACPS 













the last desired change in angular velocity. The angular velocity error, is com- 
puted by subtracting the desired angular velocity, from the estimated angular 

velocity, w, and transforming the result into control coordinates. The net estimated 
disturbance acceleration, o\j , is computed by adding the bias acceleration, to 

the estimated disturbance acceleration, Sj, and transforming the result into control 
coordinates. From these inputs the phase plane switching logic determines the 
desired angular rate changes, [Aw] , in control coordinates. Prom these desiied 
rate changes the jet selection logic determines the appropriate jets to be fired, if 
any, and the corresponding firing times, and then turns the jets on and schedules 
their turn-off . 

6.2.1 ACPS Control Axes 

This subsection presents a technique for decoupling the ACPS control of the Shuttle 
Orbiter into three essentially independent axes. This technique permits one to treat 
the rigid body rotational dynamics of the vehicle, normally represented by a 6 order 
system, as three independent 2 nd order systems, thus justifying the use of phase 
plane control logic. A similar technique was used in the design of the Lunar Module 
autopilot. The jet selection logic which has been implemented to utilize the de- 
coupled axes is limited to specified jet pairs and provides only attitude control, not transla- 
tion, in its present form. A more general scheme will be used in the future. 

The ACPS thruster group locations for the NR ATP baseline configuration are 
given in Subsection 2.2. Define a +u thruster combination as one or more thrusters 
from group 2 and/or group 12: e.g. , jet #14 and jet £31 . Also define a -u thruster 
combination as one or more thrusters from groups 1 and 11: e.g., jets £10 and 2a. 

The only restriction in selecting the jets is that the torque produced by the +u com- 
bination is opposite (or nearly opposite) in direction to that of the -u jets. Similarly, 
select a +v combination from groups 4 and 10, e.g., £9 and # 26; and a -v combina- 
tion from groups 3 and 9, e.g., #13 and #32. Finally, select a +w combination from 
groups 6 and 8, e.g. , #5 and #18; and a -w combination from groups 5 and 7, e.g., 

#6 and #17. The suggested thruster combinations give moderate rotational control 
authority while minimizing the net translation. 

The control axes are defined by the accelerations produced by the u, v and w 
jets. More specifically, let 

-BASIS = Unit vector in the direction of the +u axis 
= unit vector in the direction of the +v axis 

— BASIS 

— B'XSIS = Unit vector in the directiori the +w axis 
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J 


orbiter inertia matrix 


— P2U’ — M2U* — P2V’ control torques produced by the +u, +v, 

— M2V' — P2W' — M?W - W t ^ iruster combinations 


then 


-BASIS 

- unit 

a" 1 

( -P2U ‘ 

5?M2U> /21 

-BASIS 

= unit 

(j -1 

<Mp2V ' 

5?M2V» /2 > 



*- 1 


— M2W* 

-BASIS 

= unit 

(j 

^ — P2 W ‘ 


(The components are given in vehicle axes.) 

The control axes orientations are shown in Fig. 6-2 for the thruster 
combinations given above. Note that the +u accelerations are not exactly opposite. 
Neither are the +v accelerations. 

The transformation from control axes to vehicle (or pilot) axes is defined as: 

s}c 

-CP ' ( -BASIS' -BASIS' -BASIS* 

The transformation from vehicle (or pilot) axes to control axes is simply: 

M pc = (T&cp )* 1 

The attitude and rate errors for the autopilot are derived in vehicle axes. In 
order to utilize the decoupling technique these must be transformed from vehicle 
axes to control axes, i.e. , 



6.2.2 State Estimator for the ACPS On-Orbit Phase 

Control of the shuttle vehicle in orbit requires estimates of vehicle angular 
rate, w, and vehicle angular disturbance acceleration, . These quantities are 
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x veh 



^veh 



Fig. 6-2. u, v, w Control Axes 
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determined from measurements of the IMU gimbal angles. The procedure used 
is based on Kalman filter theory in which angular measurements are mixed with 
extrapolated state estimates to form a filtered estimate of the state. The process of 
extrapolation and measurement incorporation is described first in terms of a single 
axis. The complete three axis filter as implemented in the DECS is discussed at 
the end of this subsection. 


6.2.2. 1 Vehicle Model 


In estimation theory it is necessary to define a model of the process being 
observed. In the case of the shuttle vehicle, we will assume that it is a rigid body 
with the following equations of motion for one axis: 


where 


V 


0 

1 

— 

0 


Y 


0 


0 


0 

CO 

= 

0 

0 

1 


CO 

+ 

O' 

+ 

Of 

+ 

n 









'C 


m 


CO 

*d 


0 

0 

0 


°d 


0 


0 


n 

a 

— — 


— 


— 


— — 


— _ 


_ _ 


L — 
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vehicle altitude 


- vehicle angular velocity 
= angular disturbance acceleration 

o ~ control acceleration (produced by ACPS jets) 


n . n 

a 


known bias acceleration 
driving noise terms 


6. 2. 2. 2 Extrapolation 


The discrete filtering process can be divided into two parts: extrapolation 
and measurement incorporation. In the extrapolation process the estimate of the 
vehicle state, _$(tg), given at some time, is integrated forward to the next meas- 
urement incorporation time, t ^ . This integration includes the effect of any known 
biases or control inputs during the interval t^, t^. The result of the extrapolation 
process is a predicted estimate . In our case x'O^) is determined by integrating 

Eq . (6-1), i.e., 


A 

x 


* 

:'( t x ) = 4>(t r t 



r 


T) 


0 

°C (T) 

0 


dT 


+ 



0 


a ( t ) 
m 


d t 


0 


(6-2) 
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where 


and 


*(tj. T) = 


£(t 0 ) 


W 


“V 


“dV 


= state estimate at t 


0 


(tl - T ) 
1 
0 



= transition matrix 


(6-3) 


(6-4) 


a ,(r) is the control acceleration applied by the ACPS jets. A typical plot of a Jr ) 


C 


is shown in Fig. 6-3. It is assumed that the control acceleration can have one 
of three values I a I, 0, - ItfJ . It is also assumed that if a jet goes on, it turns 
at the beginning of the interval, t^. It is further assumed that when the jet goes off 
it remains off until the end of the interval. With this jet model we can now integrate 


1 * 

j* #(t lf t) 
l 0 


0 

“c ( t) 

dr 

|a c | SIGN(a c ) t JET (AT - ^ - -) 

|o c | SIGN(a c ) t JET 

0 


0 


where 


L 

AT = t. 


*0 ^ast 


(6-5) 


OL (t) is a known bias acceleration which is assumed to be constant over the 
m 

over the integration interval: 


“m (T) = “mV 


The third term in Eq. (6-2) then reduces to 
Before grouping terms, define 


0 “ 


“mV* T (A f- r 

O / \ 
“m T 

d t = 

“mV * T 

. o . 


0 


( 6 - 6 ) 


“WV s t W + V'o’i 1T 

Aco^(Iq) - | | SIGN ( t q ) ) 
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Combining results, we find 


A A A AT C JET 

e B (t l ) = ®B (t 0 )+ w(t 0 )AT 4 Aa) dm <t 0 ) V + Aw C (t 0 MAT 2— ' 


(6-7a) 


w'U^ = U(t 0 ) + ^w dm (t 0 ) 4 Ao^(t 0 ) 


“d (t l ) = “dv 


(6-7b) 

(6-7c) 


6. 2. 2. 3 Measurement Incorporation 

For the measurement incorporation process, the extrapolated state is combined 
with the measurement at time t in the following manner: 


W = ®B V + K e “ W ] 

«(t n ) = «'<‘ n > 4 K w ( ® (t n ) - ®B (t n )) /(t n-V 
“d (t n ) = “d (t n> + K a ( ® (t n ) ' °B (t n» 


(6- 8a) 
(6-8b) 
(6-8c) 


where 


t, -t 
Kj n 


= position gain (constant) 

K = velocity gain (constant) 

CO 

K = acceleration gain (constant) 

a 

0 (t ) = angle measurement at time t 

n n 

t = time of the last measurement incorporation 


This measurement incorporation process is similar to the linear Kalman filter 
if the measurement interval, t n - t^, is fixed. However, for the Shuttle filter, meas- 
urements will be incorporated only if 

I ® <t n ) ‘ ®B <t n ) I > 0 BT 


where 


®BT = " trap " an S le 


(6-9) 


If this condition is not true, then no measurement is processed and 

vv ■ w 

U)(t n ) = w'(t n ) 

a,(t ) = aUt ) 
d n d n 


(6- 10a) 
(6-1 Ob) 
(6-10c) 
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This additional logic helps to reduce the effect of quantization of the gimbal angle 
measurements by "trapping" small-amplitude noise before it reaches the estimate. 

6 . 2 . 2 . 4 The 3-Axis Filter as Implemented in the DFCS 

The 3-axis state estimator is basically a simple extension of the single axis 
filter described above. Due to the kinematics of Euler rotations, however, the gim 
bal angle measurements are not simply related to the vehicle body rates. We can 


circumvent this problem 
angles are defined as the 

by defining a set of "body angle" measurements, 
integral of the body rates: 

The body 


®B 1 ” J" “l dt 

(6- 1 1 a) 


®B 2 = ! w 2 dt 

(6-llb) 


9 b 3 = J* “ 3 dt 

(6- lie) 


To a good approximation we can calculate the body angles in recursive fashion as 


(n) 

1 


e D (n-l) 
B 1 


1 


0 


e GA 1 (n> © e GA 1 (n_1) 

(n) 

2 

= 

e R (n-l) 
B 2 

© 

0 

c°s<e aA >cos<e Gfl3 > 

^'•ga, 1 


e GA 2 (n) © 9 GA 2 (n ' 1) 

(n) 

3 


0T3 (n-l) 

B 3 

1 

0 

- S1N < e GA 1 )COS<9 GA 3 ) 

C°S(«g A i > 


9 GA 3 (n) © 9 GA 3 (n ' 5 ' 




( 6 - 12 ) 

where 


and 


®GA 1 

®GA 


outer gimbal angle 
inner gimbal angle 
middle gimbal angle 


© indicates modular addition 
© indicates modular subtraction 


(Modular addition and subtraction of two angles produces a result which lies in the 
range: -180° < 6 £ 180°. ) Measurement incorporation is performed as before (Eq. (6-8)) 
except that the measurements are now the computed body angles (Eq.(6-12)). For the 
x axis. 


if 


(t ) 
B n 
x 


8 b VI 

X 


> e 


BT 


(6-13) 


6-10 



then 

CD > 

a 

X 

= ®B V 

X 

© K 6[®B x (t n ) © 0 B x ( V] 

(B- 13a) 


“x (t n ) 

+ 

< 3 

n 

k PeL 

u>L B x 

V © hj 

(B- 13b) 


a dx (t n ) 

= a A (t ) + 
dx n 

K aPe 

L x 

V © »n - U ' 2 

X X 

(6-1 3c) 

otherwise , 

l f 

X 

-■ t 

n 

X 

= 9 b (t n> 

X 

(6-1 3d) 
(6- 14a) 




u (t ) 
x n 

= 

(6- 14b) 




a A (t ) 
dx n 

0i dx (t n ) 

(6- 1.4c) 


In the above 


equations 

0 = estimated value of the x body angle 

K 

x 

0 B ' ■= extrapolated value of the x body angle 

x 

~ - "measured’ 1 value of the x body angle 

0 B 



estimated x-axis body rate 


= extrapolated x-axis body rate 

a, = estimated x-axis component of disturbing acceleration 
dx 

a! = extrapolated x-axis component of disturbing acceleration 


t - current time 
n 


t 


t 


x 


time of the last measurement incorporation for the x-axis 


The y and z axis computations are similar. Note that a separate “trap" is required 
for each axis. 

After the appropriate jet control is exercised the extrapolation of the vehicle 
state, which depends to a large extent on the ACPS jet selection and timing logic used 
in the Shuttle autopilot, is performed. It is currently assumed that the jet firings 
produce accelerations about the u, v and w axes defined in Subsection 6.2.1. That is, 
the angular acceleration produced about the vehicle axes by the jets is the sum of the 
three components: 
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(6 -15a) 


-CU 


J ( -P2U ' — M2U^ 2 


^CV 


-CW 


•I ^P2V ~ — M2\T 2 


J ( -P2W -M2W J/2 


(6- 15b) 
(6- 15c) 


where J is the inertia matrix and the M's are the torque vectors defined in Subsection 
6.2. 1. It is also assumed that the restrictions on the timing of jet firings which were 
given for the single-axis filter apply. If we define 

t jet = time of firing of the u jets during the interval t Q < t < tj 

t jet = time of firing of the v jets during the interval t Q < t <t ^ 

t jet = time of firing of the w jets during the interval t^ < t < t^ 

and 

SGNC^ = sign of the desired u control acceleration 

SGNC 2 = sign of the desired v control acceleration 

SGNC^ = sign of the desired w control acceleration 

then by analogy with the results of the single axis filter: 



^uV = 

•-cir 

SGNC 1 

t JET 1 

(6- 16a) 


Aw v n 0 ) = 

l-cv 1 

SGNC 2 

“jet 2 

(6- 16b) 


A «w (t 0> = 

IScwl 

sgnc 3 

‘jETj 

(6-1 6c) 

and 







^ — dm (t 0^ = 

*-d^0^ + -m^0^ T 

(6-1 6d) 

Note that these are vector quantities 
the attitude maneuver routine as 

. The bias acceleration, a 

— m 

is computed by 



OL 

— m 

* -J ’ 


(6-17) 


where o) d is the desired maneuver rate. With these definitions, the complete three 


axis extrapolation equations can be written as 


-B !b ( V © 


tt(t Q )AT + iu,„(t n ) 


AT 


dm 0 


JET- JEj 1 n 

+ Au^ (t 0 ) (AT 2 ) + Aw v <t 0 ) (AT g— — ) 


JET, 


AffiwV <6T 


l JETo 

5 


(6- 18a) 
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(6-1 8b) 


w'Uj) = w(t 0 ) + ^w dm (t 0 ) + Aw u (t 0 ) + 6co v (t 0 ) + Aw w (t 0 ) 

a^(tj) = oj^Uq) (fi-l8c) 

Note that these quantities are vectors with components in the x, y, z vehicle (body) 
axes. 

The estimate of disturbance acceleration, , may be accurate enough for the 
powered flight phases, but it may be too noisy for use during drifting flight. If an 
estimate is needed during coasting flight to counteract such disturbances as gravity 
gradient and aerodynamic torques, then rate gyro measurements may be required. 

If not, the 0^ estimation process may be disabled by using appropriate filter gains. 
Alternatively, the gain, simply can be reduced during coasting flight. In any 

event, the bias acceleration estimate, a , should be used during coasting flight 
maneuvers to reduce fuel consumption. 

6.2.3 Control Laws for the On-Orbit ACPS Phase 

The object of the phase plane switching logic is to compute changes in angular 
velocity so that by turning on or off the appropriate jets to achieve this velocity change 
the vehicle limit cycles about a desired attitude and attitude rate. The two dimensional 
phase plane is applied to the decoupled control axes described in Subsection P.2.1. Limit 
cycling is desirable since it would be extremely costly to attempt to maintain the de- 
sired attitude and rate exactly. By designing appropriate switching lines into the 
phase plane, the vehicle can be driven into the desired limit cycle. Two separate 
logics are used depending on the size of the position and rate errors. For large er- 
rors, a switching logic has been designed to drive the state into the small-error con- 
trol region. The small-error switching logic then drives the state into a limit cycle 
about the desired attitude and attitude rate. The distinction between large and small 
errors will be quantified later in this section. 

The small-error switching logic has the following features: 

1) The design is based upon the optimal switch curves for minimizing a 
weighted combination of fuel and response time. The optimal curves 
are separated by a deadband in order to allow for limit cycling. 

2) The appropriate deadband can be chosen in order to 

a) fix the limit cycle period 

b) fix the limit cycle maximum rate 

c) fix the maximum attitude excursion while limit cycling 

3) The average attitude error is minimized by biasing the placement of 
the switch curves. 
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4) Parabolic switch curves approximate the trajectories which the state 
follows. However, hysteresis logic prevents performance degradation 
due to errors in these approximations. 

The design incorporating these features is described in the following subsections. 

In Subsection 6. 2.3.1 two idealized limit cycles are presented. One limit cycle is 
in the presence of positive disturbing accelerations while the other is in the presence 
of zero disturbing accelerations. For each limit cycle the average attitude error will 
be derived. In Subsection 6. 2.3. 2 the design of switching curves to drive the state into a 
limit cycle with zero average attitude error is described. This design is referred to 
as the small-error control logic. The large-error control logic is described in 
Subsection 6 . 2 . 3 . 3 . 

6. 2. 3.1 Idealized Limit Cycles 

An idealized limit cycle for a positive disturbing acceleration is illustrated in 
Fig. 6-4. The coordinates 8 e and tO e are chosen to be the state, where 0^ represents 
the angular position error and ciL represents the angular rate error. X,. ri is the 
attitude variation in the course of limit cycling. In the presence of a positive disturb- 
ance acceleration, , the state is driven along the trajectory from A to B. At point 
B a negative control acceleration, -U^, is applied to return the state to point A. 


e 

e 



Fig. 6-4 Idealized Limit Cycle for a Positive 
Disturbing Acceleration 
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The average attitude error, X, of the limit cycle in Fig. 6-4 is found by deter- 
mining the average attitude error for each of the trajectories OB and BO' (called 
and X ) and then 'weighting” them according to time to calculate the net average, X. 

Determination of Xp 

X. is, by definition, 

h 

x = i / ? u D t2dt (6 - 19) 

1 0 


tp the time to travel from point O to point B, can be calculated from the value of the 


9 abscissa of point B, B a , since 
e e 


B e = ? u d 1 i 2 

e 


( 6 - 20 ) 


Bg is found by first rearranging 0^ = f^u^) and 0^ - ^2^e^ suc ^ that, 


w e = W e = g 2 (0 e ) 


Then, at point B, g^©^) ~ ^2^®e^ anc1. 


and from Eq(6~20) 


Using t^ in Eq(6-10) 


B 0 ' X DB 
e 




+/ 2X DB (U C " U D ) 




X 1 X ( U C - U D ) 

X I 1 X DB \ U c ) 


( 6 - 21 ) 


(6-22) 


(6-23) 


(6-24) 


Determination of Xg: 

> 

The time average attitude error for the trajectory segment BO' can be found by 
setting t = 0, at point B. Then, 

t 2 

*2 ■ r ! [X o + V + 1 ®D - u c> ‘ 2| d > 

2 0 

o r, X t 

S 2 - ^ <U D • V c n 2 2 (6-25) 

X o is defined in Eq (6-22) and X Q can be determined from. 
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V Vi 


+ / 2X DB U D < U C " U D* 

J 


(6-26) 


To find 1 2 one can use the following relation, 


X o + < U D - V‘2 ' 0 


from which it follows. 


/ 2x l> b u d 


C ' C U D' 


(6-27) 


Substituting the results of Eq (6-22), (6-26) and (6-27) into Eq (6-25), yields: 


X 2 “ X Db( 1 " 3U / 


(6-28) 


Determination of X: 

Finally, the total average attitude error, X, can be found from: 


_ ^ Xjtj + X 2 t 2 

x " t 4 t 
C 1 C 2 


(6-29) 


Substituting the results of Eqs (6-23), (6-24), (6-27), (6-28) into Eq (6-29), yields: 


(6-30) 


This is the average attitude error for the limit cycle in Fig. 6-4. Thus, the average 
attitude error would be zero if the intercept of the curve on the right, which is the 
trajectory along which the control is applied, were shifted from point O' to have an 


abscissa of 


DR = X DB - X 


(6-31) 


A similar result can be obtained for the case of a negative disturbance acceleration. 
Equations (6-30) and (6-31) are used in the DFCS in order to minimize the average 
attitude error. 


The attitude variation, X T -. T3 , could be chosen as follows: 


1) To achieve the cycle period, T: 


set X, 


D 1 t 2 


(6-32) 


2) To achieve the limit cycle maximum rate, X 


set X, 


u r . (x r 

1 C 1 max 

2|u D | (|U C | - |U D |) 


(6-33) 


3) By simply setting X^g to the allowable attitude variation. 
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Currently, the final option is employed in the DFCS. 

An idealized limit cycle in the presence of no disturbing accelerations is illus- 
trated in Fig. 6-5. From point A the state is allowed to drift 



k 

A 

B 

f » ~ 

1 . 

• di -v 

r R 

D 

C 


Fig. 6-5 Limit Cycle in Presence of No 
Disturbing Accelerations 


to point B where a negative control acceleration, -U^, (in the form of a minimum 
impulse jet firing) is applied to drive the state to point C. From point C the state is 
then allowed to drift to point D. At point D a positive control acceleration, -*U C , is 
applied to return the state to point A. By choosing DL - DR the state limit cycles 
with zero average attitude error. The generalized limit cycle for positive disturbing 
accelerations is made to reduce to the limit cycle for zero disturbing accelerations 
by selecting 


DL = X + 7 X DB 


(6-34) 


6. 2. 3. 2 Small-Error Control Law 

Control regions for relatively small errors are illustrated in Fig. 6-6, defitjk 
tions of curves are summarized in Table 6-1, and corresponding control actions are 
given in Table 6-2. All logic in this subsection is formulated in terms of positive dis- 
turbing accelerations. In the actual implementation, (see Subsection 6.4), negative 
disturbances are treated by inverting the signs of the state variables © e and u> e . 
Biasing the right and left deadbands, DR and DL, to null the average attitude error 
in accordance with the equations derived in the preceeding subsection, defines the 
placement of the switching curves BOUNDARY-II -TO-jCOAST and BOUNDARY-IIIJTO_ 
COAST. These curves approximate the trajectory of the state during a firing. When 
the state lies in Region II, control is applied to drive the state to the curve TARGET- 
DISTURB. The vehicle is then allowed to coast. 
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BOUNDARY II TO COAST -( X + DR ) 



Fig. 6-6. SMALLERROR DISTURBANCE LOGIC 





Table 6-1: Definitions of Boundaries of Small-Error Control Regions 


Boundary / Parameter 
BOUNDARY. I_TO_Il 

BOUNDARY- iI_TO_COAST 

BOUND A R Y_ I „TO_COA ST 

BO UN DA R Y_III_TO_COA ST 

BOUNDAR Y_IV_TO_III 

BOUNDA R Y-IV-TO -COAST 


TARGET- DISTURB 


WHERE 


‘ U b- 0 - 2x ub 


DR = X DB - X 


Definition 


to 


0 

‘‘ 2(|U C | - |U D |) 

2 

8 

e 2<|U C | - |U U |) 
2 


0 = + - 

e 2CD(|U | 4 | U | ) 


^ DR ^ 


4 DR 


4 DR / 


C 1 1 D 1 

2 

U! 


- DL 


- DL ' 


9e 2(|U C | + |U U |) 

2 

W e 

0 - 4 - — 

“ 2(|U C | 4 | U (J | ) 

2 

CO 

Q = ? - DL' 

e 2CD( | U c | - |U D |) 

2 

to 


0 = 4 

e 


- X 


2 U 


D* 


DL . X + j X DB 


DR ' = DR 
DR" = (X 


+ a db i: 

+ DR') (|U D | ) 


' = DL + 4 db 

(lU c l - |U D |) + CD(|U C | 

(CD(|U C | + |U D |> - |U D | 


- IUdP 


) < |U r | - 
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Table i>-2; Control Actions in Small-Error Control Regions 


Region 

I 

II 

III 

IV 
HS_I 


HS-II 


CS..I 

CS..II 


Control Action 

Drive state to BOUNDA RY_ I__TO_COAST 
Drive state to TARGET- DISTURB 
Drive state to zero rate 
Drive state to BOUNDARY- 1 V_ TO- COAST 

If coming from LARGE ERROR CONTROL 
LOGIC or REGION III then drive state to 
zero rate, else coast 

If coming from LARGE ERROR CONTROL 
LOGIC or REGION II then drive state to 
TARGET- DISTURB, else coast 

coast 

coast 


The deadband increment A^g is added to prevent fuel-consuming overshoot when 
the disturbance acceleration is zero or very small. The points 

DR ' = DR 4 A db (6-35) 

DL ' - DL 4 (S-36) 

define the placement of the curves BOUNDAR Y_ L_TO_COAST and BOUNDAR Y__IV_ TO_ 
COAST. 

The factor CD determines the weight given to response time as opposed to fuel 
consumption, by adjusting the steepness of the switch curves that are normally crossed 
when entering the coast region. To insure that BOUNDARY_IV_TO_JCOAST is always 
less than BOUNDARY- 1 1_TO_COA ST and BOUNDARY-IJTO— COAST is always greater 
than BOUNDARY— III— TO— COAST, CD must be between 0 and 1. Currently the value 
of CD is 0. 8. 

Hysteresis has been incorporated to make the most efficient use of disturbance 
accelerations, to prevent switching line chatter, and to prevent tight limit cycling about 
the point DR. The dashed curves within the shaded area and parallel to BOUNDARY 
II - TO - COAST P r °vide hysteresis in the regions HS_I and HS_II. If the state lies 
inside either Region HS_I or HS_ II and the sign of the newly computed change in rate 
is not the same as the sign of the last computed change in rate, then no control is ap- 
plied and the vehicle is allowed to coast. However if the signs are the same, then in 
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Region HS_I the state is driven to zero and in Region HS-II the state is driven to 
TARGET-DISTURB. Whenever the state lies inside either Region CSM or CS_II the 
vehicle is simply allowed to coast. 

When the state lies to the right of BOUNDARY_I-TO_II position error is reduced 
more quickly by driving the state to BOUNDARY-I-TO-COAST and then coasting instead 
of driving the state to TA RGET—DISTURB and then coasting. This is particularly im- 
portant as the disturbance acceleration approaches zero. The curve BOUND ARY- 1. .TO. 
II is defined by the intersection of TARGET— DISTURB and BOUNDAR Y_I_TO_COAST. 

When in Region III the rate error is driven to zero, and when in Region IV the 
state is driven to BOUNDA RY_ IV—TO-COAST to more quickly reduce position and rate 
errors. 

6. 2. 3. 3 Large-Error Control Law 

Control regions for relatively large errors are illustrated in Fig. 6- 7, defini- 
tions of the boundaries are summarized in Table 6-3, and corresponding control ac- 
tions are given in Table 6-4. As in the preceeding subsection, all logic is formulated 
in terms of positive disturbing accelerations. In the actual implementation (see Sub- 
section 6.4), negative disturbances are treated by inverting the signs of the state var- 
iables, 9 e and u) e . 

The purpose of the large-error control logic essentially is to drive the state into 
the small-error control region. Bounds between the small-error and large-error 
regions are set by the values of the parameters RL and RL' (currently 4.0 and 3.0 
deg/sec respectively). The rate error bounds are equal to ± RL. The left position 
error limit, -PLL, is defined by the intersection of BOUNDAR Y_IV_TO_COAST and 
COAST-IV— BOUNDARY. Similarly, the right position error limit, +PLR, is defined 
by the intersection of BOUNDARY_I_TO_COAST and COAST— I—BOUNDARY. The 
boundaries COAST-IV- BOUNDARY and COAST— I—BOUNDARY are rate errors equal 
to ± RL'. 

TARGET-ARL is the average of RL and RL'. TARGET- ARL and -(TARGET- 
ARL) are taken to be the target rate when the state is in Regions II and I respectively, 
i.e. , when the firing is against the disturbance acceleration. When the state is in 
Regions III and IV the targets are the nearest boundary; i.e. COAST_IV_ BOUNDA RY 
and -RL respectively. 

Hysteresis is incorporated into the large-error coast regions to prevent chatter 
about 4RL and COAST-I..BOUNDARY. If the state is driven into Region HL-I from 
Region II, then it is further driven to TARGET-ARL. If the state otherwise appears 
in Region HL_I the vehicle is allowed to coast. Similarly, if the state is driven into 
Region HL II from Region I then it is further driven to -(TARGET_ARL); otherwise 
the vehicle is allowed to coast. Should the state appear in either Region CL_I or 
Region CL_II the vehicle is simply allowed to coast. 
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REGION II 
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Fig. 6- 7. Large Error Control Region 


Table 6-3: Definitions of Boundaries for the Large-Error Control Regions 


Boundary / Parameter 


Definition 


COAST- IV— BOUNDARY 
COAST-I— BOUNDARY 
TA RGET-ARL 

PLR 


PLL 


% = RL ' 

W = -RL' 
e 

« = (HL + RL')/2 

e ' 


0 

e 


e 

e 


(RL') 2 

2C1)( | U c | + |U D |) 

(RL') 2 

2CD( | U c | - |U U |) 


+ DR ' 
- DL ' 


Table 6-4; Control Actions in Large-Error Control Regions 


REGION 

I 

II 

III 

IV 

IIL_I 

HL_II 

CL_I 

CL_II 


CONTROL ACTION 
Drive state to -(TARGET-A RL) 

Drive state to TARGET— ARL 

Drive state to COAST- IV- BOUNDARY 

Drive state to -RL 

If coming from region II then drive state 
to TARGET_ARL, else coast 

If coming from region I then drive state 
to -(TARGET_A RT J, else coast 

coast 

coast 
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A final feature of the small-error control logic can be described, now that 
TARGET_ARL has been introduced. If the disturbance acceleration is very large*, 
part of the curve TARGETDISTURB may lie below - (TA RGET_ARL). In this situa- 
tion a check is made to see if a rate change to drive the state to TARGE T__DISTURB 
is greater than a rate change to drive the state to -(TARGET_ARL). If it is, the 
state is driven to - (TARGET_ARLb which means less fuel consumption with the re- 
sult of finishing in a limit cycle about the origin. 

6 . 3 On-Orbit TVC Phase: Preliminary Functional Summary 

The Orbiter Maneuvering Subsystem (OMS) is used for major on-orbit maneuvers 
and to establish the deorbit entry trajectory. During powered flight the roll, pitch 
and yaw of the space-shuttle orbiter is achieved through deflection of the two gimbaled 
engines. Computation of gimbal servo commands in response to computed errors be- 
tween commanded and measured attitudes is the function of the on-orbit TVC (Thrust 
Vector Control) phase of the DFCS. The following functional summary is an adapta- 
tion of the Apollo Command Module TVC control; it is illustrated in Fig. 6-8. 

1) The steering-loop computations generate incremental attitude commands in 
inertial coordinates and transform them into body coordinates. 

2) The DFCS back-differences the gimbal angles to obtain incremental variations 
over each sampling interval. 

3) The gimbal increments over each sampling interval are then transformed into 
body coordinates and subtracted from the commanded increments which are 
generated by the steering program. 

4) The respective attitude errors are fed to the roll, pitch and yaw compensa- 
tion filters, whose outputs contribute to the commands to the engine gimbal 
servos. These signals bear the burden of controlling rapid transients, deter- 
mining the system bandwidth and the stability of the vehicle slosh and bend- 
ing modes. 

5) As shown in Fig. 6-8b, the total command signal to each engine-gimbal 
servo is made up of the component from the associated compensation filter 
plus another component from a Thrust Misalignment Correction (TMC) loop. 
The latter component serves to bias the total command so that a zero output 
from the compensation filter will cause the thrust vector to pass exactly 
through the center of gravity (eg) when there is no eg motion and no motion 
of the thrust vector relative to the commanded angle. 
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AUTOPILOT SAMPLING PERIOD 
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Fig. 6- 8b. Block Diagram of the Thrust 





6) The two major elements of the TMC loop are its summing register, which 
supplies the bias, and its low-pass filter, DELFILTER, which tracks the 
total command signal at autopilot sampling frequencies. The difference 
between the bias and the DELFILTER output is used to increment the sum- 
ming register so as to slowly correct for the offset caused by thrust-to-cg 
misalignment. This slow reset action is roughly equivalent to inserting a 
proportional-plus- integral transfer function between the compensation fil- 
ter and the total command signal. The TMC loop is designed so that its 
dynamics have no effect on the vehicle bending and slosh modes and very 
little effect on the rigid-body stability. 

7) The operation of the TVC digital autopilot cannot be considered separately 

from that of the steering loop which interacts with the autopilot. The 
steering loop operations are shown in Fig. 6-8 for the case of External AV 
guidance. This form of guidance is based on achieving a desired acceler- 
ometer-measured velocity change, which is inserted into the steer- 

Lng program as an initial value of the computed velocity-to-be-gained vec- 
tor V . The V vector is periodically reduced during the course of the 

9 -g - g 

burn by subtracting the accelerometer-measured increments, AV. 

8) The purpose of the steering loop is to align the vehicle thrust with the 
velocity- to-be -gained. This is achieved by commanding a vehicle turning- 
rate which is proportional to the normalized cross-product between these 
two vectors. 

9) The pitch- and yaw-axis components of the vector cross-product are mul- 
tiplied by the proportionality constant K gteer to obtain the commanded 
rates about these axes - and multiplied by the autopilot sampling period, 

T 2 , to obtain attitude-command increments, which are supplied to the 
autopilot every T^ seconds. 

10) The roll control channel is designed to provide independent attitude and 

rate control about the roll axis by differential control of the OMS engines. 

6. 4 Digital Implementation (ACPS phase) 

The overall composition of the DFCS as it applies to the on-orbit phase is re- 
iterated in Fig. 6-9 in structured flow-chart form. The procedure MAJOR-TRANSIENT 
provides initial mode and flag information to be used by all mission phases. Currently, 
no initial mode or flag information unique to the on-orbit ACPS phase has been im- 
plemented in MAJOR-TRANSIENT. The procedure INPUT-INTERFACE contains, in 
general, information that is explicitly read by the DFCS from external sources. More 
specifically, when the on-orbit phase is selected, the gimbal angles are read and the 
time of their reading, t , is recorded. On the first pass of the DFCS for a given 
phase FILTEFLAND-PARAMMNIT and PARAM are both called. FILTER- AND— PA RAM — 
INIT currently does nothing in this phase. PARAM calls MASS_CALC (See Fig. 6-10) 
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Fig. 6-9. Logic Flow of DFCS 
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Fig. 6-10. MASS CALC 
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which then takes as an input the torques supplied by each ACPS thrust group and com- 
putes three more or less independent control axes and a pilot-to-control transformation 
matrix. The procedure PA RT1_ FILTER shown in Fig. 6-11 is then entered which 
takes as an input the extrapolated estimates of angular position, rate, and disturbance 
acceleration for each axis and incorporates the measurement of gimbal angles taken 
at time t n> From these updated estimates the procedure CONTROL (Fig. 6-12) calls 
upon: 

1) ERRORS (Fig. 6-13) to compute the errors in attitude and angular rate 
from the desired 

2) PHASE— PLANE (Fig. 6-14) to compute the desired change in angular rate 
to correct these errors 

3) JETXONTROL (Fig. 6-15) to select the appropriate jets and control 
their firings, according to provisional logic suggested in Subsection 6. 2. 1. 

The digital implementation of the phase plane switching logic is composed of 
three principle routines: 

PHASE-PLANE 

The procedure PHASE-PLANE selects the appropriate sign to be attached to the 
computed change in angular rate, computes position error intercepts, and also com- 
putes the boundaries between the small error control logic and the large error control 
logic, as illustrated in the logic flow diagram in Fig. 6- 14a. PHASE-PLANE then 
selects either the small-error control logic or the large-error control logic to com- 
pute the desired change in rate. 

LARGE-ER ROR-CONTROL-LOGIC 

The procedure LARGE-TERROR— CONTROL— LOGIC computes the change in rate 
to drive the state into the small error control region. As illustrated in Fig. 6- 14b, 
the region in which the state lies is singled out and then the appropriate control action 
is computed. 

SMALL-ER ROR-CONTROL-LOGIC 

The procedure SMALL-ERROR— CONTROL-LOGIC computes the change in rate 
to drive the state into a limit cycle with specified deadband about the desired position 
and rate. Current testing has used deadbands of 0.5 and 5.0 degrees. As illustrated 
in Figs. 6-14c and 6-14d, the region in which the state lies is singled out and then the 
appropriate control action is computed. 

The function BOUNDARY, illustrated in Fig. 6-14e, is implemented to compute the 
various switching curves. 
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Fig. 6-12. On-orbit Segment of CONTROL Procedure 
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Fig. 6-13. ERRORS 
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[F e * + PLR AND Ui * COAST-UBOUNDARY 


IF 0 * +PLR AND Ui * +RL 

e e 


IF 0 < -PLL AND u; < C OAST- 1 V_ BOUNDARY") 


IF 0 * -PLL AND ui * -RL 



THEN _ (COAST_I„ BOUNDARY - RL) 

r x » *«e 5 


IF a i <0.0 


IF X Dold !' 0 -° 


(COAST_IV_BOUNDARY + RL 


IF SIGNIK X D )= SIGN <X Uold >, 



Fig. 6- 14b LARGE ERROR CONTROL LOGIC 






1 



* 


A ^ 

*n * (i - c 

Fig. 6 - 14c SMALL ERROR CONTROL LOGIC 
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Fig. 6-14d. REGlON_II_TARGETTING 


B_SIGN W 2 

2 CD (UC + B_SIGN UD) + D 



Fig. 6- 14e BOUNDARY: FUNCTION (B_SIGN, a), D, UC, UD, CD) 
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Fig. 6-15a. JET_CONTROL Procedure 




Inputs: Qq, U d 


Outputs » ^JETQ* t jET > ^ 



Fig. 6-15b. UVWTIME Procedure 







IF KILL1 = ON 


THEN 


RETURN 



Fig. 6- 1 5d . OFF2UJETS Procedure 
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IF KILL1 = ON 


THEN 





Fig. 6-15e. OfflVJETS Procedure 



Fig. 6-15f. OFF2VJETS Procedure 
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Table 6-5: Glossary for On-Orbit Phase 
MNEMONIC ENGINEERING DEFINITION 



TERM 


ACC2U 

-cu 

Angular acceleration about u-axis produced by Mp^ T or M^ ?TT 

ACC2V 

IP 

o 

< 

Angular acceleration about v-axis produced by M p2v or M^ 2 y 

ACC2W 

-CW 

Angular acceleration about w-axis produced by M^ n , Tr or M™ owr 

— rZW — M l W 

AL PH AD BAR 

Sd 

Extrapolated disturbance acceleration 

ALPHADHAT 

2d • 

Estimated disturbance acceleration 

ALPHAM 


Bias acceleration 

basTsu 

-BASIS ] 

1 

basTsv 

-BASIS 

f 

Unit vectors along u, v, w control axes 

basTsw 

1 

—BASIS 

i 

[BIAS-BETA] 

^BIAS^ 

Bias body angles 

CD 

CD 

Scale factor used in phase plane switching logic 

[ CDU ] 

^ ®GA - 1 

Gimbal angles 

[ CDUER ROR] 

^GA^ 

Error in gimbal angles from the reference desired 

DELTA 

& DB 

Deadband factor 

DELTA-POSITION 

0 ' 

e 

Angular position error (within PHASE_PLANE) 

DELTA-RATE 

t 

“e 

Angular rate error (within PHASE_PLANE) 

[ DELTA _THETA_ 




E] [60^j] Incremental desired reference gimbal angles 

DELWDM 6w dm Change in angular velocity vector due to disturbance and bias 

accelerations 


DELWU Change in angular rate due to u-control jets 

DELWV Aw v Change in angular rate due to v-control jets 

DELWW Ao). w Change in angular rate due to w-control jets 

DL, DL-PRIME1, DL, DL ' 

DR, DR_PRIME1, DR, DR ' 

DR-PRIME2 DR" 

[DRATE] [ ^ ex: 3 Components of changes of angular rates 
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Table 

6-5: Glossary for On-Orbit Phase (cont) 

MNEMONIC 

ENGINEERING 

DEFINITION 


TERM 


DTC 

DTC 


Total firing time of an ACPS control jet for a maneuver 

ISTRTR 

ISTRTR 

Restart required flag (external) 

PHASER 

PHASER 

Mission phase flag (external) 

[FAST] 

[MEDIUM] 

[FAST] 
[MEDIUM] j 

Fast, medium, and slow sampling periods for each 
mission phase 

[SLOW] 

[SLOW] 

Slow sampling period for each mission phase 

GAIR 

g air 


Filter gain for air data 

INERTIA 

J 


Inertia matrix 

IM 

IM 


medium rate sampling index 

IS 

IS 


Slow rate sampling index 

ISTART 

ISTART 

Initial pass flag for mission phase 

ITM 

ITM 


Ratio of medium to fast sampling periods 

ITS 

ITS 


Ratio of slow to fast sampling periods 

[ jtorque] 

[ .JTORQUE] 

Array of ACPS jet torque vectors 

K 

K 


Sign of rate change 

RALPH A 

K a ’ 

| 


KOMEGA 

K “ ! 

> 

Filter gains 

KTHETA 

k 9 ! 

i 


KILL1 

KILL1 

ACPS jet continuity indicator 

[ MAGACC] 

[a c ] 


Array of magnitudes of control accelerations 

MAG-U„D 

I u dI 


Magnitude of disturbance acceleration 

MAGOJ-C 

|U C I 


Magnitude of control acceleration 

MP2U 

-P2U 

Torque in the +u direction by 2 u-jets 

MM2U 

-M2U 

Torque in the -u direction by 2 u-jets 

MP2V 

^ P2V 

Torque in the +v direction by 2 v-jets 

MM2V 

Mm2V 

Torque in the -v direction by 2 v-jets 
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Table 6-5: 

Glossary for On-Orbit Phase (cont) 

MNEMONIC 

ENGINEERING 

TERM 

DEFINITION 

MP2W 

M 

— P2W 

Torque in the +w direction by 2 w-jets 

MM2W 

— M2W 

Torque in the -w direction by 2 w-jets 

MCP 

^CP 

Control to pilot transformation matrix 

MFC 

M PC 

Pilot to control transformation matrix 

MIN IMP 

t IMP 

Time of a minimum impulse ACPS jet firing 

NETACC 

[a] 

Net acceleration 

OMEGA BA R 

w' 

Extrapolated angular velocity 

OMEGAHAT 

to 

Estimated angular velocity 

-PLL 

-PLL 

Left angular position error limit 

+PLR 

+PLR 

Right angular position error limit 

P-CEPT 

P 

cept 

Intermediate position error axis intercept 

PHASE 

PHASE 

Phase of the mission 

±RL 

±RL 

Upper and lower rate error limits 

POST.TION_EH ROH [9 ] 

e 

Angular position error 

RATE_ERROR 


Angular rate error 


REALSEC 

REALSEC 

Scale factor for scaling computer time [ = - ] 

\ TICK_TO_SEC / 

RUNTIME 

RUNTIME 

Runtime of the computer in computer ticks 

[ SGNC ] 

[SGNC] 

Signs of the angular rate components 

TARGET_ARL 

TARGET_ARL 

Lower rate limit average 

TF 

fc fast 

Fast sampling period 

TM 

^medium 

Medium sampling period 

TS 

^slow 

Slow sampling period 

TEMPI 

TEMPI 

Temporary vector of body angle errors in control coordinates 

TEMP2 

TEMP2 

Temporary vector of body rate errors in control coordinates 

TEMP3 

TEMP3 

Temporary vector of acceleration in control coordinates 

[T.JET] 

4 JET 

Time of an ACPS jet firing in seconds 
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Table 6-5: 

Glossary for On-Orbit Phase (cont) 

MNEMONIC 

ENGINEERING 

DEFINITION 


TERM 



THETABAR 

Ik 


Extrapolated body angles 

THETAIIAT 

Ib 


Estimated body angles 

TNOW 

t 

n 


Time of present gimbai angle measurement 

[ TJETQ] 

tjETQ 


Time of an ACPS jet firing scaled to the computer 

({A TE-INT 

-B 


"Measured" body angles 

[TLAST] 

X 

te u ] 


Time of the last gimbai angle measurement 

[TIIETA_D] 


Desired set of gimbai angles 

TICKJTOJSEC 

TICK_TO. 

.SEC 

Scale factor for scaling computer time 

[U_D] 

[u D ] 


Array of disturbance accelerations 

W-D 

«d 


Desired angular velocity 

X-BAR 

X 


Average angular position error 

X_DB 

X DB 


Deadband of limit cycle 

X_D 

X D 


One component of Aw the desired change in angular velocity 

NO-CHANGE 

X D old 


One component of the last desired change in angular velocity 

XYZERROR 

-Be 


Errors in x, y, z axis rotation angles 

NOTATION CONVENTIONS 



© 

Modular 

addition 

o 

Modular 

subtraction 


Measured quantity 

* 

Estimated quantity 

/ 

Extrapolated quantity 



Vector 




Matrix 



[ ] 

Array 
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SECTION 7 


BOOST PHASE 


7. 1 Introduction 

The baseline launch configuration control system software described in this 
section satisfies the requirement to control the booster from lift-off through main 
engine cutoff at external tank separation. The system is adaptable to various control 
laws, including attitude control, drift minimum control (with and without accelerome- 
ter feedback) and/or load relief control through proper choice of gain schedules. It 
also contains the option of using rudder and thrust vector control combined in the roll 
channel. 

The following paragraphs contain: A description of the control system with 
functional block diagrams of the pitch, yaw, and roll channels; an annotated flow dia- 
gram; and a discussion of the control gains. 

7. 2 Functional Description 

The baseline launch control system is a digital autopilot receiving body attitude 
from an inertial reference system, body rates from three sets of body mounted rate 
gyros, and translational accelerations from two sets of body mounted accelerometers. 
The current baseline includes rate gyros at the aft end of each solid rocket booster 
and the orbiter. Normal and lateral accelerometers are located in the forward avionics 
bay and the aft end of the orbiter. The autopilot also utilizes prestored engine deflec- 
tion commands and acceleration commands. The former are the nominal trim values 

for the expected center-of-gravity and aerodynamic histories, and the latter are the 
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normal linear accelerations in the nominal trajectory. The final inputs to the 

control loop are the attitude commands, which may either be prestored or come from 
a guidance system. 

From the inputs, the control system determines and sends out deflection com- 
mands for the engines and the rudder. A general optimal mixing logic has been de- 
signed and may be installed at a later date. It has been designed to distribute the de- 
sired torques efficiently among the control actuators while accounting for failures and 
saturation. ^ ^ 

Table 7-1 summarizes the control system inputs and outputs. 


7-1 


Table 7-1 INPUT/OUTPUT SUMMARY 


INPUTS 

• Inertial attitude angles 

• Body rotational rates 

• Y & Z translational accelerations 

• Prestored engine deflection commands 

• Prestored accelerations commands 

• Attitude commands from guidance system 


OUTPUTS 

• Deflection commands to orbiter engines 

• Deflection commands to SUM engines 

• Deflection commands to aerosurfaces (rudder) 

Variable control gains have been calculated to satisfy the rigid body response 
requirements. They are stored as a function of time from lift-off; the numerical 
values are presented in Subsection 7. 4. In the future, the design will be extended to 
include consideration of slosh and bending, and these prestored gains will be replaced 
by filters. All filters will be implemented as second order over second order. Fil- 
ter inclusion in two places in each sensor loop recognizes common filtering require- 
ments and reduces the order of the filter. 

Some aspects of the design presented in this section are specific to the first 
stage of boost (prior to SRB separation). However, control of the orbiter during the 
second stage will utilize basically the same control law. The major differences are 
that the rudder will not be used and there will be no acceleration feedback in the sec- 
ond phase. Furthermore, there will only be one set of body mounted rate gyros 
available. These gyros will be located in the aft section of the orbiter. All filter 
coefficients and gains will be changed at separation and vary as the mass properties, 
flexible body and fuel slosh characteristics vary. 

7. 2. 1 Roll Channel Description 

Referring to Fig. 7-1, the block in the lower right hand corner labeled ’Vehicle 
and Subsystem Dynamics” represents actuator and sensor dynamics, possible analog 
prefiltering of sensor outputs, input output voting, redundancy management, A/D and 
D/A converters. Inputs to the control system from ’’Vehicle and Subsystem Dynamics” 
are the inertial platform attitude angle and rate measurements from three sets of 
body mounted rate gyros. Outputs from the control system to the vehicle are engine 
deflection commands. 
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Fig. 7-1 Baseline Launch Configuration Control System Software Roll Channel 











Following the attitude loop in Fig. 7-1, the platform gimbal angle, 0, is differ- 
enced with the inertial attitude command, 0 c , to produce the inertial attitude error', 

0 E r Ttlis inertial attitude error is transformed to a body roll error, by adding 

the inner gimbal error times the sine of the middle gimbal angle (0„_, = 0 + 0 s i n <-). 

EB El El 

It should be noted that a limiter has been included in the attitude error loop. 

This limiter reserves some TVC authority for rate damping (even though the attitude 
error signal may be saturated) thus preventing oscillation and complete loss of con- 
trol because of large signal instability. The limiter is set to a value less than 
8-degrees. 

Returning to the "Vehicle and Subsystem Dynamics 1 ’ block in Fig. 7-1, three 
sets of roll rate gyro signals are filtered then blended and again scaled and filtered 
prior to being summed with the limited attitude error signals. The resulting signal 
is sent to the rudder and engine mixing logic for distribution to the individual engine 
and rudder actuators. Then each actuator command is limited to eight degrees to 
prevent the roll channel from usurping all pitch and yaw controllability. The rudder 
command is limited to 15 degrees, 

7. 2. 2 Yaw Channel Description 

The yaw channel block diagram in Fig. 7-2 is laid out in the same manner as 
the roll channel with the "Vehicle and Subsystem Dynamics" block in the lower right 
hand corner, the "inertial Attitude Command" in the upper left hand corner and the 
"Actuator Mixing Logic" on the right. The rudder is not used for yaw control, be- 
cause there is adequate engine control authority to handle yaw-plane torque. 

Three features which were not included in the roll channel appear in the yaw 
channel. Two accelerometer inputs have been added to the attitude loop, a forward 
loop integrator follows the limiter, and a rate command signal is differenced with 
the input from the rate gyros. 

The platform gimbal angle, <1, is differenced with the inertial attitude com- 
mand, <i c , transformed to body attitude error (<^ EB = <1^ cos 0 - 0 sin 0 cos d ) 
and scaled. The attitude error signal is then modified by summing with an accelera- 
tion error signal and limiting. The acceleration error is obtained by blending appro- 
priately filtered signals from two body mounted accelerometers and then applying gain 
and filtering compensation according to the control law or bending requirements. 

The accelerometer gain is varied with time to meet the drift minimum or load relief 
requirements. The forward loop integrator trims out steady state errors resulting 
from biases and thrust misalignments. 

The rate commapd signal is introduced because the normal attitude command 
changes at a predictable rate. The rate command signal minimizes the tendency for 
vehicle attitude to lag behind the attitude command due to rate feedback. If the 
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vehicle rolls appreciably the rate command must be applied to the yaw channel. In 
operation, the attitude command is differentiated, multiplied by the sine of 0, and 
differenced with the rate feedback signal. The rate error signal is then added to the 
attitude error signal and sent to the actuator mixing logic routine. Actuator command 
limiting is performed for the same reason stated for roll control. 

7.2.3 Pitch Channel Description 

Comparison of Fig. 7-3, which is the pitch channel block diagram, with Fig. 7-2 
(yaw channel) reveals that the only difference is the addition of prestored engine de- 
flection and acceleration commands. Both of these command schedules are the nom- 
inal values from a reference trajectory. 

The addition of the prestored acceleration command permits use of accelerom- 
eter feedback with any feasible trajectory to modify the attitude error signal to con- * 
form to the desired control law. 

The prestored engine deflection command schedule provides coarse control of 
the engine gimbal trim in the pitch plane. Center of gravity travel alone could be 
handled by the integral compensation, but aerodynamic pitching moment increases 
sharply in the vicinity of Mach one and decreases again thereafter. A high integrator 
gain would be required to follow these rapidly changing deflection requirements. 

This is undesirable because it degrades the rigid body stability margin; therefore, 
the prestored commands are used to control the nominal part of this variation. 

The rest of the pitch channel is the same as the yaw channel. All engine com- 
mands are sent to the actuator mixing logic routine. 

7. 3 Logic Flow Description 

A flow chart showing the logic flow of the control system is presented in Fig. 7-4. 
7. 4 Values of Control Gains 

The basic control system is independent of configuration, assuming that thrust 
vector control is primary for pitch and yaw, and that the principal guidance commands 
will be in the pitch plane. It includes no provision for use of elevons during boost. 

The control gain and command schedules are, of course, configuration depend- 
ent . They will be affected by modifications to either the vehicle configuration or 
the reference trajectory. This will not impact the design, however, because all con- 
trol gains and commands are continuously variable input functions of time. Some of 
them may subsequently be modified t‘o dunctions of other trajectory variables, such 
as altitude and velocity, but this will involve only minor software modifications. 

Table 7-2 presents the current numerical values of the control gains. 
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SECTION 8 


Phase 

Mnemonic 

DFCS 

Symbol in Text 

GLOSSARY 

Definition 

4 

A 1,2, 3,4 

a l* a 2’ a 3’ a 4 

Gains for aileron control 

4 

A 11,12,21,22 

a ll’ a 12’ a 21’ a 22 

Lateral fixed gains 

8 

ABC ( 1 ) 

4 roll 

Total roll channel command 

8 

ABC (2) 

6 PITCH 

Total pitch channel command 

8 

ABC (3) 

5 YAW 

Total yaw channel command 

8 

A E F ( 1 > 

0 EB 

Body roll attitude error 

8 

AB E ( 2) 

®EB 

Body pitch attitude error 

8 

ABE (3 ) 

^EB 

Body yaw attitude error 

8 

A B E F ( 1) 


Scaled roll error 

8 

ABEF (2) 

9 F 

Scaled pitch error 

3 

A B FF ( 3 ) 

4>F 

Scaled yaw ^rror 

2 r * 

ACCBCD 

~b 

a 

Linear accel measurement (body coord) 

4 

ACCG 

g 

Accelerat ion of gravity 

4 

ACCINT (1) 

Ap 1 ^ 

Roll integrated angular acceleration 

4 

ACCENT (2) 


Pitch integrated angular acceleration 

4 

ACCINT (31 

Ar^ 

Yaw integrated angular acceleration 

4 

ACCL 

a L 

Measured lift acceleration 

4 

ACCI.HN 

a LM 

Mean lift acceleration 

4 

A DOT 

a 

Derivative of angle-of- attack wrt time 

4 

ADO (1 ) 

bCL 

Angle- of-at*ack rate filter output 

4 

.ADO (2) 

Aa' 

Angle-of -at tack rate filter extrapolated value 

2,* 

AO 

a G 

Angle-of-a t tack guidance command 

4 

A GOLD 

old 

a G 

Previous-pass value of AG 

4 

AT (1) 

aS 

Angle-of-a t tack filter input 

2 

AI NT 

a INT 

Aileron integrator output 
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DFCS GLOSSARY (Cont) 


Phase 

Mnemonic 

Symbol in Text 

Definition 

2 

AI NT T 


Aileron integrator input 

2 

AINT1 


Aileron integrator noie 

2 

A I N'T f X 

max 

INT 

Aileron integrator output li mi ts ( t-ma x) 

8 

A I V (1 ) 

0 

Platform outer gimbal angle 

8 

AIV (2) 

6 

Platform inner gimbal angle 

8 

AT V (?) 


Platform middle gimbal angle 

8 

AT VC ( 1) 


Platform outer gimbal angle command 

8 

AI VC (2 ) 

e c 

Platform inner gimbal angle command 

8 

A T VC ( 3) 


Platform middle gimbal angle command 

8 

AI VE ( 1 ) 

*EI 

Platform outer gimbal angle error 

8 

AI VE ( 2) 

0 EI 

Platform inner gimbal angle error 

8 

A I V E ( 3 ) 

*EI 

Platform middle gimbal angle error 

4 

ALIM 

a lim 

Limited value of angle-of-at tack filter outout 

-> 

ah 

s 

Measured an gle-of-a t tack 

4 

AMEAN 

OCp 

Mean (trim) angle of attack 

2 

AM 

an 

Normal acceleration (-zbody, with gravity) 

2 

ANF 

an F 

Normal acceleration filter output 

2 

ANGRV 

an GRV 

Gravity corrections for AN 

2 

ANN (I) 


Normal acceleration filter nodes 

4 

AO(1) 

txOC 

Angle-of-at tack filter output 

4 

AO (2) 

Lot' 

Angle-of-at tack filter extrapolated value 

0 

APP(I) 

APP. 

L 

Gains for scaling pitch input to the pitch 




actuators 

8 

APP (I) 

ARP. 

i 

Gains for scaling roll input to the pitch 




actuators 

8 

ARR 

ARR 

Gain for scaling the roll input to the rulder 

0 

ARY (I) 

ary. 

L 

Gains for scaling roll input to the yaw 


actuators 
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DFCS GLOSSARY (Cont) 


Phase 

Mnemonic 

Symbol in Text 

Definition 

4 

AX, AY , AZ 

u ,u ,u 
x* y' z 

Jet acceleration in roll, pitch, yaw 

2 

AY* 

ay 

Measured lateral acceleration 

2 

AY MF 

ay F 

Lateral acceleration filter output 

2 

A Y * N (I) 


Lateral acceleration filter nodes 

8 

A Y Y (I) 

AYY. 

Gains for scaling yaw input to the yaw 
actuators 

2 

AZ M 

az 

Measured normal acceleration («-zbody! 

4 

B 

b 

Fixed estimator parameter 

4 

B 11,12,21 

,22 b 1 1# b 12 , b 21 , b22 

Lateral fixed gains 

4 

BANK 


Roll (bank) about air-relative velocity vector 

4 

3ANKDO (1) 

4 V 

Bank rate filter output 

4 

3ANKD0 (2) 

A 

Bank rate filter extrapolated value 

4 

BANK G 

^VG 

Bank guidance command 

4 

BANK! (1) 


Bank filter input 

4 

BANK* 

*VM 

Mean bank angle 

4 

BANKO (1) 


Rank filter output 

4 

RANKO (2) 


Bank filter extrapolated value 

4 

BDOT 

8 

Derivative of sideslip wrt time 

4 

BDO (1 ) 


Sideslip rate filter output 

4 

BDO (2) 

tb 

Sileslip rate filter extrapolated value 

4 

BD7M 


Rate limit in roll let logic 

2 

BE 


Sideslip error filter input 

2 

REF 

^EF 

Sideslip error filter output 

2 

BEN (I) 


Sideslip error filter nodes 

2 

BG 

% 

Sideslip guidance command 

4 

RI (1) 

A/9 

Sideslip angle filter input 
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Phase 

Mnemonic 

Symbol in Text 

Definition 


BKHOLD 

a old 


4 

*VG 

Previous- pass value of BANKG 

2 


0 

Measured sideslip 

4 

BO (1> 

*0 

Sideslip a nqle filter output 

4 

30 (2) 

ti' 

Sideslip angle filter extrapolated value 



• old 


4 

3 POOL D 

^ VG 

Previous- pass value of bank rate command 

4 

CAT 

COSQLp 

Cos (AHFAN) 

2 

cc 

C c 

CAS command (QC, AC, or CSTAFC) 

2 

CCAHN 


Angle-of-at t ack values for CCttN 

2 

CC A H X 


Angle-of-at tack values for CCEX 

2 

CCCMN 


CAS values (C* or q) for CC*N 

2 

cccnx 


CAS values (C* or q) for CCMX 

2 

CCMN 

min 

c C 

CAS command limit, min (working) 

2 

CC«X 

max 

c C 

CAS command limit, max (working) 

2 

CE 

C E 

CAS error filter input 

2 

CEF 

C EF 

CAS error filter output 

2 

CEN (I) 


CAS error filter nodes 

2 

CF 

C F 

CAS feedback (OF, AF, or CSTAFF) 

2 

CO 

C G 

CAS guidance command (QG or CSTAPG) 

4 

Ct FP 

c ' 

! P 

Scheduled estimator parameter 

4 

CLDAP 


Scheduled estimator parameter 



5a 


4 

CNFP 

C' 

% 

Scheduled estimator parameter 

4 

CNDAP 

c' 

n 

Scheduled estimator parameter 


CR^B 

fiBL 


2,4 

c b 

r 

Ref erence- to- body direction cosine matrix 

4 

CRUD 

C RUD 

Rudder crossfeed gain 

2 

CSTAPG 

C G 

CSTAR guidance command 

4 

DM , DA2 

8 al' ®a2 

Deadbands for pitch jet logic 

- 

DAC 

6a c 

Aileron command 
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DFCS GLOSSARY (Cont) 


Phase Mnemonic Symbol in Text 


Definition 


2 

DACF 

8a CF 

2,“ 

0 A CMX 

. max 

8a c 

2 

DACN (I) 


2 

DAG 

« a G 

8 

DA I VC 

®C 

- 

P A MAX 

§a max 

4 

DAT 

6a t 

2 

DATRTM 


4 

DB 1 , DB2 


2 

DDG 

•o 

- 

DEC 

5e c 

2,u 

DECMN 

.min 

6e c 

2,4 

DEC IX 

. max 
fie c 

2 

DECO 


2 

DECOP 


2 

DEG 

fie G 

- 

DEI.C 

6EL C 

R 

DELL 4,5 

6 4L’ 6 5L 

8 

DELOFF 

6 orb 

8 

DELP(I) 

6 iP 

8 

DEL P ( I ) 

6 iP 

8 

DELR 4,5 

6 4R’ 6 5R 

8 

DELSFM 

6 srm 

8 

DELY (I) 

6 iY 

8 

DELY (X) 

6 iY 

- 

DEMIN, DEMAX 

e min c max 
6 , o 

e e 

- 

DERC 

6ER^, 


Mleron command filler output 
Aileron command limi t ( ♦- max) 

Aileron command filter nodes 
Ai leron guidance command 
Platform inner gimbal rate command 
Physical aileron limit ( + -max) 

Aileron trim estimate 

Effective aileron trim, for mole change 

Deadbands for roll let logic 

TAEtt guidance interface variable 

Elevator command 

Elevator command limit, min 

Elevator command limit, max 

Elevator command from speed loop crossfeels 

BCE version of DECH 

Elevator guidance command 

Left elevon command 

Commands to S P "I left gimbals 

Prestored orbiter engine deflection commanl 
Commands to orbiter pitch gimbals (i = 1,2,D 
Commands to SRM pseudo pitch gimbals (i = u,5) 
Commands to SRM right qimbals 
Prestored SRM engine deflection command 
Commands to orbiter yaw gimbals (i=1 f 2,3) 
Commands to SRM pseudo yaw gimbals (i^a,^) 
Phy sica 1 elevator lower (upper) limit 
Right elevon command 
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DFCS GLOSSARY (Cont) 


Phase 

Mnemonic 

Symbol in Text 

Definition 

4 

DET 

6e t 

Elevator trim estimate 

2 

DSTRIK 


Effective elevator trim, tor mode change 

2 

DFC 

6f c 

Body-flap command 

2 

DFCF 

5f CF 

Flap command filter output 

2 

OFCfl 

5f CM 

Flap command from Mach profile 

2 

DFCMCH (I , J ) 


2 x N matrix of (DFCP1 , flACH) pairs 

2 

DFCMN 

6f ™ ln 

Flap command limit, min 

2 

DFCMNB 


BCS version of DFCNN 

2 

DFCMX 

fif m ax 

Plap commanl limit, max 

- 

DFMIN, nm ax 

,.,miri . .max 
oi , 6 f 

Physical flap lover (upper) limit 

4 

DP1STF 

q|= 

Blending parameter for elevator use 

4 

DP2ST R 

^2 

Blending parameter for aileron use 

4 

DP 3STR 

^3 

Blending parameter for rudder use 

4 

PR1,PF2 

®0l* *02 

Deadbands for yaw jet logic 

- 

DRC 

6r C 

Rudder command 

2 

D p CMX 

«r™ ax 

Rudder command limit (f-max) 

2 

DP F 

6r F 

Rudder loop feedback 

2 

DRG 

® r G 

Rudder guidance command 

- 

D p M A X 

„ max 
6r 

Rudder physical limit (♦-max) 

8 

DP C LI 


Platform outer gimbal rate command 

2 

DR m RIK 


Rudder effective trim, for mode change 

- 

DSC 

® s c 

Speed brake command 

2 

DSCF 

6s cf 

Speed brake command filter output 

2 

DSCHK 

, mi n 
6 S c 

Speed brake command limit, min 

2 

DSC9IKB 


BCS version of DSCMV 

2 

DSCrtX 

. max 
® s C ' 

Speed brake command limit, max 

2 

DSCP 


Past value of DSC 
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Phase 

Mnemonic 

Symbol in Text 

Definition 

2 

DS DM X 

• max 
6s 

Speed brake (command) rate limit 

2 

DSG 

« s G 

Speed brake guidance command 

- 

DSF1IN, DS1 AX 

6s min , 6s maX 

Speedbrake physical upper (lower) limit 

2 

DTC 

5t c 

Throttle command 

2 

DTCF 

8t CF 

Throttle command filter output 

2 

DTC1N 

cl min 

6t c 

Throttle command limit, min 

2 

DTC'lf! B 


BCS version of DTC M N 

2 

DTCMX 

*tr x 

Throttle command limit, max 

2 

DTCP 


Past value of DTC 

2 

DTD*X 

8t maX 

Throttle (command) rate limit 

2 

DTD*1XB 


PCS version of DTDPX 

2 

P T G 

St G 

Throttle guidance command 

- 

DT^ir , DT* A X 


Physical throttle lower (upper) limit 

2,4 

DYNP 

q 

Dynamic pressure estimate 

4 

DYNP 1,2,3 

^1* ^ 2' ^ 3 

Dynamic pressure switch levels for elevator. 




aileron, roll let mode III 

2 

BINT 

e INT 

Elevator integrator output 

2 

EINTI 


Elevator integrator input 

2 

EINTMN 

min 

e INT 

Elevator integrator output limit, min 

2 

EINTf X 

max 

e INT 

Elevator integrator output limit, max 

2 

EINTI 


Elevator integrator node 

4 

KL 1,2, 3, 4 

^1* ^2* ^3* ^4 

Temporary blending parameters 

2 

ET 

e 

Elevator trim tracker output 

4 

ETA 

r) 

lateral scheduled gain (late entry) 

2 

ET DP 

e TDB 

Elevator trim tracker input deadband 

2 

ETDMX 

• max 
e T 

Elevator trim tracker rate limit 

2 

ET F 

e TF 

Elevator trim tracker filter output 
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Phase 

Mnemonic 

Symbol in Text 

Definition 

2 

ETMN 

min 

e T 

Elevator trim tracker input limit, min 

2 

ETMX 

e™ ax 

Elevator trim tracker input limit, max 

a 

FI , F2 

f l* f 2 

Longitudinal blending parameters 

u 

F3, F4 

f f 
3 J 4 

Lateral blending parameters 

4 

FILTFF 

FILTER 

State estimation flag: o=no filter (use measured 

values of state) 

1=filter 

4 

FKLAT 

K 

lat 

Scheduled estimator gains: lateral charnel 

4 

FK LON 

K lon 

Scheduled estimator gains: longitudinal chanrel 

4 

G 1 ,G2 

G l* G 2 

Lateral fixed gains 

2 

GAIN? 


CAS aileron integral sain 

2 

G A I P 


CAS aileron integrator gain ( C A P *G A T NT ) 

4 

GATR 

g a 

Filter qain for air data (=TF/TS) 

4 

G A M 1 1 , 12, 

21,2 y iv r 12 » >21* y 22 

Lateral scheduled gains (late °ntry) 

4 

GAMA 

y 

Flight path angle 

2 

GAN (T) 


Normal acceleration filter coefficients 

2 

GAP 


Roll rate to aileron CAS loop gain 

2 

GAP* 


Roll rate to aileron SAS gain (working) 

2 

GAPMB 


BCS version of G A PM 

2 

GAPffS (I) 


Gain-switch-sensitive SAS values for GAP* 1 

2 

G AX 


SAS roll stick gearing gain 

2 

GAXB 


BCS version of G A X 

2 

GAYK (J) 


Lateral acceleration filter coefficients 

2 

GBE(I) 


Sideslip error filter coefficients 

2 

GCF (I) 


CAS error filter coefficients (working) 

2 

GCEA (I) 


A n gle-of-a t.tack values for GCF (I) 

2 

GCFC (I) 


CAS values (C* or q) for GCE(I) 
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Mnemonic 
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Definition 

2 

GCFQF 


C S T A F q feedback gain 

2 

GCH 


Altitude loop gain ( H to C) 

2 

GCTH 


Theta loop gain (TH to C) 

2 

gcx 


CAS pitch stick gearing gain 

2 

GDAC(I) 


Aileron cmnd filter coef* (A to R crossf? 

2 

GDFC 


Flap command filter coefficient 

2 

GPSC 


Speed brake to elevator crossfeed filter 

2 

GPTC 


Throttle to elevator crossfeed coefficient 

2 

GEC 


CAS loop gain 

2 

GSCA 


Angle- of-attack value for GEC 

2 

GECC 


CAS value (::* or g) for GEC 

2 

GE F 


Flap to elevator crossfeed gain 

2 

GFF3 


BCS version of GEF 

2 

GEIC 


CAS elevator integration qain (GEC*GEINT) 

2 

GETCA 


Angle-of-at t ack value for GFIC 

2 

GFTCC 


CAS value <C* or q) for GEIC 

2 

GEINT 


CAS elevator integral gain (working) 

2 

GEINT A 


Angle-of-at tack value for GEINT 

2 

GEINTC 


CAS value (C* or g) for GEINT 

2 

GEQF 


Elevator damper loop gain (working) 

2 

GEOEA 


Angle-of-attack value for GEQF 

2 

GEQFP 


BCS version of GEOF 

2 

GEOFC 


CAS value (C* or q) for G E OF 

2 

GEQMS (I) 


Gain-switch-sensitive SAS values for GEO F 

2 

GES 


Speed brake to elevator crossfeed gain 

2 

GESB 


BCS version of GES 

2 

GET 


Throttle to elevator crossfeed gain 


d> 


oef . 
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2 

GETB 


BCS version of GET 

2 

GETF 


Elevator trim tracker filter coefficient 

2 

GEX 


SAS pitch stick gearing gain 

2 

OFXB 


PCS version of GFX 

2 

GFX 


Flap controller gearinq gain 

2 

GFXB 


BCS version of GFX 

2 

GHS(I) 


Altitude error filter coefficients 

4 

GNTJi , GNU2 


Lateral scheduled gains (late entry) 

2 

GPF (I) 


Roll rate error filter coefficients 

2 

r.pps (i) 


Bank angle error filter coefficients 

2 

G PH PS 


Psi-loop gain (psi to phi) 

2 

GPHT (I) 


Phi-turn filter coefficients 

2 

GP (T ) 


Roll rate feedback filler coefficients 

2 

GPflB (I) 


BCS version of GPM(I) 

2 

GPPH 


Phi-loop gain (Dhi to p) 

2 

GPSE (I) 


Psi error filter coefficients 

2 

GPX 


CAS roll stick gearing gain 

2 

GO* (I) 


Pitch rate feedback filter coefficients 

2 

GQ P*B ( T ) 


BC c version of GQW (I) 

2 

GQW 


Pitch rate washout filter coefficient 

2 

GRAY 


Lateral acceleration feedback gain 

2 

GPP 


Sideslip feedback gain 

2 

GPCPH 


Phi to RC gain r turn coordination 

2 

GRDAC 


Aileron to rudder crossfeed gain 

2 

GFDACB 


BCS version of GRDAC 

2 

GPF (I) 


Yaw rate error filter coefficients 

2 

GPEB (I) 


BCS version of GRF(I) 
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2 

GPM 


Yaw rate filter coef. (elevator loops) 

2 

GP PE 


Yaw rate feedback gain 

2 

GP PW 


Yaw rate washout feedback gain (working) 

2 

GRPWA 


AUTO version of GREW 

2 

GRPWB 


BCS version of GFPW 

2 

GPPWS (I) 


Gain-switch-sensi t ive SAS values for GRRW 

2 

GRW 


Yaw rate washout filter coefficient 

2 

GP WB 


BCS version of GPW 

2 

GPX 


Pedal gearing gain (working) 

2 

GFXAHT 


RAHT mode value for GRX 

2 

GPXB 


BCS version of GFX 

2 

GPX GAS 


RS AS mode value for GRX 

2 

GSPDG 


TAE^ quidance DDG to speed brake gain 

2 

GST! 


Speed brake velocity loop gain 

2 

G5 HE (I) 


Speed brake values for G n E (I) 

2 

GSX 


Speed brake manual controller gearing 

2 

GSXB 


BCS version of GSX 

2 

GTPDG 


TAEM guidance DDG to throttle loop gain 

2 

GTHE (I) 


Theta error filter coefficients 

2 

GTH 


Throttle velocity loop gain 

2 

G T UE ( I) 


Throttle values for GHE(I) 

2 

GTX 


Throttle manual controller gearing 

2 

GTXB 


BCS version of GTX 

2 

GIJE (I) 


Velocity error filter coefficents (working) 

2 

GX A 


Roll stick filter coefficient 

2 

G X AX 


Roll stick parabolic shaping gain 

2 

GXF 


Pitch stick filter coefficient 
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2 

GXEX 


Pitch stick parabolic shaping gain 

2 

GXF 


Pedal filter coefficient 

2 

GXRX 


Pelal parabolic shaping gain 

2 

GXS X 


Speedbrake controller parabolic shaping gain 

2 

GXSXP 


BCS version of GXSX 

2 

GXTX 


Throttle controller parabolic shaping gain 

2 

GXTXB 


PCS version of GXTX 

2 

HE 

h e 

Altitude error 

2 

HEF 

h ef 

Altitude error filter output 

2 

MEN (I) 


Altitude ^rror filter nodes 

2 

HG 

h g 

Altitude guidance command 

2 

H M 

H 

Measured altitude 

2 

HP EF 

h ref 

EAL' 7 ' mode reference altitude command 

2 

I ABES 


ABES present (1) ; ABES absent (0) 

2 

TAGAIN 


AGAIN switch position (1 thru 4) 

2 

IAIN? 


Aileron integrator active (1); inactive p) 

2 

IANGPV 


Gravity corrections: UNITY(I), C0S T H(2}, 

COSTH * COSPH(l) 

2 

IATNCK 


Pilot requested values for IAINT 

2 

XA^R I M 


Roll stick trim activity: 0, ♦ 1 (rgt ) , - 1 (lef t ) 

2 

IBFDRK 


Beta feedback (1); lat. accel. feedback (n> 

2 

reel 


CAS command switch 11 options (0=off) : 


1 = ECAS mode or control stick steerinq 

2 = E ATT or SALT mode 

3 = thg or HG mole 

a = CG mode (QG or CSTARG) 

5 = AG mode 
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Phase Mnemonic 
2 ICC2 

2 ICF 

2 IC ST/. P 

2 IEGATN 

2 TETNCK 

2 IEINT 

2 IETRItf 

2 IGFAPX 

2 IHC 

U IIAT 

1*1 

TNGUIP 


Symbol in Text 


t lat 


IM 


INGUID 


INSENS INSENS 


2 IPC 


Definition 

CAS command switch *2 options (0 = off) : 

1 = theta loop (TH3EF or THG) 

2 = altitude loop (HPEF or HG) 

CAS feedback: 0(1), CSTAR(2), ALPHA ( S) 

CAS mode options: CSTAR(I), FCAH(O) 

EGAIN switch position (1 thru 4) 

Pilot- requested values for IEINT 

Elevator inteqrator active (1); inactive (0) 
Pitch stick trim activity: 0 , «■ 1 (f vd) , - 1 (a f t) 

Extend gear (1) f retract gear (0) 

Altitude command: HREF (1), HG (2) 

Mode flag for entry phase lateral control 

0 = early entry, "reverse" aileron 

1 * late entry, "proverse" aileron 
Hedium-sampling-ra te index: 

0 = perform calculations this pass 
Flag specifies guidance command set 

neg=command s not updated this pass 
+ /- 1 A NKG , ALPHG specified 

♦ /- 2=30110, QG specified 
Flag specifies source of measurements 

0=provided by simulator (hybrid sim.) 
1=sensors or environment (SSFS) 

2=Gf>N and sensors or environment 
(SSFS) 

Roll rate command options (0=off) : 

1 = ACAS mode or control stick steering 
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IPHC 


iphu 

IPHPC 


2 IPSC 

2 IQTRN 


2 

2 

2 

2 

2 

2 


TRCHAI 

IRCHFI 

TRGAIN 

IRTRTK 

TRW 

TR WP 

rs 


IS 


ISTAFT 


ISTART 


Def init ion 

2 = AA-T or A KDG mode 

3 = PHG or PSG node 
^ = PG mode 

Bank angle command options (0=off ) ; 

1 = A A^T mode 

2 = A HI)G mode 

3 = PSG mode 
*4 = PHG mode 

PP*1 valid p), invalid (0) 

Phi sourcfi for yaw rate command options: 

1 = zero (beta /r-wash . fdhk) 

2 * P H 'l (measured) 

3 = PHC (commanded, if available) 

Psi command: PSFEF(I), PSG (2) 

Pitch rate turn coordination option: 

1 = off, use 0 feedback 

2 = use q washout feedback 

3 = use Q=5*TANPHI feedback 

Request change IAIN? state (1), no request (0 ) 
Request change IEINT state (1), no request (0) 
RGAIN switch position (1 thru 4 ) 

Pedal trim activity: 0, ♦ 1 (lft ) , - 1 ( r gt) 

Taw rate washout act ive ( 1 ) , i nac ti ve ( 0 ) 

Pal load indicator for IRW 
Slow-sampling-rate index: 

0 = perform calculations this pass 
Initialization pass(1), not in itializati on (0) 
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_ 

TSTRTP 

ISTRTR 

Request initialization ( 1) r no request (0) 

2 

IT HC 


Theta command: THFEFp) , THO(2) 

- 

ITM, ITS 

ITM, ITS 

Number of fast passes per medium (slow) pass 

4 

IXX 

r xx 

Fixed estimator parameter 

4 

TZZ 

*zz 

Fixed estimator parameter 

2 

.TACAIN 


local past-value copy of IAGAIN 

2 

JFGAIN 


Local past-value copy of IFOAIN 

2 

ip; e af x 


Local copy of IGF A R X 

2 

JFCHAI 


Local copy of IRCHAT 

2 

JRCHFI 


Local copy of IRCHEI 

2 

JPGAIV 


Local past-value copy of IRGATN 

4 

KA, KE 

K , K 
a J e 

Aileron (elevator) deflection limit parameter 

4 

K A I.F 1 

K i 

al 

Elevator trim q a i n (early entry) 

4 

K A LF2 

K 0 
0.2 

Elevator trim gain (late entry) 

U 

KC1,KC2 

KC l, KC2 

Low-dynamic- pressure constant filter qairs 

4 

KPCU 

K * 

Aileron trim gain 

8 

KX (I) 

K Yi 

Individual y accelerometer gain 

8 

KZ (I) 

k zl 

Individual z accelerometer gain 

4 

LA M DAA,LAHDAE X a * 

Switch set to 1 if pitch (lateral) A CPS jets 
are needed 

4 

MA, MF 

m , m 
a # e 

Slope of KA(KE) curve 

2/4 

MACH 


Mach number 

2 

MOD A 


Current aileron loop mode: 

1 = ARCS 2 = A S A S 3 = ACAS 
4 = A AMT 5 = A ATT f = A HDG 

2 

MOD ACH 


Indicate aileron loop mode change 
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2 MODAG 

2 MGDAGR 
2 MOD AF 

2 MOTR 

2 MOD3P 

2 MOTE 

2 1 0 D EC H 

2 10 DEG 

2 MODEGP 

2 MODEF 

2 10 DR 

2 10DRCH 

2 10DRG 

2 MOD RG R 


3 - 1 0 D A R chanqe 4 - first pass 
Current aileron guidance mode (M0DA=4): 
1 = PG 2 - PHG 3 = PSG 
Requested values for MODAG 
Requested values for 1CDA 
BCS mode (1), not RCS mode (0) 

Requested values for mqdb 
C urrent elevator loop mode: 

1 = EBCS 2 = ES AS 3 = EGAS 
U = EAUT 5 = EA T T 6 = EALT 
Indicate elevator loop mode change 

1 = no change 2 = MOD EG E chanqe 
3 - MODFR change 4 = first pass 
Current elevator guidance mode (MODK-U): 
1 = QG 2 = THG 3 = AG H = HG 
Requested values for MODEG 
Requested values for MODE 
Current rudder loop mode: 

1 = RBCS 2 = RSAS 
3 = PC AS (=2=FSAS) 4 = RAUT 

Indicate rudder loop mode chanqe: 

1 = no chanqe 2 = MODRGE change 
3 = MODFR change 4 = first pass 
Current rudder guidance mode (10DR=4) : 

1 = yaw rate feedback 

2 = yaw rate washout feedback 
Requested values for 10DRG 
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Definition 


2 10IPP 

2 MODS 

2 MOPT 

2 MO DU 


2 MODTJCH 
2 MODUFP 

2 MOPIIFT 

2 MODUG 

2 MODUGR 

2 MODUR 

4 MU 

4 P 11, 12,21 ,22 

2 PC 


Requested values for MOPP 
Current speed brake mode: 

1 = PCSSP 2 = MANSB 3 = IASSB 
4 = ATITCSB S - DSGSB 6 = TAEKS 
Current throttle mode: 

1 = BCSTH 2 - MANTH 3 = I A STH 

4 = AUXOTH 5 = DTGTH 6 = TAE M TH 
Current, speed control mode: 

1 = UBCS 2 = (JMAN (IAS-man) 

3 = UA3B (IAS-sb) 4 = U A 'IT 

5 = HATH (lAS-th) 

Indicate speed control loop mode change: 

1 = no change 2 = MODUGF change 
3 = MODUR change 4 = first pass 
Flag for FIXUP subroutine sped control 

1 = MANUAL or AUTO-dir 2 - lAS-hold 

3 = A UTO (UG- quid) 4 = AUXO ( T AEM-gil i H 

Flag for velocity filter moding: 

0 = UFILT inactive 1 = SB-ufilt 

2 = TH-ufilt 

Current guidance speed control mode: 

(1-10, see speed control flag table) 
Requested values for MODUG 
Requested values for MODI! 

Roll jet switch parameter 

p ll* p 12* p 21* p 22 Lateral derived qains (late entry) 

Roll rate command 
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8 

PCFCK 

k pa 

Pitch channel error gain 

2 


max 

P C 

Roll rate command limit, ♦-maximum 

2 

PE 

P E 

Roll rate error / error filter input 

2 

PFF 

P EF 

Roll rate error filter output 

8 

PECK 

K 9 

Pitch error gain 

2 

PEN (T) 


Roll rate error filter nodes 

3 

PEP 

P ER 

Limited pitch channel error 

8 

PEP I NT 

P ERINT 

Weighted inteqral of PER 

2 

PC 

P G 

Roll rate guidance command 

- 

PH AS F 

PHASE 

Flight, phase flag (internal copy of PHASER) 

- 

PHASER 

PHASER 

Flight phase flag (external) : 1 ^rollout phase 


2=cruising phase 
4 = e n t r y phase 
B = orbi t TVC phase 
f = orhit. RC5 phase 
7=insert ion TVC phase 





8=booster TVC phase 

2 

PHC 

0 C 

Bank angle command 

2 

PHCrtX 

. max 
0 C 

Bank angle command limit, ♦-maximum 

2 

PHE 

0 E 

Bank angle error / error filter input 

2 

PHEF 

0 EF 

Bank angle error filter output 

2 

P HFN (I) 


Bank angle ^rror filter nodes 

2 

PHC 

0 G 

Bank angle guidance command 

4 

PHILAT 

c 

Estimator state transition matrix: lateral 




channel 

4 

PHILIP 

0 lim 

Limited value of bank filter output 

4 

PHILCN 

.CL 

*ion 

Estimator state transition matrix: longitudinal 


channel 
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2 

PH M 

0 

Measured bank angle 

2 

phpef 

0 REF 

AATT mode reference bank angle command 

2 

PHSEI 

^ SEL 

Pilot- selected values for PH BE" 

2 

PHTF 

0 tf 

Phi-turn filter output and/or 

bank anqle for turn-coord ination (e S 

2 

PH7N { T ) 


Phi-turn filter nodes 

4 

PT (D 

Ap 

Roll rate filter input 

8 

PT. TP' 

PLIM 

Pitch error software limit 

2,* 

pm 

P 

Measured roll rate 

2 

{' TS 

p dis 

Poll rate for * A^T or A.HDG reengage after 
ACSS 

2 

?M p 

P F 

Roll rate filter output 

2 

P M F P 


BCS version of PMF 

T 

A 

P M ND (I) 


BCS version of PMN(I) 

2 

PAN (T) 


Roll rate filter nodes 

4 

PN 

p N 

Nominal roll rate 

8 

PP ATFC 

PRATE C 

Pitch rate command 

8 

P P A T F F 

PRATEp 

Scaled pitch rate error 

8 

pp a teh 

PRATE m 

Blended pitch rate 

8 

PPGG (I) 

k ql 

Individual pitch rate gyro gain 

8 

PRGK 

k q 

Pitch rate gain 

■> 

PSE 

¥ 

E 

Heading error 

2 

PS FF 

^EF 

Heading error filter output 

2 

PSTN (I) 


Heading error filter nodes 

2 

PSG 

*G 

Heading guidance command 

2 

PS M 

* 

Measured heading 

2 

PSFEF 

\P 

kEF 

AHDG mode reference heading command 
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2 

pqc-p] 

'F 

SEL 

P i lot-solectpi values for PSFLE 

-> 

i- 

oc r 

q CT 

Pitch turn coordination g-command 

2 

QCT^X 

q max 

q CT 

OCT maximum (0. - min) 

2 

Q C, 

q G 

Pitch rate guidance command 

4 

oi P) 

Aq 

Pitch rate filter input 

2,4 

0* 

q 

Measured pi^ch rate 

2 

o m r i , c 

q DIS 

Pitch rate f or or HALT renngage 




ECSS 

2 

om p 

q F 

Pi^ch rate filter output 

2 

or m 


BCS version of UTF 

2 

orN (i) 


Pitch rate filter nodes 

2 

0*NB (T) 


PCS version of Q*N(T) 

4 

ON 

q N 

Nominal pitch rate 

2 

QWF 

q WF 

Pitch rate washout filter output 

2 

0 w i 


Pitch rate washout filter node 

4 

P 1 

l 'l 

Angle-of-attack error limits 

4 

R8 

r 8 

Poll angle q rror limits 

4 

? 1 1, 12,21 ,22 

r ll» r 12' r 21’ r 22 

Lateral derived gains (late entry) 

8 

R A ('TV K 

K 0 

Roll error gain 

2 

PC 

r c 

Yaw rate command 

4 



Derivative of bank angle 

4 

PPOTC 

0 C (0 C ) 

Bank rate command 

2 

? 2 

r E 

Yaw rate error 

2 

p p p 

r EF 

Yaw rate error filter output 

n 

PEN (I) 


Yaw rate error filter nodes 

2 

a E N B ( I ) 


BCS version of FEN(2) 

8 

PEP 

r er 

Limited roll error for TVC 
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2 

PG 

r G 

Yaw ra to guidance command 

2,4 

P»0 


Reference air density 

4 

RT (1) 

Ar 

Yaw rate filter input 

8 

FT.IM 

RLIM 

Roll error software limit 

2,4 

PI 

r 

Measured yaw rate 

2 

PMf 

Ff 

Yaw rate filter output (elevator loops) 

4 

P N 

r N 

Nominal yaw rate 

2,4 

POLL 

0 

Computed roll angle 

2 

ROT.LG 

0 G 

Guidance commanded roll anqle 

4 

PPPT 1 , 2, H 

r 0l' r 02’ r 03 

Parameters for maneuver rate calculation 

4 

HP 

— rel 

Position vector in reference coordinates 

8 

RFATFF 

RRATEp 

Scaled roll rate for T VC 

8 

R p ATE r 

rrate m 

Blended roll rate signals 

8 

P?f GO 

K Pi 

Individual roll rate gyro qain 

8 

p RG^ V * 

K p 

T VC roll rate gain 

2 

R WEB 


PCS version of PWF 

2 

R W 1 B 


BCS version of R W 1 

2 

PWF 

r WF 

Yaw rate washout filter output. 

2 

PW 1 


Yaw rate washout filter node 

4 

:S 

S 

Fixed estimator parameter 

4 

SA T 

sina T 

Sin (AMEAN) 

4 

SIGMA A,B,P 

Q 

R. 

Q 

Q 

^6 

Gains for let logic - pitch, roll, yaw 

8 

T A 01 NT 

T 

Integrator jain 

2 

TTGDFX 


TF * speed brake rate limit 

2 

TPTD^X 


TF * throttle rate limit 

2 

TE M P 


Temporary local variable(s) 

2,4 

TF* 1 P F p 

t a 

Air temperature, deq F and dimensionless 

2,4 

TEMPE K 

tk a 

Input air temperature in deg. Kelvin 
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2 

T FTD* X 


TF * elevator trim tracker rate limit 

- 

TF 

t f 

Fast sampling interval 

2 



/ 2. 

2 

TFP 


BCS version of TF2 

? 

T G S P P 0 


* TAFM/S3 loop gain 

2 

TGTDDG 


TF * T A EM / T M loop gain 

2 

T f l E 

®E 

Theta error 

2 

THEF 

e 

EF 

Theta error filter output 

2 

THEN (I) 


Theta error filter nodes 

2 

THG 

e 

G 

Theta guidance command 

o 

i. 

THM 

V 

Measured theta 

2 

TH3EF 

0 REF 

EATT mode reference pitch attitude command 

2 

THSE l 

®SEL 

Pilot- selected values for T H F E p 

- 

TM 

t m 

Medium sampling interval 

- 

T N EXT 

TNEXT 

Next PFCS entry time 

- 

T S 

T 

S 

Slow sampling interval 

2 

TXADMX 


TF * roll stick r a t ^ limit 

2 

TX A TP 


TF * roll stick trim rate 

2 

T X A TF B 


BCS version of T X A T P 

2 

TXFDMX 


TF * pitch stick rate limit 

2 

TXE^F 


TF * pitch stick trim rate 

2 

TXFTFF 


BCS version of TXFtp 

2 

TXFDMX 


TF * pedal rate limit 

2 

TX FTP 


TF * pedal trim rate 

2 

TXPTFB 


BCS version of TX P T R 

2 

TXTDMX 


TF * throttle controller rate limit 

2 

T2G AIP 


TF2 * roll integrator gain 
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2 

T2GE t C 


TF2 * pitch inteqra*or gain 

4 

U1 ,02 

U l' U 2 

Roll let logi c gains 

2 

u f 

U E 

Velocity error 

2 

!JEF 

U EF 

Velocity error filter output 

2 

OEF (7) 


Velocity error filter nodes 

2 

II G 

U G 

Velocity quidance command 

4 

OJ P ) , (2) , 

(3) u x , u y . u z 

ACPS torque commands 

4 

UL (T) 

a-L 

Vector in longitudinal plane normal to velocity 

2 

on 

U 

Measured velocity 

2 

TJF FF 

U REF 

IAS-hold mode reference velocity command 

2 

OS EL 

U SEL 

Pi lot- selected values for tJREf- 

2,4 

UX V,TJYV,M7, V 

uv , uv , uv 

X* — y — z 

Unit vectors along the velocity coordinate 
system, in reference coordinates 

4 

VBPDY 

-body 

Air-relative velocity in body coordinates 

4 

VM ACM 


Lateral switch parameter 

4 

VR 

Yrel 

Air-relative velocity in reference coordinates 

4 

VI ,W2 

Wl» w 2 

Cains for aileron control logic 

8 

wn m.d 

rrate l 

Individual roll rate qyro signal 

8 

WB (2,1) 

PRATE l 

Individual pitch rate gyro signal 

8 

WR (3, I) 

YRATE l 

Individual yaw rate qyro signal 

4 

WD 

w d 

Gain for elevator control 

2 

XA 

xa 

Local variable for XA* processing 

2 

X A R 

xa B 

Roll stick detent bias 

2 

XAD9 

xa DB 

Roll stick deadband 

2 

XAOMX 

• max 
xa 

Roll stick rate limit 

2 

X A F 

xa F 

XA filter output 

2 

X A* 

xa 

Measured roll stick deflection 
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X» rDIS 

Xa DIS 

2 

XAMX 

max 

xa 

2 

X A P 

xa p 

2 

X A T 

xa T 

2 

X A TB 


2 

XATtfX 

max 

xa T 

2 

XATHXP 


2 

X A TR 

xa T 

2 

XATRB 


2 

XE 

xe 

”> 

XE B 

Xe B 

2 

XEPB 

xe DB 

2 

XEP*X 

max 

xe 

2 

XEF 

xe F 

2 

XEM 

xe 

2 

XF'IDTS 

Xe DIS 

2 

XEtfN 

min 

xe 

2 

XEM X 

max 

xe 

2 

XEP 

xe p 

2 

X E T 

xe^ 

2 

XFTB 


2 

XETMN 

min 

xe p 

2 

XETMNP 


2 

X E'^’IX 

max 

xe p 

2 

XETUXP 


2 

X ETR 

xe p 

2 

XFTRB 



Definition 

Roll stick disengage level for A CSS 
Roll stick limits (♦-max) 

Processed roll stick 

Poll stick trim 

BCS version of X AT 

Roll stick ♦rim limit (♦-max) 

PCS version of XA Ttt V 

Roll stick ♦rim rate 

BCS version of XATP 

Local variable for XE'P processing 

Pitch stick detent bias 

Pitch stick deadband 

Pitch stick rate limit 

XE filter output 

Pleasured pitch stick deflection 

Pitch stick disengage level for ECSS 

Pitch stick minimum (aft) 

Pitch stick maximum (fvd) 

Processed pitch stick 
Pitch stick trim 
BCE version of XF^ 

Pitch stick trim minimum (aft) 

BCS version of XETfIN 

Pitch stick trim maximum (fwd) 

BCS version of XE m flX 
Pitch stick trim rate 
BCS version of XE m R 
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DFCS GLOSSARY (Cont) 


Phase 

Mnemonic 

Symbol in Text 

Definition 

2 

X F R 


Flap controller bias 

2 

X F M 

£t 

Measured flap cont roller deflection 

2 

XFMMN 

~ r min 

xf 

Flap controller minimum 

4 

XI 1 


Scheduled control qain: longitudinal 

4 

XT 2,3, 4, 

^ 62* ^3* ^4* ^5* ^6 

Scheduled control gains: lateral, early entry 

2 

XR 

xr 

Local pedal variable 

2 

XPB 

xr b 

Pedal detent bias 

2 

XPDB 

xr DB 

Pedal deadband 

2 

XFDMX 

* max 
xr 

Pedal rate limit 

2 

X° F 

xr F 

XR filter output 

2 

XFM 

xr 

Measured pedal deflection 

2 

XPMX 

max 

xr 

Pedal limits (t-max) 

2 

XRP 

xr p 

Processed pedal 

2 

X? rr 

xr p 

Pedal trim 

2 

X P TB 


BCS version of XPT 

2 

XPTMX 

max 

xr p 

Pedal trim limits (♦-max) 

2 

XRTMXB 


BCS version of XR T MX 

2 

X ? TH 

xr T 

Pedal trim rate 

2 

TFTRP 


BCS version of XRTR 

2 

xs 

xs 

Local variable for XSM processing 

2 

XSM 

xs 

Measured speed brake controller deflection 

2 

XSMDIS 

XS DIS 

XSM disengage level for IAS-sb 

2 

XSMMN 

~ min 
xs 

XSM minimum (retract) 

2 

XSXREF 

XS REF 

XSM reference for T A S - s b disengage 

2 

XSP 

xs p 

Processed XSM 

2 

XT 

xt 

Local variable for XT M processing 

2 

xTr 

xt 

Measured throttle controller deflection 
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DFCS GLOSSARY (Cont) 


Phase 

Mnemonic 

Symbol in Text 

2 

XTflDIS 

xt 

DIS 

2 

XTY1N 

~min 

xt 

2 

XTMPFF 


2 

XTP 

xt p 

2 

X^PB 


2 

XTF P 


2 

XTPP8 


8 

YACCFP 

y f 

8 

YACIF (I) 

Y. 

i 

8 

YACSIG 

*M 

8 

Y A OK 

k ay 

8 

YCFGK 

k ya 

8 

YEGK 

K * 

6 

YEP 

y er 

8 

YEPINT 

y erint 

8 

YLIM 

YLIM 

8 

YP ATFC 

YRATE C 

8 

Y R ATE F 

YRATE_ 

F 

8 

YR ATF W 

YRATE m 

8 

YRGG (I) 

k rl 

8 

YFGK 

Kr 

8 

ZACCEP 

• • 

Z F 

8 

Z A CD 

• • 

Z D 

8 

ZACIN (I) 

Z. 

L 

8 

ZACSIG 

Z M 

8 

ZAGK 

k az 

4 

ZETAI 


4 

ZFTA 1,2 

4 1' ^2 


Definition 

XT« disenqaqe level for IAS-th 
XT* minimum (idle) 

XTW reference for IAS-th disengage 

Processed XTtf 

BCS version of XTP 

Past value of XTP 

BCS version of X^PP 

Scaled y acceleration error 

Individual y accelerometer signal 

Blended y accelerometer signals 

Y acceleration gain 

Yav channel error gain 

Yaw error gain 
Limited yaw channel error 
Weighted integral of YER 
Yaw error software limit 
Yaw rate command 
Scaled yaw rate error 
Blended rat^ gyro signals 
Individual vaw rate gyro gain 
Yaw rate gain 

Scaled z acceleration error 
Z acceleration command 
Individual z accelerometer signal 
Blended z accelerometer signals 
Z acceleration gain 
Gain for elevator control 
Gains for aileron control logic 
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PERIPHERAL ROUTINES 


The routines presented in this section were part of the Apollo / 
Skylab flight programs. They are a baseline for the Space Shuttle 
pointing and maneuvering routines. 
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SECTION 1 


AUTOMATIC TRACKING 


The automatic tracking function applies to any mission situation 
requiring the orbiter to track a designated body (or target vehicle) by 
pointing a specified body fixed vector at the target or to any mission situation 
requiring the orbiter to rotate at a constant rate around a specified axis. 

Typically, automatic tracking can be utilized for rendezvous navigation, local 
vertical tracking, landmark tracking, and passive thermal control. In this mode 
of operation, the autopilot (DAP) is utilized as a rate command system. 

Attitude control, if required, is maintained by the Tracking Attitude Routine (R61) 
with attitude feedback from the navigation functions. RBI operates at a variable 
computation period. At , which is generally much longer than the autopilot 
sampling period of 0. 1 seconds. Initial alignment of the vehicle to the tracking 
attitude is generally performed by R61 using the automatic maneuver routine (RBO). 
In the following TV denotes the target vehicle and ORB denotes the orbiter vehicle. 

1 . 1 Tracking Options 

There are five automatic tracking options as defined below: 

Option Purpose 

0 Point a specified spacecraft vector at the 
TV. Do not constrain rotational attitude 
about this pointing vector. (Use VECPOINT) 

1 Point a specified spacecraft vector at a 
celestial body. Do not constrain attitude 
about this vector (Use VECPOINT) 

2 Rotate at a specified rate about a specified 
spacecraft vector. 

4 Same as Option 0, but constraining attitude 
about the pointing vector to a specified 
value (OMICRON). 

5 Same as Option 1, but constraining attitude about 
pointing vector to a specified value (OMICRON). 
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1 • 2 Option 0: Orient the Spacecraft Such that a Specified Body- Fixed Point- 
ing Vector Aligns with the Estimated Line of Sight to the TV. 

The pointing vector (SCAXIS) is defined by two angles: 

1) A specified positive angle (7) from the spacecraft 

X axis about the spacecraft Z axis (0S7< 360 deg) 

2) A specified positive angle (p ) about the Y'-axis 
produced in the XY plane as a result of y( -90 5 p 5 
90 deg) . 
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The LOS is determined by advancing the estimated ORB and TV state 
vectors to the current time. Once the coincidence of SCAXIS and LOS has been 
established, R61 establishes a spacecraft angular rate to maintain the coin- 
cidence. For Option 0, R61 provides the RCS DAP with the following quantities: 

1 ) desired LOS rate plus LOS correction 

rate in control-axis coordinates 

2) IGA d , MGA d , OGA d DAP reference angles, 

i.e., center of the dead- 
band . 

3 ) ZGA a vector defining the desired incremental 

changes in the IMU gimbal angles every 
0 . 1 second . 

The desired LOS rate (y^) is computed as follows: 

- LOS = ‘ -C 

-correction =<[NBSMl SCAXIS X [REFSMMAT] UNIT <U LOS ))/At C 

where NBSM uses the DAP reference angles THETADX, Y, and Z for the trans- 
formation from navigation-base to stable-member coordinates. 

If I u ,■ I exceeds MAXRATE, distribute as follows: 

I - correction | 


-correction 


= MAXRATE [UNIT (u> 


correction 


-LOS = [REFSMMAT] 


UNIT ( — LOS > V -C ] ' 


- correct ion 


tJ CA = [ NBCA] [SMNB] U 2 LOg , 


where is the LOS rate in stable member coordinates. [ SMNB| is the trans- 

formation from stable member axes to navigation base coordinates. [NBSM) is 
the transformation from navigation base coordinates to stable member coordinates. 

[ REFSMMAT 1 is the transformation from reference inertial coordinates to 
stable member axes. (NBCA] is the matrix for transforming a vector from 
navigation base to control-axis coordinates. j Y w > > and Y q are the position and 

velocity vectors of the TV and ORB in basic-reference coordinates. SCAXIS is 
the specified axis ( 7 , p ) . MAXRATE is the maximum correction rate and is 

a conservative estimate of the R61 computation interval (thus preventing overshoot 
in vehicle attitude by the time of the next pass) . 
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The vector AGA is computed as follows: 


AGA = 


cos IGA sec MG A 
-cos IGA tan MGA 
sin IGA 


0 -sin IGA sec MGA 

1 sin IGA tan MGA 
0 cos IGA 


(0.1 ca 


LOS 


whe re 

IGA = inner gimbal angle 
MGA = middle gimbal angle 
OGA = outer gimbal angle 


where ^q S is the LOS rate in stable-member coordinates computed previously. 


1 - Option 1: Orient the Spacecraft Such that SCAXIS Aligns with the LOS to a 
Specified Celestial Body. 

For Option 1, the LOS is determined by advancing the estimated ORB state 
vector to the current time. Once the coincidence to SCAXIS and LOS has been 
established, R61 establishes a spacecraft angular rate to maintain the coincidence. 
R61 quantities provided to the RCS DAP under tracking-option 1 are as described 
for Option 0 with the following exceptions: 

~LOS = ^ ne ° ^ s i£ht direction, including aberration correction 


and. 


if the target is a star or planet, 
~LOS “correction' 


if the target is the center of the earth. 


“LOS = [ ref SMMAT] 


UNIT (r c )x(v c ) 


£ c 


4 W 

“correction 


1*4 Option 2: Rotate the Spacecraft About SCAXIS 

Option 2 is effected by a special routine known as the Automatic Rotation 
Routine (R67) . SCAXIS is as defined in Option 0. R67 initializes COF , BRATE , 

[ DEL] , [ MIS) , and ADB in KALCMANU and then executes the KALCMANU coding 
to issue the following quantities to the RCS DAP: 
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1 ) u; .... desired rotation rate 

2) IGA d , MGA d , OGA d DAP reference angles, 

i.e., center of deadband. 

3 ) AGA .... a vector defining the desired incremental 

changes in the IMU gimbal every 0.1 seconds. 

These calculations are repeated once per second. 

1 1 5 Option 4: Orient the Spacecraft as in Option 0, but Further Specify an 
Angle (OMICRON) of Rotation About SCAXIS 

SCAXIS is as defined in Option 0; OMICRON is the positive-sense angle 
between planes A and B - where plane A is defined by the negative orbital angular 
momentum vector (-H) and the LOS, and plane B is defined by the spacecraft 
+ Y _ axis and the LOS. 

Alignment of the LOS with either the spacecraft Y axis or the momentum 
vector is proscribed since the necessary two-plane condition for specifying 
OMICRON would not exist. 



PLANE B 
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For Option 4, R61 provides the following quantities to the RCS DAP: 

— ca * * • * desired LOS rate in control-axis coordinates 

2) IGA dj MGA d , OGA d> .... DAP reference angles, i.e., center of the 

deadband . 

3) AGA .... a vector defining the desired incremental changes in the IMU 

gimbal angles every 0.1 second. 

The desired LOS rate (u? rA ) in control-axis coordinates is obtained as 
follows: ' ' 


“LOS r W "-C 

( UNIT (u T ) x (v w - v ) 

= [REFSMMAT ] LOS - W - C 

IjUo i i . 

u 

l-LOSl 

w ca =[nbca] [smnb] u los , 

W ^ ere -LOS * s rate in stable member cooridnates; r^, v^, and 

are the position and velocity vectors of the TV and ORB in basic reference coordi- 
nates . 


The vector (AGA) is computed as follows: 


AGA 


cos IGA sec MGA 

-cos IGA tan MGA 

sin IGA 

+ u , . 

—correction, 


0 -sin IGA sec MGA 

1 sin IGA tan MGA 
0 cos IGA 


(0. 1) u, 


LOS 


where 


m 

correction 



-OGA 



DY 



•IGA 

MGA 


0 . 1 

At^ 


^DX * ® DY J ^ DZ des ^ red gimbal angles ?s 
At c = conservative estimate of R61 computation interval 


The desired gimbal angles ® dy * ^DZ are c° m Puted as follows: 

First, four unit vectors are computed! 

s A = cos (OMICRON) unit (u LOg X [v c X r c ] ) 

+ Sin ( ° MICR ° N) unit [“ L OS X (—LOS X <^C X V)] 
s' A = unit ( SC AXIS X ^ sc ) 
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-B ‘ -LOS 


- 1 B = vSCAXIS 

Given these unit vectors, the unit vectors x, y, and z are calculated defining 
the desired directions of the x ( y, and z axes of the navigation base with respect 
to the Basic Reference Coordinate System. The unit vectors x, and z and the vectors 

-SM' ^SM* -SM from the matrix REFSMMAT are then used to compute the desired 
gimbal angles 0 DX , ^DZ* 

1 . 6 Option 5: Orient the Spacecraft as in Option 1, but Further Specify an Angle 
(OMICRON) of Rotation About SC AXIS . 

SCAXIS is as defined in Option 0; 

OMICRON is as defined in Option 4. 

For Option 5, R6 quantities provided to the RCS DAP are as in Option 4, 
except as follows: 


and. 


-LOS = line of sight direction including aberration correction 


if the target is a star or planet 


-LOS = 0; 


if the target is the center of the earth. 


-LOS ' [REFSMMAT] 


UNIT (r c ) X (y c ) 

iTi 
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SECTION 2 


ON-ORBIT ATTITUDE MANEUVER ROUTINE 


2 . 1 Introduction 

The purpose of the attitude maneuver routine (ATTITUDE_MANEUVER) is to 
automatically reorient the spacecraft during coasting flight. ATTITUDE^ MANEUVER 
is an adaptation of the R60 routine from the Apollo* 5 * flight programs to the Space 
Shuttle with the assumption that a 3-gimbal IMU is used. This implementation is 
coded in HAL and described in accordance with preliminary structured program- 
ming techniques. The general scope of variables and structure of ATTITUDE 

MANEUVER is illustrated in Fig. 2-1. Each of the modules and variables used in 
the routine are defined in Tables 2-1 and 2-2 respectively (Subsection 2.7). 


The desired attitude to which the spacecraft is to be aligned can be specified 
in one of two ways: 

1) By a set of terminal desired gimbal angles, ^ c j. 

2) By two unit vectors, and ^ANB’ 

|#Jis the set of desired (or commanded) gimbal angles defined by 


where 


*c 

Q 


desired outer gimbal angle 

e 

c. 

- 

desired inner gimbal angle 

e 


desired middle gimbal angle 

m 


_ - 


u A T is a body-fixed axis vector, u A , with components in navigation-base- 
coordinates, u pg]yi is a desired direction vector, u^, with components in stable- 
member coordinates and y D g M specifies the direction in inertial space to which 
u is to be aligned . A flag switch, EUAG^3AXIS is provided to specify which 

of the two types of inputs is to be used. 

If the vector type of input is chosen, the routine known as VECPOINT is used 
to compute a set of desired terminal gimbal angles, [# C L based on the values of 
-DSM' anc * -ANB J,anc ^ current g irri fr a l angles, [0]. Regardless of the type of 
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COM POOL USED BY ATTITUDE MANEUVER 



Fig. 2-1 ATTITUDE MANEUVER Struct 












input, KALCMANU is called if the automatic option is selected. Otherwise control 
returns to the executive system. KALCMANU computes a set of parameters (u^, A) 
that determine how the reorientation is to be performed. The maneuver is per- 
formed by rotating the spacecraft about the axis of single equivalent rotation, u^, 
by the angle A^ . The vector u r and the angle A^ are completely determined by 
the initial orientation, [0], and the desired terminal orientation, [0^]. Once and 
A^ are computed, KALCMANU calls the procedure STEERING. 

The procedure, STEERING, cycles every At c seconds and computes, on each 
pass, five sets of quantities to be used by the digital autopilot in following the angu- 
lar path defined by u^, A. The first of these quantities is a set of autopilot refer- 
ence angles, [0^](or intermediate desired CDU angles), which is updated regularly 
during the maneuver. To achieve a smoother sequence of commands between suc- 
cessive updates, the program generates a set of incremental CDU angles, 

[A0^],to be added to [0^] by the digital autopilot. The steering routine also com- 
putes the component maneuver rates, (in spacecraft coordinates), a set of bias 
accelerations, g , and a set of bias angles, The bias angles, are 

to be used by the autopilot to prevent fuel-consuming overshoot when starting and 
stopping an automatic maneuver. 

The preceding method allows the autopilot to establish and maintain a space- 
craft rotation about the vector u at a fixed rotation rate w The maneuvers are 

-r -d 

timed in open-loop fashion so that after a predetermined interval [A0^], ck and 
^BIAS^ are Set zero anc * the autopilot reference angles, [0^], are set equal to the 
terminal angles, [0 ]. Thus, upon completion of. the maneuver, the spacecraft 
will finish up in a limit cycle about [0 ] . 

2 . 2 ATTITUDEJVIANEUVER 

A logic flow diagram of ATTITUDEJVIANEUVER is illustrated in Fig. 2-2. 

If the vector option is selected by the calling program, the routine VECPOINT is 
used to compute the terminal desired gimbal angles [0 ]from an< ^ * 

If the automatic maneuver is to proceed the automatic mode switch, AUTOJVIODE, 
must be set "ON". If it is ”ON lf , KALCMANU is called to calculate the rotation 
vector, u^, and the angle of rotation, A^, and to initiate the cycling of STEERING. 

2.3 VECPOINT 

The desired gimbal angles, [ 0 ] # computed by VECPOINT are determined on 
the basis of minimizing the rotation angle required to align with in 

general, this procedure will reduce the fuel required for the maneuver by not 
constraining the orientation of the spacecraft about the pointing axis, defined by 
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-DSIVr Additionall y» the attitude about U D g M will be left unchanged. Consequently, 
the crew may specify the attitude about the pointing axis to satisfy additional com- 
munication, visibility, or thermal constraints, if desired. This is done by manually 
reorienting the spacecraft about the pointing axis either prior to or after the maneu- 
ver. If this procedure is done after the automatic maneuver, a final trim will usually 
be required to assure that is aligned with u D gjyp 

To compute the desired gimbal angles, [0j, from y DgM , u^ NB> and the cur- 
rent gimbal angles, [ 8 ] , several service routines are required. These routines are 
described below. 


2.3.1 ANGLES_TO_DCM 

This procedure converts three gimbal angles, [0],to a direction-cosine matrix, 
C, that relates the corresponding spacecraft navigation-base orientation to the 
stable-member frame. The formulas for this conversion are: 




c i,i 

c 

1,2 

C 

1,3 

* 







c 

= 

C 2, 

1 

c 

2, 2 

c 

2, 3 



c 

_ 3, 1 

c 

.3, 2 

c 

3, 3 

1, 1 

_ 

cos 

8 . cos 8 





i 

m 



1, 2 

- 

-cos 8. 

sin 8 

cos 8 



i 

m 


o 

1, 3 

- 

cos 

8. sin 8 sin 

8 + 



l 

m 


o 

2, 1 

- 

sin 

8 






m 




'2, 2 

= 

cos 

8 

cos 8 





m 

0 



2, 3 

= 

-cos 8 

sin 8 





m o 


3, 1 

= 

-sin 0 . 

cos 8 





l 

m 



3, 2 


sin 

8. sin 8 cos 

8 + 



l 

m 


0 

3, 3 

- 

-sin 6 . 

sin 8 

sin 

i 8 



l 

m 


o 


2-5 


where 


0 - outer gimbal angle 

0. = inner gimbal angle 

- middle gimbal angle 

Also see Fig. 2-3 for the logic flow diagram of ANGLES_TO DCM . 

The interpretation of this matrix is as follows; If a , a , a represent the 

components of a vector in navigation-base axes, then the components of the same 

vector in stable-member axes (b , b . b ) are 

x* y* z 


( b 


a 

1 x 


x 

{ b 

= C(9) 

i j 

a 

j y 

y 

1 b 


a 

1 y 

I 

i z 


2.3.2 DCM_TO_CDU ANGLES 

This procedure extracts the CDU angles from a direction-cosine matrix, 
that relates the navigation-base coordinates to the stable-member coordinates. 
The formulas for this conversion are 


arcsin C 


2 , 1 


0. - arcsin (-C„ - /cos 0 ) 

l 3,1 m 

[ If C l, 1 is ne g ative > is replaced by ir SGN (0 ) - 9. 


■.] 


arcsin (-C„ „/cos 8 ) 
2, 3 m 


[l f C 2 2 is negative, Q q is replaced by v SGN (0 q ) - 
The logic flow diagram is illustrated in Fig. 2-4. 


2.3.3 ROTATION_MATRIX 

This procedure computes the direction-cosine matrix, R, that relates one co 
ordinate frame to another that is rotated with respect to the first by an angle A 
about a unit vector u . The formula for this matrix is 

# # y ^ 

R A) = I cos A + u r u^ (1-cos A) + UX r sin A 
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C 1 ^ = COS (0 2 ) COS (0 3 ) 



Fig. 2-3 ANGLES TO DCM : Procedure ([#]) Assign (C) 
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Fig. 2-4 DCM_TO_ANGLES : Procedure (C) Assign ([#]) 
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where 


i 


1 0 0 
0 1 0 
0 0 1 



2 

u 

rx 

u 

rx 

u 

ry 

u u 
rx rz 

T _ 

u u 
-r - r 

u u 
ry rx 

2 

u 

ry 

u u 
ry rz 


u u 
rz rx 

U rz 

u 

ry 

2 

u 

rz 


0 

-u 

rx 

u 

ry 


* 

UX = 
r 

U rz 

0 

-u 

rx 



-u 

ry 

u 

rx 

0 



u = unit rotation vector resolved into spacecraft 
-r 

axes 

A = rotation angle 


The logic flow diagram is illustrated in Fig. 2-5. 

The interpretation of R is as follows: If a x< a , a z represent the components 

of a vector in the rotated frame, then the components of the same vector in the 

original frame (b , b , b ) are 
x y z, 


1 b \ I 

a 

1 x I I 

X 

1 I >Jc 


| b ) = R (u , A) < 

a 

) y ( “ r 

y 

(hi 

a 

\ z ) ' 

[ z 


The process used by VECPOINT in determining the terminal desired angles, 
[0 ], is illustrated schematically in Fig. 2-6 and in logic flow form in Fig. 2-7. 
VECPOINT uses two frames of reference: the L frame and the M frame - as well 
as the stable-member (SM) frame. The L frame is simply a base frame used as 
the starting point for subsequent rotations. In VECPOINT computations, the L 
frame is the initial navigation-base frame computed from the initial gimbal angles. 
The vector is resolved from the stable-member frame into this frame as 
U DL us * n £ the t ransformation C^SM The ^ f rame is then rotated so that 
and -DL are coin cident. The r °tation vector, u p , is defined by the cross product 
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Fig. 2-5 ROT AT ION_M A T RIX : 

Procedure (u r ,A ) Assign (R) 


I R„, <u . A) 





Fig. 2-6 The Coordinate Frames and 

Transformations used in VECPOINT 




(c) The SM frame 










between y AL and y DL . Rotation about the cross-product vector will assure that 
spacecraft maneuvering will be minimized and that the rotation vector will have no 
component in the or u DI direction. 

If the two vectors, y AL and y DL , are either parallel or antiparallel, a special 
computation of u^ is performed. When y A ^ and u^^ are almost parallel then [ 6 J 
is set equal to [#^]. If the two vectors are almost antiparallel, then a 180° rotation 
is required for alignment. In this case, the rotation vector is chosen so that u 
is perpendicular to y AL and lies in the plane that contains u ysM(L) and u yL . If, 
however, y A ^ happens to be perpendicular to the plane, then any vector in the plane 
may be used for In this event the x-axis is chosen. From y r and A the rotation 

matrix, between the M- and L- frames is computed. Multiplying R by 

^LSM ^ e ^ errnines the direction - cosine matrix ^ JVIS IV1 ^ rom w ^ich the terminal desired 
gimbal angles, [e c ]. are extracted. 

2 . 4 KALCMANU Maneuver- Parameters Calculation 

The procedure KALCMANU is designed to generate commands for the auto- 
pilot to reorient the spacecraft from an initial attitude to some desired attitude 
specified by [<g. Although a number of authors have discussed optimal solutions 
to this problem, the task is a formidable one because of the nonlinear nature of 
rigid- body dynamics and the kinematics of rotational motion. Optimal solutions 
for a particular class of spacecraft configurations are treated in Reference 1. In 
this reference, the criterion for optimality is based upon the product of maneuver 
duration and fuel expenditure. Using the same criterion, previous design studies 
based upon a number of suboptimal approaches* 2 ^ have shown that a reasonable 
design compromise is to rotate the spacecraft about a single inertially fixed axis. 

This technique is employed in KALCMANU. 

The process used by KALCMANU in computing the rotation axis, u , and the 

rotation angle, A^, is illustrated in Fig. 2-8. From the current gimbal angles 

KALCMANU computes the transformation matrix between the initial spacecraft 

axes and the stable-member axes, This procedure also computes the matrix 

relating the desired terminal spacecraft attitude and the stable-member axes, 

^CSM ’ The transformation between the desired and the initial attitude is R , where 
. CB J 

where 


* 5jC rj. # 

R CB ■ C BSM C CSM 

* 

■^CB can partitioned into its symmetric and antisymmetric components as follows. 
For the symmetric part. 


❖ 



1 ,* 

2 (R CB * 


# t 

r cb> 
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By comparing this with the equations described in ROTATION_MATRIX it can *be 


seen that 


>}t ije rp 

R 0 = I cos A^ + (1-cos A m > 


M 


For the antisymmetric part. 


* 1 * 

R A = ■ 2 )(R CB 


# '"T 

- r cb ) 


-u 
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rz 
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rz 
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r y 
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rx | 
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sin A 
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r y 


The maneuver angle can easily be obtained from Rq B as 


/ Trace(R CB )-l j 


A M = cos 

The rotation vector, u r> is normally extracted from the antisymmetric part 


of R that is, 


unit (- R , 


> Ry 


-R , 


3, 3 


1,3 


1,2 


As sin A^ approaches zero, however, the unit operation becomes indeterminate 
and an alternative calculation must be performed. If the maneuver angle is less than 
0.25 deg, then no complicated automatic maneuvering is necessary and therefore 
one can simply set the autopilot reference angle, [0^], equal to the terminal desired 
angles [ 0 ] and exit. This is done by calling HOLD_ATTITUDE . If A^ >170 deg, 
a method employing the symmetric part of R^ B is used to determine u f , as illus- 
trated in Fig. 2-8. The magnitudes of the components of u r are easily determined 
as 



To obtain the signs of each component, it is first necessary to determine the 

* 

sign of the largest component by referring to the values of R^. This is done to 
avoid the problem caused by truncation as sin A^ approaches zero. If A^ is equal 
to 180 deg, then the sign of the largest component may be arbitrarily chosen. The 

signs of the remaining components may be determined by examining the off-diagonal 

* 

terms in Rg. 
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2. 5 Steering Through a Maneuver 

The KALCMANU parameter calculations described in the preceding subsection 
solve the problem of what angular path to take in order to reorient the spacecraft 
without excessive fuel expenditure. However, the problem of how to steer the space- 
craft along this path must now be considered. One solution is to iteratively solve for 
u r and A m explicitly as a function of the current spacecraft attitude and to control 
the spacecraft angular velocity, by firing the RCS jets so that is parallel to u f . 
This approach, however, leads to several problems. First of all, because of the 
deadband requirements of the autopilot, and control inaccuracies, explicit guidance 
can lead to excessive changes in control commands as the spacecraft approaches the 
terminal attitude. This problem of terminal guidance is characteristic of explicit 
guidance schemes unless special consideration is given to the terminal-control 
problem. Secondly, there is a timing problem since the computation of and A^ 
is somewhat lengthy and would require an excessive amount of computer time. 

Instead, it is better to formulate the problem in terms of perturbations from a smooth 
reference trajectory that satisfies the end-point constraints. The method employed 
is to solve the global problem once via KALCMANU in order to establish a reference 
trajectory (defined by A^) and then have the autopilot solve the micro-problem 
of controlling the spacecraft deviations from this reference trajectory, using small- 
angle approximations. This approach not only avoids the terminal-guidance problem 
and the timing problem but also provides "tighter” path control for satisfying other 
boundary conditions that may be placed upon the selection of the reference trajectory. 

In order to determine the required steering interface between KALCMANU 
and the autopilot, it is first necessary to define the nature of the autopilot control 
variables. To simplify the discussion, the problem of rotation about a single axis 
is considered first. There are two state variables of interest; namely, the attitude 
error, 0 defined by 

0 = 6 - 9 , 

e d 

where 

0 = actual vehicle attitude 

6 , = desired vehicle attitude (constant) 
d 

and the angular velocity of the vehicle, w. The "state" equations are 

0 e (t) = « (t) 

and 

w(t) = f (0 , w) 
e 
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where f (6 & , to) is the control acceleration (applied by the jets) as determined by non- 
linear phase-plane switching logic. The control law is defined so that 6 and to will 

e 

be nulled to within the control-law deadbands so as to maintain attitude hold about 

the reference attitude, 6,. 

d 

New state variables can be defined as follows: 

V*> * 8,, » - V‘> * <w 

and 

“ (t) = w(t) - 0) 
e d 

where 

*d“> ■ *<,'*„> + <“d> ft -*o> 

= a constant 

^BIAS = a cons t an t 
^d^o^ = B fto> = constant 

The state equations for the new variables are 


< f ) - OJ 

e e 

w = f { 0 , oj ) 
e e e 

With these control variables, the autopilot will attempt to null the rate error w 

* g* 

thus forcing the spacecraft to rotate at a constant rate. Note that at the beginning 
of the maneuver the desired angle 0 d (t Q ) is set equal to the initial angle 0(t ) . The 
constant, added to the attitude error to prevent overshoot when starting and 

stopping the maneuver, 0 BI is computed as 


where 


jS 


BIAS 


J d hi 1 

-rnrr 


= the control acceleration 

To interpret the effect of refer to Fig, 2~9a, This figure shows a 

simplified phase-plane logic and a typical phase-point trajectory. The maneuver 
starts at point A with # e = 0 and w = 0. With the application of steering commands, 
the phase-point instantaneously jumps to a biased condition at B where <t> = 0^ 

andc^ = -w d {to d is negative for the case illustrated). The phase-point drifts to 
point C, where the control jets are turned on to reduce the rate error toward zero 
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at point D. The trajectory then traverses a limit cycle DEFG until the maneuver 
finishes at point H, where the biases $ BIAS and “ d are removed. At this point, the 
trajectory jumps to an unbiased condition at point I. Note that the actual rate is 
approximately equal to the desired rate and the phase-point drifts to point J, where 
the jets decelerate the vehicle back into a limit cycle around the desired attitude. 
The situation without angular bias, & mAS > is illustrated in Fig. 2-9b. Note the 
fuel-consuming overshoot in this case. The inclusion of the bias term is based on 
the assumption that the initial rates are small and that the desired rates are 
achievable during the maneuver. For the three-dimensional problem, it must also 
be assumed that 3 T AC is small enough so that the small-angle assumptions are not 

Dino 

violated. 

A bias acceleration, or , is also computed and sent to the autopilot to better 
approximate the switching curves in the phase-plane design. In the formulation of 
the phase-plane design a simplification is made to the equation of motion 

• * _ i # _ i * 

w(t) = J M-J (wx Jw) 

where 

M = the torque supplied by the jets 

J = a , the control acceleration 

— - c 

i * 

-J (o> x JoA = a , the acceleration due to the 
- — -m 

non-diagonal moment of 
inertia J. 

The bias acceleration, ot is assumed small and neglected to arrive at the state 
equations 


— — 
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— - 
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m - 0 d <t) 
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0 1 
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~ e 
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ot 

w<t) - “ d (t) 
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and allow the design of three independent phase-planes. To refine the determination 
of the switching curves, is sent to the autopilot and is used in estimating the net 
acceleration . 

The above steering logic has been implemented in the procedure STEERING 
which is described below. 

2.6 STEERING 

The procedure STEERING computes five sets of steering parameters to be 
sent to the autopilot. These parameters are 
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Fig* 2-9 Maneuver Phase-Plane Trajectories 


2-18 


1 ) 


[0^], a set of intermediate desired gimbal angles 

2) [A 0^], a set of incremental desired gimbal angles 

3) a Set Si^kal an g* es 

^ the desired angular velocity vector 

5) a , a set of bias accelerations 

-m 

STEERING operates on a cyclic basis with period At c . The parameters[0gj 
and are computed once on the first pass whereas [0^]and [A0^]are computed 
every At seconds until the maneuver is complete. 

A logic flow diagram is illustrated^ in Fig. 2-10. The initial gimbal angles 
are read and the FIRST_PASS switch is set "ON". The incremental angle of rotation 
is computed from the requested magnitude of angular velocity and the time increment 
At . From this incremental rotation angle and the rotation vector u^, an incre- 
mental rotation matrix AR is determined by calling ROTATlON_MATRIX. Two timing 
parameters are computed. These are the duration of the maneuver, t^. 


M 


= A-./co 
M r 


and the time at which to stop the maneuver, t , 


t = t, ff + t , 
s M cycle 

where t c ^ c ^ e has been set equal to the current time. 

Cyclic computation of the steering parameters takes place in open-loop fashion 
until the computed time to stop the maneuver, t , is reached (Fig. 2-10). If the 
AUTO_MODE switch is "OFF" then STOPRATE is called which sets the desired 
reference gimbal angles [0^] to the present gimbal angles and nulls the remaining 
four sets of parameters [A0^], [0-gj^g], and as illustrated in Fig. 2-11. 
Control then returns to the executive system. However, if the AUTO__MODE switch 
is "ON" the direction-cosine matrix between the navigation-base and stable-member 
coordinate system^ is determined by 

* * * 

CbSM N+1 = Cbsm n AR 

From this matrix, the gimbal angles projected to the time of the next pass, 
are calculated. 

On the first pass STEERING sets the FIRST_PASS flag to "OFF" and com* 
putes « [0 ], and a as illustrated in Fig. 2-10. The desired angular velocity 

Q DinO m 

vector, is computed from the requested magnitude of angular velocity, 0J r > and the 
rotation vector, u^, supplied by KALCMANU. 
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* 



Fig. 2-11 STOPRATE : Procedure 
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8 BIAS(U) 


“dc(U) 1 “dc(U) 1 / (2 1 ff c (U) 1 } 

S BIAS(V) 

— 

“dc(V) 1 “dc(V) 1 / <2 ,<a, c (V) 1 } 

0 BIAS(W) 


“dc(W) l u dc(W)l / (2l «c (W)l) 


where are the control accelerations as determined by the auto- 

pilot . 


On every pass STEERING also updates the incremental desired gimbal angles, 
[A0 d ], by 


[A» d ] - <[<g - [ej) (t fast /At c ) 

[0^] = desired reference gimbal angles projected 
At c seconds in advance 

[fl B ] - present desired reference gimbal angles 
At c = period of the STEERING cycle 

t fast = sam P lin g period of the autopilot 


These incremental angles are subsequently added to [0 rf ] by the autopilot at the start 
of each autopilot pass. [9j]is then updated to [flgjand [0 R ]is replaced by [0 ]for 
use on the subsequent pass. The time of the next proposed pass, t le< is computed 
by adding At c to the present value of t cycle> If t cycle is less than the termination 
time of the maneuver, t g , another pass is performed. Otherwise, HOLD ATTITUDE 
is called to update the desired reference gimbal angles, [0 ], to the commanded 
terminal set of gimbal angles [0^Jand nulls the remaining steering parameters 
[A0 d ], [/3 mAS ], <f d , and <* m * This updating is illustrated in Fig. 2-12. Control is 
then returned to ATTITUDE__MANEUVER which in turn transfers control to the 
executive system. 
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Fig. 2-12 HOLD ATTITUDE : Procedure 
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Subsection 2.7 Table 2-1 ATTITUDE MANEUVER MODULE DICTIONARY 
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COMPOOI. heing on level 0 
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Table 2-2 ON-ORBIT ATTITUDE MANEUVER VARIABLE-CONSTANT DICTIONARY (Cont) 
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Table 2-2 ON-ORBIT ATTITUDE MANEUVER VARIABLE-CONSTANT DICTIONARY (Cont) 
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Table 2-2 ON-ORBIT ATTITUDE MANEUVER VARIABLE-CONSTANT DICTIONARY (Cont) 
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